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CHAPTER  1 


INTRODUCTION 


1.1  GENERAL 

Located  at  Holloman  Air  Force  Base,  ten  miles 
southwest  of  Alamogordo,  New  Mexico,  is  the  United  States 
Air  Force's  6585th  Test  Group,  whose  prime  responsibility 
is  the  testing  and  evaluation  of  new  weapon  systems 
prior  to  their  being  introduced  into  the  Air  Force 
inventory.  These  systems  to  be  tested  range  from 
guidance  systems  and  inertial  components,  to  blast 
interaction  tests  on  various  air  frames,  to  aircraft  seat 
ejection  tests.  This  testing  is  accomplished  on  a  large 
number  of  test  fixtures,  among  which  a  few  include  a 
precison  high  acceleration  (g)  centrifuge,  precision  rate 
tables,  vibration  tables,  three-axis  flight  simulators, 
and  the  high  speed  test  track. 

1.2  DESCRIPTION  OF  THE  HIGH  SPEED  TEST  TRACK 

The  high  speed  test  track  is  a  facility  where 
captive  land-based  vehicles  may  be  tested  to  speeds  up  to 
9,000  feet  per  second.  These  captive  vehicles  are  called 
sleds  and  may  be  of  the  configurations  of  dual  rail, 
monorail,  wide  gage,  narrow  gage,  or  of  an  outrigger 
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form.  See  Figure  1.1  for  a  pictorial  representation  of 
these  sled  types . 

The  track  facility  itself  consists  of  50,788  feet  of 
continuously  welded  track.  There  are  two  rails  spaced  84 
inches  between  centerlines  (dual  rail  spacing)  running 
the  entire  50,788  feet.  The  last  15,200  feet  contains  a 
third  rail  which  is  26.32  inches  east  of  the  centerline 
of  the  east  rail  of  the  dual  rail  section  (narrow  gage 
spacing)  . 

The  rail  is  a  heavy  duty  crane  rail  of  171  pounds 

per  linear  yard  versus  65  pounds  per  yard  for  standard 

railroad  rail,  and  comes  in  39  foot  lengths.  The  rail  is 

tied  down  in  both  the  lateral  and  vertical  directions 

every  52  inches.  The  rail  is  stretched  before  being 

secured  by  tie  downs  such  that  a  tensile  load  will  be 

maintained  in  the  rail  at  temperatures  up  to  120  degrees 

Faherenheit.  The  master,  or  west,  rail  is  aligned  to  a 

reference  line  to  within  +_0.005  inches,  with  the  other 

two  rails  aligned  to  within  +_0.010  inches  of  the  true 

gage  spacing  from  the  master  rail.  The  track  does  not 

follow  earth  curvature,  but  consists  of  straight  segments 

connected  by  transition  arcs  of  one  million  foot  radii. 

See  Figure  1.2  for  a  representation  of  the  track 
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details. 


TRACK  DETATI S 


1.3  BACKGROUND  ON  SLED  DESIGN 


Present  day  track  test  programs  are  requiring  sleds 
which  operate  in  the  hypersonic  velocity  range  for 
monorail  sleds  and  up  to  Mach  3.0  for  dual  rail  sleds. 
The  design  of  these  sleds  demands  load  prediction 
techniques  that  insure  sled  integrity  while  minimizing 
sled  weight. 

Design  loads  are  based  first  on  the  "quasi-steady- 
state"  forces.  The  qua s  i -s  t  e  a d y- s t ate  forces  are  slowly 
varying  forces  and  are  defined  as  the  forces  acting  on 
the  sled  due  to  inertia,  thrust,  friction,  braking,  and 
aerodynamics.  Inertia,  thrust,  braking,  and  the  friction 
forces  are  adequately  estimated  based  on  weight  analysis, 
momentum  exchange  equations,  and  past  experimental  tests. 
The  aerodynamic  forces  are  basically  derived  from  wind 
tunnel  testing.  However,  the  results  of  the  wind  tunnel 
tests  are  used  with  some  reservation  due  to  the 
approximations  required  to  implement  the  wind  tunnel 
tests,  i.e.,  scaling  of  sleds,  boundary  layer  effects  of 
the  ground  plane-sled  interaction  surface,  and  the  fact 
that  the  results  are  only  good  in  the  supersonic  velocity 
regime.  The  equations  of  motion  are  solved  to  estimate 
the  q  u  a  s  i - s t e a d y - s t a t e  forces  where  tables  of  sled 
weight,  thrust,  and  moment  arms  versus  time,  as  well  as 
aerodynamic  forces  versus  velocity,  are  interpolated  by 


means  of  a  computer. 

Besides  the  q u a s  i - s t e a d y - s t a t e  forces,  the  sled 
structure  is  also  subjected  to  significant  oscillatory 
forces  caused  by  surface  roughness  of  the  rails, 
propulsion  transients,  aerodynamic  buffeting,  acoustics, 
and  sled-track  interactions.  These  "dynamic”  forces  have 
been  handled  in  various  ways  over  the  years,  none  of 
which  have  been  totally  successful  in  accurately 
predicting  sled  dynamic  forces.  These  methods  have  been 
successful  from  the  standpoint  they  yielded  a 
conservative  sled  design.  The  history  of  these  sled 
dynamic  estimations  will  now  be  briefly  discussed. 

1.3.1  LAMBDA  METHOD3 

This  method,  developed  in  early  track  testing, 
amounts  to  adding  an  additional  steady-state  force  on  the 
sled.  This  additional  force  was  defined  in  terms  of  a 
factor  called  lambda  (X  ),  which  was  based  on  sled 
velocity.  This  factor  was  applied  at  the  sled’s  center 
of  mass  and  acted  as  an  inertial  loading  in  the  vertical 
and  lateral  directions  on  the  sled.  Figure  1.3  shows  the 
lambda  factors  used  in  previous  dual  rail  and  monorail 
sled  designs. 
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LA*  DA  (A)  FACTOR 


Figure  1.3  LAMBDA  FACTOR  VERSUS  VELOCITY 
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1.3.2  BETA  METHOD 


Wind  tunnel  tests  were  conducted  on  dual  rail  sleds 
to  measure  aerodynamic  loads  (lift  and  drag)  on  the 
slipper  beams.  These  were  added  to  the  computer- 
predicted  thrust  and  acceleration  forces  to  obtain  the 
q u a s i - s t e a d y - s t a t e  slipper  loads.  These  results  were 
subtracted  from  actual  sled  run  data  and  left  only  the 
dynamic  forces.  A  beta  (8  )  proportionality  factor  was 
derived  which  related  the  dynamic  forces  to  the 
predicted  qua s i -s t e ad y-s t a t e  forces.  The  vertical  forces 
were  then  defined  as  the  product  of  the  beta  factor  times 
the  predicted  qua s i -st e ad y-st a t e  forces. 

1.3.3  MONORAIL  SLED  DESIGN  TECHNIQUES 

Neither  the  lambda  no-  beta  method  are  acceptably 
accurate  for  use  on  monorail  sled  designs.  This  is 
mainly  due  to  the  fact  that  oscillatory  forces  in 
monorail  sleds  are,  in  general,  an  order  of  magnitude 
greater  than  in  dual  rail  sleds.  Also,  due  to  the  very 
high  stiffnesses,  the  vibrational  amplitudes  of  the 
monorail  sleds  are  a  function  of  sled  velocity  as  well  as 
the  sled  structural  response,  whereas  the  dual  rail  sled 
vibration  strain  energies  are  located  at  the  lowest 
natural  modes  of  vibration  (in  the  order  of  15-20  Hz.). 


Many  theoretical  studies  have  been  performed  trying 


to  correlate  the  predicted  vibration  results  of  monorail 

sleds  to  the  measured  results.  These  have  ranged  from 

subjecting  a  mathematical  model  of  the  sled  of  varying 

complexity  to  a  nonlinear  forcing  function  of  one  form  or 
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another.  Mixon  provides  an  excellent  synopsis  of  all 
the  theoretical  attempts,  which  in  the  final  analysis 
yielded  only  marginal  results. 

Mixon^  then  went  through  an  extensive  design  and 
test  program  to  validate  force  measuring  transducers 
(Figure  1.H)  for  use  with  primarily  monorail  sleds.  All 
forces  were  resolved  into  their  respective  components 
within  three  percent  of  the  applied  loading.  By  the 
proper  addition  or  subtraction  of  components,  total  force 
in  a  given  direction  could  be  defined.  When  multiplied 
by  the  proper  lever  arm,  the  appropriate  moments  acting 
on  the  sled  and  slipper  assemblies  are  determined. 

Mixon  then  mathematically  predicted  the  resulting 
dynamic  forces  due  to  a  sled  traversing  a  "rough"  rail. 
He  did  this  by  developing  a  nonlinear  model  of  a  monorail 
sled.  This  model  consisted  of  a  rigid  body  approximation 
of  the  vehicle  with  up  to  six  (6)  degrees  of  freedom 
associated  with  the  center  of  mass.  This  resulted  in  a 
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sled  model  similar  to  that  shown  in  Figure  1.5. 


FIGURE  1.5  RIGID  BODY  MATH  MODEL 


The  supports  for  the  sled,  represented  as 
spring/damper  assemblies,  mathematically  allow  for  the 
clearance  that  exists  between  the  slipper  and  the  rail. 
The  slipper  gap  introduces  the  nonlinearity  to  the 
problem — the  classic  case  of  a  nonlinear  spring  system, 
one  with  clearances  (Figure  1.6). 

The  problem  is  attacked  by  solving  the  equations  of 
motion  about  the  sled  center  of  mass  to  find  the  slipper 
deflections.  If  the  resultant  slipper  deflections  exceed 
the  specified  slipper  gap,  »"ai  1  —  slipper  contact  occurs. 


This  contact  induces  a  dynamic  force  into  the  sled 


through 
a  sled 
vehicle 


the  slipper.  By  iterating 
run,  the  maximum  dynamic 
can  be  predicted. 


on  this  concept  during 
forces  acting  on  the 


FIGURE  1.6  CLASSICAL  EXAMPLE  OF  SPRING  WITH  GAP 

This  approach  is  totally  contingent  upon  an  accurate 
forcing  function  (rail  roughness)  being  applied  to  the 
model.  Original  attempts  at  defining  rail  roughness  were 
based  on  a  400  foot  section  of  track  that  was  surveyed  on 
the  top  surface  only  at  ten  inch  intervals. 

Correlation  by  Mixon  showed  that  the  rigid  body 
approach  with  a  nonlinear  forcing  function  and  a  five  (5) 
percent  structural  damping  assumption  nverpredicted  the 
vertical  dynamic  forces  by  20  to  50  percent  of  measured. 
However,  lateral  dynamic  forces  were  consistantly  under¬ 
estimated  by  50  to  60  percent.  The  lateral  force 
discrepancies  were  assumed  to  be  attributable  to  the 
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assumptions  that  (a)  the  lateral  rail  forcing  function 
was  identical  to  the  vertical  rail  survey,  (b)  to 
neglecting  any  cross-track  wind  buffeting  effects,  or  the 
lateral  rail  roughness,  and  (c)  to  misalignment  of  the 
center  of  gravity  in  the  lateral  direction. 

1.3.4  SLEDYNE--PHASE  I7 

Based  on  the  results  and  recommendations  of  Mixon’s 
monorail  work,  a  contract  was  let  to  Mechanics  Research 
Incorporated  (MRI)  to  develop  a  nonlinear  math  model  for 
use  with  dual  rail  sleds,  as  well  as  monorails.  What 
resulted  was  a  design  guide  and  the  development  of  the 
SLEDYNE  computer  program  for  predicting  dynamic  force 
resultants  on  the  slippers  of  a  traveling  sled. 

The  development  of  SLEDYNE  was  based  on  the  earlier 
rigid  body  monorail  work  of  Mixon.  Instead  of 
restricting  their  solution  of  the  equations  of  motion  to 
rigid  body  response,  the  equations  allowed  for  some  of 
the  impact  energy  to  be  absorbed  by  the  sled’s 
vibrational  characteristics.  The  limiting  factor  in  the 
derivation  was  that  only  pitch  plane  motions  (vertical 
translations  and  pitch  rotations)  were  allowed.  This  did 
not  allow  for  any  lateral  ( c r o s s - t r a c k )  inputs.  These 
cross-track  dynamic  forces  were  considered  to  be  no 
greater  than  the  vertical  forces  to  be  estimated,  and 
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could  be  accounted  for  by  the  superposition  of  a  properly 
scaled  vertical  force. 

The  one  overriding  assumption  inherent  in  this 
derivation  was  that  it  be  possible  to  isolate  two  sled 
properties  and  their  respective  responses.  These  were 
the  slipper  stiffnesses,  which  transmit  only  vertical 
dynamic  forces  into  the  remainder  of  the  sled  structure, 
and  the  sled  body,  which  can  be  defined  by  a  set  of 
normal  vibrational  modes.  By  the  proper  use  of  this 
assumption  all  of  the  sled's  behavior  is  preserved. 

The  approach  then  taken  was  to  iterate  the  equations 
of  motion  for  the  given  sled  trajectory  over  the 
prescribed  forcing  function.  Again  the  forcing  function 
was  the  original  400-foot  section  of  surveyed  rail. 

The  results  of  SLEDYNE  compared  reasonably  well  with 
Mixon's  monorail  measurements.  The  actual  forces  were 
dependent  on  the  number  of  modes  and  damping  percentages 
used.  Five  percent  damping  and  use  of  six  modes 
predicted  vertical  slipper  forces  within  _+ 1  0  percent  for 
maximum  velocity.  No  comparison  could  be  made  for 
lateral  dynamics.  The  results  also  correlated  within  +10 
percent  of  the  previously  measured  vertical  dynamic 
forces  for  dual  rail  sleds.  The  limitations  with  these 
results  were  that  they  were  only  considered  accurate  when 
predicting  forces  during  the  coast  phase  of  the  profile. 
The  simulation  was  never  attmpted  when  a  thrusting  vector 


was  applied  to  the  model  due  to  engine  transients. 

One  important  development  for  dual  rail  sleds  came 
from  this  study.  A  parameter  called  the  SIMP  (Simulated 
IMpact  Parameter)  factor  was  fortuitously  found.  The 
forward  and  aft  slipper  beams,  for  conventional  dual  rail 
sleds,  are  considerably  more  flexible  than  the  sled  body 
so  the  structure  could  be  represented  by  two  spring-ma3S 
systems.  A  simple  equation  was  developed  which  gave  the 
SIMP  factor.  This  factor  could  be  used  to  enter  a  chart 
and  determine  a  resultant  slipper  beam  force  based  on 
actual  measurements  of  the  dynamic  slipper  beam  force.  A 
distributed  inertial  loading  could  be  applied  to  the 
sleds  from  a  proper  combination  of  these  SIMP  factor  beam 
loadings.  This  yielded  an  intermediate  design  stage  that 
allowed  for  model  refinement  before  proceeding  to  the 
SLEDYNE  design  phase. 

1.3.5  SLEDYNE — PHASE  II8 

! 

Several  years  after  the  MRI  development,  the  Air 
Force  Flight  Dynamics  Lab  (FIBRA)  at  Wr ight-Patter son  Air 
Force  Base  was  contracted  to  make  improvements  in  the  use 
of  the  SLEDYNE  program.  This  resulted  in  adapting  the 
present  structural  response  program,  NASTRAN,  for  use 
with  SLEDYNE.  SLEDYNE  was  originally  developed  to  be 


compatible  with  MRI’s  STARDYNE 


In  addition  to  program  adaptation,  FIBRA  did  an 
extensive  study  of  the  forcing  function  being  used.  In 
past  analyses,  the  rail  roughness  sample  set  was  looped 
end  to  end  to  generate  sufficient  rail  length  for  the 
forcing  function.  FIBRA  studied  generating  a  rail  from 
the  given  data  set  but  randomly  defining  the  sequence  of 
the  data  points.  FIBRA  considered  using  a  Monte  Carlo 
method  of  generation,  with  each  slipper  being  independent 
of  the  others,  i.e.,  for  dual  rail  sleds  four  rail 
profiles  were  generated.  Track  analysts  did  not  accept 
this  approach  as  giving  reasonable  results,  so  currently 
the  rail  forcing  function  is  generated  by  means  of  a 
random  generating  process  to  define  the  start  point  and 
then  looping  the  data  end  to  end. 

1.4  LIMITATIONS  IN  DYNAMIC  ANALYSIS 

Present  limitations  in  the  use  of  SLEDYNE  are  due  to 
the  derivation  of  the  equations  only  allowing  for  pitch 
plane  reactions  to  determine  vertical  dynamics.  However, 
the  roll  moment  of  the  sled  also  inputs  forces  in  the 
vertical  direction.  Of  prime  importance  in  the  final 
design  is  the  present  assumption  that  the  lateral  rail 
dynamics  are  a  proportion  of  the  vertical  dynamic  forces, 
which  has  been  verified  only  for  monorail  sleds. 

Another  limitation  of  the  program  is  the  definition 


of  the  rail  forcing  function  that  is  to  be  subjected  to 

the  sled.  This  particular  problem  has  been  addressed. 

q 

Mixon,  Evans,  and  Gilliam7  had  additional  rail  roughness 
measurements  taken  and  then  perfomed  the  statistical 
analysis.  This  survey  was  accomplished  by  taking 
fifty-one  (51)  measurements,  one  foot  apart,  in  a  2,000 
foot  section  of  track.  This  was  repeated  each  2,000  feet 
for  the  entire  50,000  feet  of  track.  The  survey  included 
roughness  measurements  on  the  rail's  top,  side,  and 
bottom  of  the  west  rail. 

Results  of  the  analysis  indicated  the  rail  could  be 
grouped  into  four  (4)  distinct  statistical  roughness 
categories.  Further  analysis  showed  some  surprising 
results:  each  39  foot  section  of  rail  showed  itself  in  a 
power  spectral  analysis,  as  also  did  the  52  inch  tie  down 
spacing.  The  final  conclusion  was  that  not  enough  data 
points  were  available  to  accurately  define  a  frequency 
dependent  rail  forcing  function  for  dual  rail  sleds. 
These  additional  measurements  have  been  made  and 
analyzed  for  the  final  definition  of  the  rail  forcing 
function.  These  are  discussed  in  Chapter  Four. 

Other  limitations  of  the  program  are  the  inability 
of  the  user  to  define  a  varying  rail-slipper  gap,  a 
varying  velocity  profile,  and  the  requirement  for  model 
adaptation  for  use  with  SLEDYNE  versus  NASTRAN. 


1 .5  AREAS  OF  STUDY 


It  is  the  intent  of  this  study  to  define  new 
equations  of  motion  that  will  predict  dynamic  forces  in 
both  the  lateral  and  vertical  directions.  In  addition, 
the  equations  are  to  be  defined  such  that  dynamic 
transmitted  thrust  between  two  coupled  sleds  can  be 
predicted.  See  Figure  1.7  for  a  pictorial  representation 
of  this  concept.  The  derivations  of  the  sled  equations  of 
motion  are  to  be  carried  out  in  matrix  form  so  that  if 
only  one  sled  is  to  be  analyzed,  the  second  sled  and  its 
coupling  effect  on  the  first  can  be  ignored.  Finally,  it 
was  intended  to  reduce  the  equation  matrices  to  analyze 
monorail,  narrow  gage,  and  dual  rail  sleds,  either 
coupled  or  singularly. 

Another  improvement  intended  to  be  made  is  the 
ability  to  define  a  velocity  profile  for  input  into  the 
forcing  function  analysis.  Additionally,  variable 
slipper  clearance  definition  for  special  sleds  and  a  new 
slipper  stiffness  definition  to  eliminate  the  need  for 
model  adaptation  between  NASTRAN  and  the  existing  SLEDYNE 
are  planned. 

The  final  direction  of  study  will  be  to  correlate 
the  equation  developments  versus  real  time  force 
measurements.  This  correlation  will  be  accomplished  by 
comparing  predicted  forces  against  existing  monorail  data 
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FIGURE  1.7  PICTORIAL  REPRESENTATION  OF  COUPLED  SLEDS 


CHAPTER  2 


DERIVATION  OF  BASIC  EQUATIONS  OF  MOTION 

2.1  GENERAL 

The  basic  sled,  whether  dual  rail,  monorail,  'narrow 
gage,  or  outrigger,  has  five  (5)  modes  of  constrained 
rigid  body  motion  in  addition  to  the  downtrack  movement 
of  the  sled.  Based  on  Figure  2.1  and  the  given 
coordinate  system,  the  sled  has  the  following  rigid  body 
capabil ities  : 

1)  vertical  bounce— z 

2)  lateral  bounce - y 

3  )  pitch - 0 

y 

4)  yaw - 0 

z 

5)  roll-- - 0 

x 

Under  present  design  standards  it  has  been  assumed 
that  vertical  bounce  (z)  and  pitch  (e  )  are  the  only 
critical  design  factors  to  account  for,  and  the  other 
three  modes  of  excitation,  combined,  at  most  result  in 
loads  equal  to  this  form  of  motion.  As  discussed  in 
Chapter  One,  several  attempts  at  defining  rail  roughness 
forcing  functions  have  been  tried  using  only  this 
assumption.  However,  rat  is  factor y  verification  has  been 
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performing  a  modal  reduction  to  simulate  the  dynamic 
motion  of  the  sled. 

The  excitation  of  the  sled  is  due  to  the  five 
possible  rigid  body  modes  plus  the  elastic  modes  of  free 
vibration  participating  in  the  direction  of  these  rigid 
body  modes.  The  modal  reduction  is  based  on  using  only 
the  primary  modes  of  free  vibration  (i.e.,  the  first  two 
or  three  modes  in  each  rigid  body  direction)  for  each 
time  step  instead  of  solving  the  modal  equations  each 
time. 

2.2  DEVELOPMENT  OF  SIMULATION  EQUATIONS — ONE  SLED 

2.2.1  DEFINITION  OF  SLED  DISPLACEMENTS 

2. 2. 1.1  BASIC  VIBRATION  THEORY 

A  finite  element  model  of  a  sled  consists  of  a 
number  of  nodes  connected  by  elements.  Each  node  is 
assumed  to  have  six  (6)  degrees  of  freedom  (three 
translational  and  three  rotational).  The  dynamic 
equations  of  motion  for  free  vibrations  of  a  sled 
constrained  at  the  slippers  can  be  expressed  in  matrix 
form  as 

[m]{u}+[k]{u}=0  (1) 
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where  {u}  represents  the  vector  of  displacements  of  the 
degrees  of  freedom  (dof)  associated  with  the  sled.  The 
mass  is  represented  by  [m]  and  is  a  diagonal  matrix  of 
the  lumped  parameter  model  such  that  is  the  lumped 

mass  at  the  ith  node.  Finally,  [  k  ]  is  the  sled's 
stiffness  matrix. 

From  basic  vibration  theory^,  the  solution  of  the 
harmonic  equation  represented  by  equation  1  is  given  by 

{ u  }  =  {  $  }  cos  ( ojt  +  ip  )  (  2 ) 

where  {  $  }  is  the  matrix  of  eigenvalues,  co  is  the 
natural  frequency,  and  ip  represents  a  phase  angle.  If 
equation  2  is  differentiated  twice  and  substituted  into 
equation  1,  the  standard  eigenvalue  problem  results  and 
is  given  by 

uj2Cm]  {0  >  =  [k]  {4,  }  (3) 

Solution  of  this  equation  gives  the  appropriate 
eigenvalues  and  eigenvectors. 

2.2. 1.2  DERIVATION  OF  RIGID  BODY  DISPLACEMENTS 


The  mass  and  stiffness  matrices  shown  in  equation  3 


were  derived  under  the  assumption  that  the  sled  was 
supported  against  vertical  and  lateral  motion  by  the 
slippers.  For  a  true  r e p r e s e n t a t i o n  of  a  sled  tra¬ 
jectory,  this  motion  must  be  modified  to  include  the 
effects  of  rigid  body  movements.  This  enhancement  of  the 
sled  motion  can  be  expressed  as  the  sum  of  a  set  of 
displacement  functions.  These  functions  being  defined  are 
the  displacements  due  to  the  natural  sled  body 
frequencies  and  those  associated  with  the  rigid  body 
motions. 

In  order  to  define  rigid  body  movements,  consider 
Figure  2.2.  For  any  particle  P  that  is  offset  from  the 
sled’s  cen t er-o f-gr a v i ty  (denoted  as  CG  in  the  figure)  by 
direction  vectors  7,7,z,  there  will  be  displacements  due 
to  rigid  body  motions  of  the  CG.  This  can  be  expressed 
by  the  summation  of  the  translational  component  and  the 
translational  effect  of  the  rotary  term.  In  other  words, 
the  sled  displacement  in  the  X  direction  (designated  as 
u)  is  given  by 

u  =  x '  -  y  0  2  +  z  .  y  (  4  ) 

where  x'  denotes  the  translational  rigid  body  motion  in 
the  X  direction  and  1^  ^  represents  the  rotary  effect 
based  on  small  angle  theory.  it  must  be  noted  that  since 


the  sled  is  traveling  down  track  (X  direction),  the  rigid 


Consequently  the  rigid  body  displacements  can  be 
symbolized  in  matrix  form  as  shown  in  equation  7. 

2. 2. 1.3  DERIVATION  OF  DISPLACEMENT  EQUATIONS 

Superimposed  on  the  rigid  body  motion  displacements 
will  be  the  associated  deformations  caused  by  the  elastic 
motion  of  the  sled  under  natural  vibration.  Hence  the 
total  deflection  of  a  point  with  respect  to  time  can  be 
written  as  shown  in  equation  8,  where  q^(t),  r^(t),  and 
s^(t)  are  the  amplitude  components  of  £  ( x , y , z )  in  the 

X,Y,  or  Z  directions,  respectively. 

Since  the  sled  has  been  divided  into  a  finite  number 
of  elements,  the  final  displacement  matrix  can  be  written 
such  that  [D]  is  a  Nx3  matrix  where  N  represents  the 
number  of  nodes  in  the  original  finite  element  model  as 
shown  in  equation  9,  where  k  is  the  number  of  elastic 
modal  vibrations  being  used  and  must  be  less  than  the 
number  of  nodes,  N. 

2.2.2  EQUATIONS  OF  MOTION 

As  is  the  case  for  all  damped  linear  systems,  the 
equations  of  motion  can  be  written  as 


mx  +  cx  +  kx 


F 


(10) 


I 

where  m  is  the  system  mass,  c  is  the  damping  constant,  k 
is  the  stiffness  inherent  in  the  system,  and  F  is  any 
external  forcing  function  applied  to  the  system,  and  x 
represents  the  displacement  of  the  system. 

In  order  to  perform  a  dynamic  stress  analysis  on  a 
sled,  the  external  forcing  function,  F,  must  be  known. 
Part  of  the  forcing  function  is  known,  such  as  the 
aerodynamic  loads,  drag  forces,  thrust  vectors,  and 
accelerations  (the  q u a s  i -s t ead y-s t a t e  loads).  However, 
the  impact  forces  into  the  slippers  of  the  sled  are 
unknown  and  must  be  determined. 

To  solve  for  these  unknown  slipper  impact  forces, 
the  equations  of  motion  will  be  derived  using  Lagrange's 
equation.11  This  equation  may  be  written  as 


d _  3T 

dt  3X  . 


3W  ,  9R  ,  3  U 

-  +  - - -  +  - 

3X.  3X .  3X. 

111 


(11) 


where  T  is  the  kinetic  energy  of  the  system,  W  is  the 
work  done  by  the  system,  U  is  the  potential  energy,  R  is 
the  Rayleigh  dissipation  function,  and  Q  ^  are  the 
n on c o n s e r v a t i v e  forces  acting  on  the  system.  The  Xi 
represent  the  independent  variables  of  motion  y',  z',0  , 


0...  0,  .  and  q.  , 


r,,  s.  for  i  =  1  to  k. 


2.2.2. 1  KINETIC  ENERGY  DERIVATION 


The  kinetic  energy  of  a  system  under  displacement  i 
given  by 


1  /  2 

T  =  2  m v ^  (12) 

where  v  is  the  velocity  of  the  displacement.  In  terms  of 
matrix  notation  this  can  be  written  as 

•  • 

T  =  V?  t  D  ]  T  [  m  ]  [  D  ]  (13) 

where  [D]  is  the  velocity  of  the  system  and  is  the 
derivative  of  [D]  with  respect  to  time.  So  [ D ( t )  ]  can  be 
w-itten  as  shown  by  equation  14. 

Hence  the  kinetic  energy  term,  when  expanded  in 
matrix  form,  can  be  expressed  as 

T  :  7,  [  A  ]  T  [  m  ]  [  A  ]  (15) 

T 

w  h  e  •"  e  [A]  is  a  3xN  matrix  and  defined  in  equation  16. 

T 

The  differential  terms  of  [A]  can  be  factored  out 
T 

to  give  [A]  as  the  product  of  a  differential  matrix  and 

*  T 

a  constant  matrix  as  shown  in  equation  17,  where  [ n  ]  is 
a  1  x  (jN+6)  matrix  and  [  C  ]  ^  is  a  (3N+6)  x  3N  matrix. 

So  now  equation  15  can  be  rewritten  as 
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T  =  Vatn]  [C]l[m][C][r;] 


(  18) 


If  a  mass  property  matrix  [M]  is  defined  to  be 


[M]  =  [C]T[m][C]  (  19) 


then  the  kinetic  energy  of  the  system  can  be  expressed  as 


T  s  V2  [ h  ]  T  t  M  ]  [  *1  ]  (20) 


with  [M]  being  a  matrix  of  order  ( 3 N  +  6  )  x ( 3 N  +  6  )  .  For 


con f or m ab i 1 i t y  [m]  must  be  a  3Nx3N  matrix  and  of  the  form 


shown  in  equation  21. 


Now  the  elements  of  [M]  can  be  defined  as 
3N 

M11  =  ^  =  m  =  sled  mass  in  X  direction 


3N  2 

M  =  Z  m..(1)  =  m  =  sled  mass  in  Y  direction 

i  =  1  11 

3N  2 

M  =  z  m-.(1)  =  m  i  sled  mass  in  Z  direction 

35  i=1  11 

3  N  2  2 

M  =  z  m..(z  +y  )  =  mass  moment  of  inertia  about  X  axis 
i  =  1 

or  roll  inertia  =  I 

r 

3  N  2  2 

M  =  z  m..(x  +2  )  =  mass  moment  of  inertia  about  Y  axis 
55  i=1  ii 

or  pitch  inertia  =  I 

P 

3N  2  2 

^66  =  £  m  i  i  ( x  +y  ^  =  mass  rn°'nent  of  inertia  about  Z  axis 


or  yaw  inertia 


Or,  in  a  more  general  form 


Mpl  \!11CKp-6)]T[miiHCi(p-6)1  f<,r  PO>6  (22) 
Now  for  the  special  case  rfhere  p=l  this  reduces  to 


3N 

M  =  y  [C 
PP 


i  ( p-6 ) ]  [mii][Cl( p-6 ) ] 


(  23) 


which  is  the  generalized  mass  of  the  (p-6)th  mode  for 
p>6  . 

For  the  case  p  and  1 >6  and  p^l  equation  22  reduces 
to  zero  since  the  modes  are  orthogonal  to  each  other,  or 


Mpl  =  0  for  p/l,p>6,l>6  (24) 

For  the  case  where  p f 1  and  either  p  or  1<6  results 
in  the  equation  shown  by  equation  25. 

Finally  for  the  case  p^l,  and  p  and  1  are  less  than 
6,  results  in  Mp^  =  0  due  to  the  orthogonality  of  rigid 
body  motion. 

Thus  [M]  can  finally  be  written  as  indicated  by 
equation  26 . 

2. 2. 2. 2  POTENTIAL  ENERGY  DERIVATION 


The  potential,  or  strain,  energy  (U)  of  a  system  is 


given  by 


U  =  V2  kx2  (27) 

where  k  is  the  stiffness  and  x  is  the  generalized 
displacement  experienced  bj  the  system.  Using  previously 
defined  displacements  and  matrix  notation,  this  becomes 

U  =  V2  [D]T[k] [D]  (28) 

It  was  shown  earlier  that  the  velocity  matrix  was  a 
product  of  a  differentiated  general  displacement  times  a 
constant  matrix.  By  integrating  the  velocity  matrix  the 
displacement  matrix  is  obtained,  and  it  similarly  may  be 
expressed  as  the  product  of  a  general  displacement  matrix 
times  the  constant  matrix,  or 

[D]  =  [n  J[C]  (29) 

Therefore  the  potential  energy  may  be  expressed  as 
u  =  y,  [n  ]T[  K]  [n  ]  (30) 

where 

[K]  =  [ C]T[k]  [ C]  (31  ) 


(6 


and  [K]  is  a  ( 3N+6 ) x ( 3N+6 )  matrix 


This  matrix  has  elements  based  on  the  following  cases: 

1.  For  p  and  1<7»  rigid  body  motion  is  occurring  and  no 


strain  energy  is  associated  with  this  case,  K  =  0. 


2.  For  the  case  p  and  1>6  and  p^l,  orthogonality 


controls  and  again  K  ^  =  0 


3.  Finally,  for  p  and  1>6  and  p=l 


K 

PP 


tCi(p-6>)T[kii1CCl(p-6)1 


o  r 


K 

PP 


i=1.N  (32) 


2. 2. 2. 3  DAMPING  FORCE  DERIVATION 


The  Rayleigh  dissipation  function  governs  the  form 
of  the  damping  forces.  If  it  is  assumed  that  damping 
forces  are  proportional  to  the  velocity,  i.e.,  R=cV,  the 
Rayleigh  function  is  given  by 


R  =  V2  [r'  ]T[CD][n  ]  (33) 


where  [CD]  is  the  system's  damping  constants.  For 
simplification,  this  damping  matrix  is  assumed  to  be  a 
diagonal  matrix  of  order  ( 3 N  + 6 ) x ( 3 N  +  6  )  and  has  the  form 
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shown  in  equation  3*3. 


2. 2. 2. 4  DERIVATION  OF  WORK 


The  work  of  a  system  is  defined  as  the  product  of 
the  force  times  displacement.  For  this  problem,  work 


will  be  the  product  of  the  slipper  forces  times  the 
displacement  the  slipper  undergoes.  By  the  problem 
definition  the  structure  has  been  pinned  against  modal 
displacements  at  the  slippers.  Therefore,  only  the  rigid 
displacements  can  cause  the  slipper  "spring"  to  do  work. 

Earlier  the  rigid  body  displacements  for  any  point 
of  the  sled  were  derived,  equations  5  and  6.  Now  the 
points  on  the  sled  are  specified  to  be  the  slippers,  so 
the  values  of  y^,  and  zi  are  now  defined.  Let  1  and 
1  represent  the  forward  and  aft  distances  from  the  sled 
CG  to  the  slippers.  Similarly,  lr  and  1^  will  represent 
the  cross  track  distances  for  dual  rail  type  sleds. 
Finally  z  is  the  vertical  distance  that  the  CG  is  above 
the  slippers. 

From  these  definitions  the  vertical  slipper 
deflections  are  given  by 


=*  -  ]f  v  '  *1  x 


•.Vip,,, 


1 r  v  +  >i  x 


•V'Pn; 


z '  +  1  0  -  1  0 

ay  ax 


&VgaPz 


z '  +  1  ■ 1  +1  o 

ay  r  x 


v 

A  gaPj 


where  the  A  Vg  a  p  represents  the  change  in  slipper  gap 

due  to  varying  rail  properties.  The  subscripts  on  d  and 
v 

A  gap  represent  forward  (  F  )  ,  aft  (A),  left  (L),  or  right 
(R)  slippers  and  the  superscript  stands  for  vertical  (V) 
deflections.  Also  from  equation  6  the  lateral 
deflections  can  be  defined  as 


6fl 

-  y'  +  lf 

0*20  + 
2  X 

ALgapFI 

5fr 

=  y  ♦  if 

+  ’  i  + 

2  X 

aVpfr 

r  36) 

4 

=  V'  -  1 

a 

+  Z  o  + 

2  X 

■LsaPAi. 

4 

=  y  *  -  1 

7  a 

o  +  7  0  + 

2  X 

a4par 

From 

simple 

spring  theory,  the 

force  reacting 

spring  is 

given  by 

F  =  K 

A 

(  37  ) 

where  K  is  the  spring  stiffness  and  A  is  the  distance  the 
spring  displaces.  For  the  sled,  the  slipper  spring 
forces  occur  only  when  the  slipper  and  rail  have 
contacted.  This  happens  whenever  any  slipper 
displacement  exceeds  the  allowable  slipper  gap,'  .  Then 
the  slipper  forces  can  be  specified  as 
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F  (  i  )  =  K  (  i  )  iS  (  i  ) 


(  38  ) 


where  i  is  a  counter  of  one  to  eight  defined  in  Table  2-1 
and  7  (i)  is  the  difference  in  slipper  displacements  and 
slipper  gap,  or 


6  ( i  )  =  6  (  i  )  -  e  ( i  )  (39) 


TABLE  2-1  SLIPPER  COUNTER  DEFINITION 

i 

Definition 

1 

Forward  left  slipper,  vertical  force 

2 

Forward  right  slipper,  vertical  force 

3 

Aft  left  slipper,  vertical  force 

4 

Aft  right  slipper,  vertical  force 

5 

Forward  left  slipper,  lateral  force 

6 

Forward  right  slipper,  lateral  force 

7 

Aft  left  slipper,  lateral  force 

8 

Aft  right  slipper,  lateral  force 

When  rail  slipper  contact  does  occur,  damping  forces 
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are  associated  with  that  contact  so  the  slipper  force, 
equation  38,  must  be  modified  to 

F ( i )  =  K(i)6(i)  -  D(i)  (40) 

where  D(i)  represents  the  appropriate  damping  force. 

In  order  to  solve  for  the  slipper  reaction  forces 
two  things  must  be  known:  the  slipper  gap  and  its  change, 
and  the  damping  forces.  The  slipper  gap  is  simply  the 
nominal  clearance  built  into  the  slipper  assembly. 
Normally  this  clearance, e  ,  is  one-eighth  (0.125)  inch. 
To  perform  the  analysis  an  initialized  condition  must  be 
defined;  therefore  all  slipper  gaps  will  be  assumed  to  be 

e/2. 

The  change  in  slipper  gap  is  dependent  on  two  items. 
The  first  is  the  position  downtrack  and  the  rail 
alignment  at  that  position,  hereafter  referred  to  as  rail 
roughness.  The  second  item  is  the.  physical  change  in  the 
rail  cross-sectional  properties.  So  the  change  in 
3lipper  gap  is  written 

A  gap ( i )  =  Y  (  i  )  [  x  ]  +  y  (  i  )  C  x  ]  (41) 

So  the  slipper  force  dependent  on  the  slipper  gap 
will  be  zero  unles  the  following  relationship  holds 
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|  6(1)  -  e/2 1  >  0 


(  42) 


5  (  i ) 


By  maintaining  the  proper  sign  for  6  C  i )  ,  it  will 
establish  whether  rail  contact  was  made  on  the  upper, 
lower,  inner,  or  outer  slipper  surfaces. 

From  basic  structural  dynamics,  damping  forces  are 
generally  assumed  to  be  proportional  to  the  velocity,  or 
F  =  CV  For  the  travel  of  the  sled  there  is  the 
possibility  of  four  damping  forces  per  slipper,  which  can 
act  on  one  of  the  four  above-mentioned  surfaces. 

The  velocity  needed  for  the  damping  forces  will  be 
the  velocity  at  which  the  slipper  impacts  the  rail.  This 
velocity  is  dependent  on  the  direction  of  travel  of  the 
slipper,  and  has  been  previously  defined  in  terms  of 
displacements,  (  i  )  .  To  find  the  velocity  at  impact  it 
is  a  matter  of  differentiating  6(i)  with  respect  to  time, 
o  r 


V  (  i  )  =  ^(Mi)) 


(43) 


So  for  the  velocity  in  direction  one  (1)  this  is 
d 


V  ( 1 )  = 


dt  (Z'  -  Vy  -  Ll9x  +  Y(1)(X>  +  Y(l)(x) 


-  e/2) 


z'  1f°y  ~  1 1 0 x  +  Y(l)(x)  +  Y(l)(x) 
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(44) 


m 


o 


o 


L* 


» 


Since  1 f  and  1^  are  constants  no  contribution  to  the 
velocity  term  is  made. 

The  problem  with  equation  44  stems  from  defining  Y 
and  <$  .  This  is  due  to  both  Y  and  6  being  functions  of 
downtrack  position  and  not  time.  However,  these 
functions  can  be  made  t i m  e  -  d  e  p e n d e n t  since  X  =  Vt,  where 
V  represents  the  downtrack  instantaneous  velocity.  So 
Y ( x )  =  Y ( Vt  )  and 

V(x)  -  k-(Y(x>>  -  ^(Y(vt)>  -  iiHii  +  (45) 


Since  constant  velocity  is  assumed  for  a  given  time 
period,  the  differential  of  the  velocity  is  zero.  Now 
"  d  t  ^  ^  is  the  change  in  track  height  for  a  given  time 
increment,  or  it  can  be  called  the  local  slope.  Hence 

Y(x)  =  VY(x)  (46) 

where  x  represents  the  two  track  stations  encompassed  by 
the  time  interval,  t. 

Similarly  y(x)  can  be  shown  to  be 

Y(x)  =  Vy(iT)  (47) 


Finally,  the  slipper  dampl/ig  force  for  direction  one  (1) 


The  work  done  by  the  spring,  W,  may  be  given  as  the 
integral  over  the  displacement,  or 

T( i)8  _  7(1)8 

W  =  i  Z  [  F  ( i )  +  D(i)]d6(i)  =  f  I  K(i)  f(i)dc(i) 

0  i  =  1  0  i  =  1 

8  _ 

w  =  V2  K(i)id)2  (52) 

l  =  i 

2. 2. 2. 5  EVALUATION  OF  DERIVATIVES 

Now  that  all  of  the  terms  associated  with  Lagrange's 
equation  of  motion  (equation  11)  have  been  defined,  it  is 
possible  to  evaluate  the  partial  derivatives  and 
substitute  them  back  into  the  original  equation. 

2. 2. 2. 5.1  KINETIC  ENERGY 

From  equations  18  and  19  the  kinetic  energy  was 
shown  to  be 

T  =  V,  [  -*  ]  T  [  M  ]  [  n  ]  (53) 

or  in  non-matrix  form 

T  =  V;  M  r*  2  (  5  H) 
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2. 2. 2. 5. 2  POTENTIAL  ENERGY 

From  equation  30  in  a  non-matrix  form,  the  potential 
energy  was  shown  to  be 

U  =  I/,  K  n  2  (57) 

So 

—  =  Kn  =  | K 1 [ n ]  (58) 

2 . 2 . 2 .5 . 3  DAMPING  FORCE 


From  equation  31*  the  damping  force  was  derived  as 


but  K(i)  (i)  is  the  force  in  spring  i,  so 


;)w 


8 

Z  F 


i  =  1 


(i) 


3^(1) 

On 


(64) 


By  the  proper  substitution  of  <s  (i)  the  following  values 
for  c,  W  /  0  n  are  obtained  and  or  shown  in  equation  65. 


2. 2. 2. 5. 5  NON-CONSERVATIVE  FORCES 


The  non-conservative  forces  are  those  defined  by  the 
damping  forces.  These  forces  will  be  derived  by  finding 
the  virtual  work  done  by  the  damping.  In  a  mathematical 
expression  this  becomes 

N 

*wnc  =  .^1D(i)6ni(i)  =  Q(i)  (66) 

So  by  the  proper  substitution  of  6  (i)  the  matrix  shown 

in  equation  67  is  obtained  for  Q(i). 

2. 2. 2. 5. 6  FORCING  FUNCTION 

If  equation  11  is  rewritten,  a  forcing  function  F(t) 
can  be  defined  to  be 


{  F  { t )  } 


♦  Q  (  i  ) 

■1  x  . 

l 


(68) 


4 
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This  forcing  matrix  then  would  look  like  that  shown  in 


equation  69. 


2.2.3  SUMMARY 


If  equations  56,  58,  60,  and  68  are 

substituted  into  equation  11,  the  following 


combined 

results 


and 


[ M ]  {  fi  }  +  [  c ]  t  n  }  +  [  K ]  {  n  >  =  { F }  (70) 

where  matrices  [ M  ]  ,  [C],  [K],  and  { F }  are  defined  in 

equations  25,  3 4 ,  32,  and  69.  respectively. 

The  system  of  equations  has  been  developed  to 
include  rigid  body  motion.  It  is  now  necessary  to 
define  the  elements  of  the  damping  matrix,  [C],  of 
equation  70  corresponding  to  the  elastic  modes  of 
vibration.  Since  the  rigid  body  modes  are  uncoupled  from 
the  elastic  modes,  equation  70  can  be  written  as  the  sum 
of  the  rigid  motions  plus  the  elastic.  So  the  equation 
can  be  written  for  the  elastic  modes  as  shown  in  equation 
71,  where  equation  71  represents  3N  uncoupled 
differential  equations.  For  any  i  this  equation  can  be 
rewritten  as 


M  .  .  x  ,  . 
ii  (  i 


•ft) 


+  L i  i  X (  i  -  6  ) 


+  M 


(  i  -  6  )  '  (  i  -h) 


=  0  i  =  7  ,  3 


N  (72) 


and  x,.  ,  (  represents  either  q  ,  r  .  ,  nr  s  .  . 


From  damping  theory,  it  is  assumed 

[C]  =  2Y[M]  (73) 


so  equation  72  becomes 


M  .  .  x  ,  .  ,  .  + 

11  ( i-6 ) 


2 YM .  .  x  ,  .  , . + 

n  (i-6) 


M  .  .  to' 
1 1 


(i-6)X(i-6)“° 


i  =  7  ,8 


This  differential  equation  has  the  charact 
solution  of 

x, .  , .  =  AeSt  (75) 

(  l- 6  ) 

So  equation  74  can  be  written  as 

[s2  +  2  Y  s  +  w2( i.6) ]MiiAeSt  =  °  (76) 


..N  (74) 


r  i  s  t  i  c 


This  equation,  for  the  non-trivial  case,  has  roots 


For  the  critically  damped  case 


Y 


2 


2 

*  (  i-6  ) 


0 


(  78  ) 


where  Y  =  TQ  or  uj  .  =  Y  Q  and  Y  Q  is  the  value  of  damping 

for  the  critically  damped  system.  For  the  structural 
response  of  the  sled  the  underdamped  case  is  the  only  one 

of  interest,  or  Y / Y  <  1.  If  it  is  assumed  that  the 

o 

system  damping  ratio  is  the  same  as  the  viscous  damping 
ratio,  i.e.,  f>(i_6)  =  ^  *  then  the  damping 

coefficient  becomes 


i  i 


=  2  iM 


1 1 


=  2  v 


/  V  w 
o 


(i-6)  Mii 


(79) 


o  r 


Cii  =  2  (  i-6)  a,(  i-6)  Mii 


(80) 


Now  all  elements  of  the  generalized  coordinate 
differential  equation  have  been  defined.  The  only 
additional  modification  necessary  is  to  include  the 
effects  of  any  q u a s i - s t e a d y - s t a t e  forces  acting  on  the 
structure  at  the  time  of  interest.  These  forces  are 
defined  in  Appendix  B,  and  are  added  to  the  force 
definition  as  indicated  in  equation  81. 
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2.3  DERIVATION  OF  SIMULATION  EQUATIONS — COUPLED  SLEDS 

Due  to  the  nature  of  the  types  of  sled  couplings 
used  at  the  Test  Track,  the  equations  of  motion  for  the 
second  sled  are  identical  to  those  of  the  first  sled, 
provided  the  conditions  of  compatibility  are  maintained 
at  the  coupling  joint.  Appendix  B  contains  the 
derivation  for  the  q  u  a  s  i  -  s  t  e  a  d  y  -  s  t  a  t  e  forces  acting  on 
either  sled  that  are  necessary  for  defining  the 
differential  equations  of  motion  in  terms  of  the 
generalized  coordinates.  To  define  the  quasi-steady- 
state  force  terms,  the  coupling  compatibility  conditions 
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are  necessary.  These  have  also  been  derived  in  Appendix 

B. 

Using  the  results  of  coupling  compatibility,  the 
dynamic  simulation  of  sled  two  is  identical  to  that  just 
discussed  in  Chapter  2.2,  provided  that  the  proper  modal 
displacements  and  masses  are  used.  When  push  pad  or 
monoball  couplings  are  used  between  sleds,  the  type  of 
coupling  allows  for  decoupling  of  the  sleds  during  the 
modal  analysis.  This  means  that  individual  modal 
analyses  for  sleds  one  and  two  are  allowed. 

Preliminary  thought  indicated  that  using  a  Coleman 

coupling  between  sleds  would  not  allow  the  decoupling  of 

the  sleds  into  their  respective  modal  analyses.  This 

inability  to  decouple  the  sleds  led  to  the  investigation 

of  transmitted  thrust  transients  between  sleds  joined  in 

1  2 

this  manner.  Mixon  proceeded  to  investigate  this 
phenomenon.  He  concluded  that  the  were  due  to  motor 
burnout  spikes  exciting  the  downtrack  mode  of  the  pusher 
sled  and  this  being  transmitted  into  the  forebody  sled. 
In  other  words,  the  only  coupling  effects  between  sleds 
is  in  a  downtrack  mode.  Therefore,  individual  modal 
analyses  of  each  sled  may  also  be  accomplished  for  this 
t  ype  of  coupling. 

The  net  result  is  that  the  downtrack  equations  of 
motion  for  the  f  o  r  e  b  o  j  y  sled  a  o  e  dependent  on  the 
downtrack  motions  of  the  pushes  sled.  This  dependency 


between  sleds  means  the  equations  of  motion  of  the  front 
sled  will  have  to  be  modified.  Also,  the  equations  of 
motion  for  the  back  sled  will  have  to  be  solved  prior  to 
solving  them  for  the  front  sled. 


CHAPTER  3 


DASTARR 


3.1  GENERAL 


In  Chapter  Two,  the  equations  of  motion  for  a  system 
of  coupled  sleds  were  derived.  This  derivation  used  the 
normal  modes  of  free  vibration  for  the  sled,  the  natural 
frequencies,  the  generalized  masses,  and  the  diagonal 


matrix 

o  f 

the 

lumped 

masses  . 

These  quantities 

were 

assumed 

t  o 

b  e 

known  , 

and,  in 

fact,  they  are. 

These 

quantities  are  obtained  by  conducting  a  free  vibration 
analysis  of  a  finite  element  model  of  the  sled  that  has 
been  restrained  against  translation  at  the  slipper 
support  points.  At  present  the  eigenvalue  extraction  is 
accomplished  using  the  NASTRAN  computer  program.  These 
parameters  are  then  reformatted  into  data  acceptable  to  a 
program  called  DASTARR1^  (Dynamic  Analysis  of  a  Sled 
Traveling  Along  a  Rough  Rail). 

3.2  DASTARR  DESCRIPTION 

DASTARR  is  the  computer  program  written  as  part  of 
this  study  to  assemble  and  accurately  apply  the  equations 
of  motion.  DASTARR  computes  the  time  history  response  of 
a  sled  that  is  traversing  a  track  that  has  a  defined 


characteristic  rail  roughness  derived  from  sample  survey 
measurements  . 

DASTARR  uses  a  fourth-order  Ad ams-Mo ul t o n  predictor 
corrector  method  based  on  the  classical  fourth-order 
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Runga-Kutta  method  used  for  time  integration.  *  The 
rationale  for  using  this  integration  approach  is  to  use 
the  identical  method  as  in  the  existing  SLEDYNE.  The 
overall  intent  of  DASTARR  is  to  have  the  ability  to 
collapse  from  a  six  degree  of  freedom  system  into  the 
pitch  plane  system  of  SLEDYNE  and  yield  the  same  results. 

The  results  of  jASTARR  are  a  set  of  force  vectors 
and  accelerations  associated  with  each  slipper.  These 
forces  and  accelerations  are  then  converted  into 
equivalent  inertial  forces  acting  at  all  the  individual 
mass  points  of  the  sled.  These  inertial  forces  are  then 
recombined  with  the  original  finite  element  model  and  a 
static  analysis  can  be  performed. 

3.3  DASTARR  USER-SUPPLIED  INPUT 

DASTARR  was  written  to  be  a  general  purpose  program 
in  what  it  would  allow  for  simulation  variables.  For 
instance,  it  will  allow  for  either  a  variable  velocity  or 
a  constant  velocity  simulation.  It  will  examine  the 
"rigid  body"  system  response  for  gross  order  of  magnitude 
forces  and  then  accept  modal  participation  for  refinment 
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of  the  results.  It  allows  for  constant  quasi-steady- 
state  force  input  (simulating  constant  velocity);  or  time 
or  velocity  dependent  force  input  for  the  variable 
analysis.  The  variable  force  input  can  be  done  one  of 
two  ways.  It  can  be  input  as  a  tabular  function  or  an 
actual  trajectory  parameter  file  containing  the 
information  may  be  used.  Also,  the  slipper  force  output 
may  be  specified  as  either  q u a s i - s t e a d y - s t a t e  forces, 
dynamic  forces,  or  as  peak  slipper  forces. 

The  qua s i -s t ead y- s t a t e  force  inputs  are  actually  the 
nontransient  force  components  acting  on  the  sled.  These 
are  the  aerodynamic  drag  and  lift  forces  as  well  as  the 
thrust,  transmitted  thrust,  and  braking  forces  acting  on 
the  sled  at  that  given  velocity.  In  order  to  handle  the 
simulation,  the  program  must  know  the  magnitude  and 
location  of  each  of  these  forces.  Figure  shows  the 

graphical  representation  of  these  forces  and  associated 
locations  for  program  input. 

All  required  inputs  into  the  program  amount  to  an 

extensive  data  file.  An  interactive  prompt  program  was 

written  to  facilitate  the  data  input.  This  program 

prompts  by  asking  questions  about  necessary  data  and 

properly  formats  the  data  for  use  with  DASTARR.  This 

program  is  called  SLEDINP  and  is  documented  in  its  user's 
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manual  . 


3.4  DASTARR  OUTPUT 


DASTARR,  in  the  modal  participation  mode,  allows  the 
user  to  specify  up  to  five  displacement  degrees  of 
freedom  per  sled  to  be  monitored.  The  program  will 
provide  the  maximum  displacements  and  accelerations 
associated  with  these  degrees  of  freedom  that  were 
calculated  during  the  simulation.  This  output  is  very 
important  for  sled  design  as  some  sled  components,  such 
as  Velocity  Measuring  System  (VMS)  heads,  are  limited  in 
allowable  motion  and  this  feature  allows  for  monitoring 
the  effects  of  dynamics  on  such  an  item. 

As  previously  mentioned,  DASTARR  allows  the  user  to 
specify  slipper  force  ouput  type  versus  time.  Slipper 
forces  may  be  defined  as  dynamic  predicted  forces,  peak 
forces,  or  quasi-steady-state  forces.  Based  on  the 
specified  type  of  slipper  force,  DASTARR  will  create  a 
histogram  of  all  calculated  forces  during  the  simulation. 
The  histogram  contains  15  bins  for  positive  forces  and  15 
bins  for  negative  forces  p  e  r  slipper.  The  bin  size  is 
based  on  the  calculated  maximum  slipper  force  seen  in  the 
simulation.  The  histogram  feature  was  added  for  two 
reasons.  The  first  reason  was  to  identify  the  number  of 
maximum  peak  occurrences,  and  the  second  reason  to  1 1  low 
for  testing  of  data  distribution  in  the  bins.  Mixon 
proved^  that  dynamic  forces  w  e  r  e  normally  distributed  and 


hypothesized  that  peak  forces  were  Rayleigh  distributed. 
The  histogram  feature  will  allow  for  testing  of  these 
hypotheses.  This  is  discussed  further  in  Chapters  Five 
and  Six  for  dual  nail  and  monorail  sleds,  respectively. 

Finally,  DASTARR  creates  a  slipper  force  versus  time 
output  file  for  the  simulation  that  may  be  used  for 
external  data  man i pu 1  a t i on  .  Particular  examples  of  this 
usage  is  for  force  versus  time  plotting  and  for  creating 
Power  Spectral  Density  (PSD)  plots  of  the  analysis. 

3.5  PRELIMINARY  DASTARR  VERIFICATION 

In  order  to  verify  the  DASTARR  program  output,  a 
preliminary  study  was  done  to  compare  the  pitch  plane 
version  of  DASTARR  to  existing  SLEDYNE  results.  The  data 
base  for  this  comparison  was  the  Cruise  Missile  Sled 
model  used  on  the  ABIT  (Aircraft  Blast  Interaction  Test) 
series.  This  particular  sled  had  been  subjected  to 
extensive  simulation  on  SLEDYNE.  This  effort  was  to 
repeat  existing  test  parameters  for  comparison  purposes. 

Results  from  one  particular  situation  using  a  rough, 
periodic  rail  (Track  Station  3^000  top)  for  the  forcing 
function  input  are  summarized  in  Table  3-1.  Rail  forcing 
function  inputs  are  discussed  in  Chapter  Four.  In  the 
table,  there  are  four  cases  investigated  for  the  DASTARR 
These  may  be  described  as  follows.  Case 


hr> 


simulat ion . 


ounr 


number  one  is  the  derived  equations  of  motion  being 
exactly  as  in  the  existing  SLEDYNE;  however,  the  sign 
convention  for  pitch  plane  accelerations  are  in  the 
opposite  sense  for  DASTARR.  Case  two  is  the  DASTARR- 
derived  equations  of  motion  including  the  instantaneous 
downtrack  (x-axis)  acceleration  vector.  Case  three  is  the 
same  as  case  two,  only  no  modal  participation  is 
involved.  Finally,  case  four  is  the  DASTARR-derived 
equations  of  motion  with  matching  sign  convention  to 
SLEDYNE  and  no  downtrack  instantaneous  acceleration. 

Notice  from  the  table  that  slippers  one  and  three 
a^e  nearly  identical  between  SLEDYNE  and  the  DASTARR  case 
number  4.  The  slight  differences  are  attributable  to 
r-ound-off  error  j_n  the  integration  of  the  equations  of 
motion  between  computers  used.  The  SLEDYNE  version  is 
run  on  a  64-bit  machine,  whereas  the  DASTARR  version  is 
run  on  a  16-bit  machine.  However,  for  rails  two  and  four 
the  results  are  significantly  different.  Upon 
investigation  of  this  discrepancy,  it  was  found  that  the 
phasing  of  the  two  rail  forcing  function  start  points  is 
critical.  Under  the  existing  SLEDYNE  simulation,  this 
phasing  was  unknown  and  incapable  of  being  repeated  due 
to  the  random  start  point  generator  program  used  to  start 
the  functions.  At  this  point  a  mini-study  was  undertaken 
to  determine  exactly  what  the  effect  was  of  the  rail 
functions'  respective  start  points. 


It  was  decided  to  start  with  zero  phase  separation 
between  rails  and  proceed  by  phasing  the  second  rail  to 
start  in  20  feet  increments  behind  the  first  rail  until 
200  foot  of  separtion  had  been  achieved.  The  r  e  s  u  ►  s  of 
the  study  are  shown  in  Figure  3-2  for  the  sled's  front 
two  slippers.  Marked  on  the  figure  are  the  forces 
predicted  by  SLEDYNE.  Note  that  these  marks  are  placed 
in  both  the  positive  and  negative  directions  since 
SLEDYNE  only  looks  at  the  absolute  maximum  and  does  n  °  t 
maintain  signs.  As  can  be  seen  in  the  figure,  it  is 
possible  to  have  any  number  of  phasing  start  points  that 
will  meet  SLEDYNE  predictions  for  these  two  slippers. 
Upon  examining  all  four  slippers,  it  was  found  that 
slippers  one,  two,  and  three  could  be  adequately 
simulated  by  a  phasing  separation  of  approximately  150 
feet.  However,  slipper  four  could  not  be  simulated  in 
any  manner. 

Further  study  of  the  model's  modal  analysis 
indicated  that  slipper  four  was  very  sensitive  to  any 
kind  of  structural  input.  For  this  velocity,  the  rail  is 
imparting  a  fundamental  frequency  of  14.4  hertz  due  to 
the  long  period  undulation.  Actual  vibration  data  taken 
on  this  sled  showed  slipper  four  having  a  fundamental 
resonance  at  14.7  hertz.  From  this  it  was  concluded  that 
the  uncertainty  due  to  round-off  error  coupled  with 
slipper  four's  modal  sensitivity  at  this  f  r  e  q  u  e  :i  :  y 


DASTARR  Simulation  of  Phasing  Effects 


contributed  to  being  unable  to  completely  duplicate 
results  between  the  two  simulation  programs. 

3.6  CONCLUSIONS 

Based  on  the  apparent  adaquecy  of  simulation  between 

DASTARR  and  SLEDYNE  for  pitch  plane  analysis,  it  was 

concluded  that  DASTARR  provided  sufficiently  accurate 

results  to  continue  with  further  studies.  At  this  same 

time  period  another  dual  rail  sled's  finite  element  model 
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had  been  verified  through  modal  test  studies.  This 

sled  is  highly  instrumented  to  measure  slipper  forces  and 
a  sled  run  had  just  been  accomplished  with  verifiable 
slipper  data.  It  was  decided  to  replace  the  overly 
sensitive  Cruise  Missile  Sled  with  the  more  reliable 
Ramjet  sled  model  for  the  balance  of  the  simulation 
studies  for  dual  rail  sleds. 


CHAPTER  4 


RAIL  ROUGHNESS 

4.1  GENERAL 

Earlier  in  Chapter  One  the  limitations  of  the  rail 
forcing  function,  otherwise  known  as  rail  roughness,  were 
discussed.  Results  of  this  rail  profile  study  indicated 
that  not  enough  rail  roughness  measurements  had  been 
taken  (51  samples  spaced  one  foot  apart  every  2000  feet 
of  track).  The  study  indicated  that  if  a  low  frequency 
rail  undulation  (a  track  period  in  the  neighborhood  of 
200  feet  in  length)  existed  for  dual  rail  sleds  traveling 
at  3000  feet  per  second,  the  first  sled  structural 
frequency  (i.e.,  20  hz.)  might  be  excited.  Based  on  this 
Knowledge,  further  profile  measurements  were  directed. 
The  results  of  this  survey  indicated  there  was  one 
section  of  track  that  did  have  a  fundamental  period  of 
208  feet,  and  another  of  270  feet. 

From  this  rail  data  it  was  necessary  to  determine  if 
this  profile  could  excite  the  fundamental  resonances  of 
the  sled.  This  resulted  in  a  directive  to  perform 
additional  computer  simulations  with  the  new  rail 
profiles  to  define  the  potential  for  damaging  sleds  based 
on  this  profile.  The  objectives  of  this  special  study 


we»-e  four-fold,  being 


a.  To  establish  a  final  and  fixed  rail 
roughness  simulation,  if  possible,  for  SLEDYNE  or 
its  newer  version. 

b.  To  establish  a  fixed  rail  length  to 
optimize  the  existing  SLEDYNE  prediction  forces  for 
a  constant  velocity  run> 

c.  To  establish  if  long  period  undulations 
(track  periods  of  approximately  200  feet  in  length) 
introduce  significant  increases  in  impact  forces. 

d.  To  establish  if  rail  grinding  is  effective 
in  reducing  the  sled  impact  forces. 

4.2  RAIL  ROUGHNESS  STUDY 

These  requirements  resulted  in  developing  a  test 
matrix  for  performing  this  study  (Table  4-1).  This 
matrix  was  originally  set  up  to  use  the  existing  S L F. DYNE 
program.  The  tests  were  to  be  conducted  on  both  dual 
rail  and  monorail  sleds  traveling  at  a  constant  velocity 
for  that  particular  simulation,  but  varying  the  velocity 
between  simulations.  The  rail  profiles  to  be  used  are 
presented  in  Figures  4.1  through  4.8.  These  are  the 
results  after  removing  the  local  mean  from  the  data.  The 
data  passed  a  quality  control  process,  and  the  resultant 
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Figure  4.4  North  15000  Rail  Bottom 


Figure  4.8  TS  44000  Rail  Bottom 


s  t  a  ri  J  ti  r  of  the  estimates  (  5  )  for  each  sample 

y  ,  X 

set  f  fill  data  a  r  e  jh"wn  in  Table  4  -  1  . 

T  n  e  s  ••  fail  profiles  were  chosen  t  >•>  best  repr  e  s  <>  n  t 
the  actual  rail  conditions  and  still  answer  the  basic 
questions  pr*sed  by  the  study.  Track  Stations  (  T  S  )  340  j  0 

and  4  4  0 '!  0  ,  both  t«p  and  b^tt^m  of  the  rail,  contain  long 
period  undulations  on  the  o  r  J  ^  r  of  2  00  feet.  T5  4  4  0  00  is 
smooth  compared  to  the  rough  TO  j  4  0  0  0  r  a f  \  {  a  n  S  of 

y  > x 

0.012  inches  versus  0.043  inches).  In  addition,  the  south 
15  0  0  0  feet  of  r a i 1  top  had  the  same  standard  error  as  T3 
44000  rail  top,  which  has  a  periodic  content.  3y 
comparing  these  two  rail  functions,  the  effect  of 
periodicity  can  be  examined.  By  comparing  the  periodic 
rail  functions,  the  cumulative  effect  of  roughness 
combined  with  periodicity  can  be  investigated. 

By  comparing  the  north  and  south  15000  feet  of 
track,  rail  grinding  effectiveness  can  be  determined. 
The  north  15000  feet  has  been  extensively  ground  over  the 
last  decade,  whereas  the  south  15000  feet  has  been 
continually  aligned  and  straightened  over  the  past  <■' 
years  . 

4.3  SLEDYNE/CRUIoE  Mi  TO  ILK  BLED  SIMULATIONS 

The  previously  mentioned  0  r  u  i  s  e  Missile  Sled'  3  L  E  l1  V  N  hi 
s  t  u  d  y  o  f  C  h  a  p  t  e  r  Three'  w  a  s  i  dir  0  *  ,• ,,  ...  u  1  4  of  f  h  i  s  *-  a  i  1 


roughness  study  directive.  Shortly  into  the  simulation 
sequence  it  was  found  tnat  the  existing  SLEDYNE  rail 
forcing  function  (referred  to  hereafter  as  the  baseline 
rail  function)  was  the  controlling  factor  for  rail 
length.  Slipper  dynamic  forces  peaked  at  a  D  o  u  t  4000  feet 
of  track  simulation  for  the  baseline  function,  and  at 
lengths  slightly  less  for  the  other  rail  functions. 

At  this  point  in  the  study  it  became  apparent  that 
the  user-supplied  quasi-steady  lift  forces  would 
drastically  affect  the  slipper  results.  The  results  of  a 
mini-lift  study  showed  that  using  the  rough  periodic 
rail,  TS  3  4  0  0  0  top  of  rail,  slipper  forces  peaked  in 
p  n  a  s  e  at  a  lift  of  3  0  0  0  pounds  for  the  3  000  feet  per 
second  velocity.  F o  r  a  velocity  of  1  p  0  0  feet  per  second, 
lift  forces  of  3000  pounds  minimized  slipper  forces  and 
again  all  were  in  phase  with  each  other.  For  a 
n  ■  •  n  periodic  r a i  1  function  (the  baseline  function),  the 
lift  did  not  appear  to  affect  slipper  forces  as 
!  r  , a  tic  ally.  The  slipper  forces  stayed  in  an  operating 
e n  v  e ’  o  p  e  of  about  +1000  pounds  regardless  of  lift  or 
v  •  •  !  o  o  i  t.  y  .  No  discernible  phasing  of  slipper  forces  were 
1  lent.  Figures  4.9  and  4.10  show  the  graphical  results 
-  1  !i  "  Lift  study  on  the  Cruise  Mi  $s  i  l  e  C>  1  e  1  . 

T  n  e  b  1  lance  '  f  the  study  :''ncentratj  d  on  the  effect 
of  '  n  g  per  i<>  l  r ,,  undulations  and  used  T  S  3  4  0  0  0  and  T  S 
4  a  i  1  0  rail  f  u  n  i' t  1  o  n  3  .  >  juluc-il  re:.d  1  is  of  ►  a  e  3  (>  studies 
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FIGURE  4.10  LIFT  STUDY  VELOCITY  AT  1500  FPS 


are  shown  in  Figures  4.11  and  4.12.  Upon  examination  of 
the  figures,  several  inferences  may  be  made.  These  are 
summarized  below. 

a.  The  rougher  the  rail,  the  larger  the  slipper 
forces  . 

b.  The  higher  the  velocity,  the  larger  the  slipper 
forces . 

c.  Low  amplitude,  or  "smooth",  long  period 
undulation  (i.e.,  200  feet)  does  not  appear  to 
excite  the  sled . 

d.  Large  amplitude,  or  "rough",  long  period 
undulation  does  appear  to  excite  the  sled 
considerably  . 

4.4  CONCLUSIONS  AND  RECOMMENDATIONS 

No  conclusions  were  possible  at  this  time  concerning 
rail  grinding  effectiveness  nor  the  extent  of  periodic 
content  on  slipper  forces.  Computer  simulation  costs 
curtailed  the  rail  roughness  study  at  this  point  without 
answering  the  posed  questions  outside  of  optimal  rail 
length.  Here  it  was  decided  to  suspend  the  study  until 
DASTARR  was  operational  and  then  repeat  the  study  to 
answer  questions  left  unanswered. 

The  approach  to  be  used  for  further  study  will  be  to 
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use  the  same  rail  function  for  all  rail  surfaces  during 
the  simulation.  To  insure  repeatability  the  start  point 
will  be  controlled  and  not  randomly  generated.  The  start 
point  will  be  controlled  for  the  lateral  surface  of  rail 
number  one.  Each  successive  rail  surface  (number  one 
vertical,  number  two  lateral,  and  number  two  vertical) 
will  start  in  increments  of  20  feet  beyond  the  start 
point  of  its  predecessor.  To  generate  a  rail  of  4000 
feet  in  length,  the  rail  functions  will  be  looped  end  to 
end  until  this  distance  of  simulated  track  has  occurred. 

The  results  of  the  study  as  well  as  program 
verification  will  be  presented  in  Chapter  Five  for  dual 
rail  sleds  (the  Ramjet  model)  and  in  Chapter  Six  for 


monorail  sleds 


CHAPTER  5 


DUAL  RAIL  SLED  SIMULATION  — RAM  JET 


5.1  GENERAL 

5.1.1  RAMJET  SLED  DESCRIPTION 

The  Ramjet  Mission  Sled  is  a  special  purpose  sled 
designed  to  carry  either  a  ramjet  or  scramjet  engine  in  a 
simulated  environment  at  speeds  up  to  Mach  3.5.  The 
engine  is  isolated  from  the  sled  structure  through  a 
series  of  coiled  wire  isolators  that  are  designed  to 
operate  at  a  10  hertz  frequency.  This  frequency  was 
chosen  to  simulate  the  aircraft  environment  and  also  to 
safely  fall  into  an  open  envelope  of  sled  operating 
frequencies. 

Due  to  the  unique  geometry  of  the  sled  (see  Figure 
5.1),  the  use  of  conventional  slipper  beams  was  not 
possible.  The  sled  was  mounted  on  belleville  springs  to 
keep  the  structural  response  of  the  sled  within 
historical  dual  rail  sled  response  (i.e.,  20  hz).  These 

springs  are  designed  such  that  the  slipper  hangars  are 
connected  to  the  slipper  through  the  springs  via  a 
special  instrumented  slipper  pin.  The  springs  were  sized 
to  impart  a  20  hertz  heave  frequency  and  be  able  to  carry 
loads  up  to  65  kips  each  on  the  front  slippers  and  up  to 

88 


35  kips  each  on  the  rear. 

Using  the  known  sled  geometry  a  math  model  was 

created  to  define  eigenvectors,  eigenvalues,  and  mass 

properties  for  use  in  the  simulation  program  DASTARR. 

The  math  model  was  defined  to  run  on  MSC  NASTRAN  version 

63-  The  resultant  model  is  depicted  in  Figure  5.2.  The 

validity  of  the  math  model  has  been  tested  through 

vibration  testing  to  verify  the  predicted  eigenvectors 

and  eigenvalues.  The  results  are  summarized  in  an 

1  7 

in-house  test  report. 

Special  instrumented  slipper  pins  were  used  to  mount 
the  sled  to  the  sprung  slippers  to  measure  actual  slipper 
forces.  Data  from  these  pins  were  verified  by  first  cal¬ 
ibrating  the  pins  in  a  tensile  test  machine  using  the 
actual  sled-borne  instrumentation.  The  pins  were 
calibrated  to  within  one  (1)  percent  of  known  load 
applications  in  a  lab  setting.  The  pins  were  recalibrated 
on  the  track  prior  to  running  the  sled  to  further  insure 
the  slipper  pin  response.  The  calibration  was 
accomplished  at  the  receiving  end  of  the  telemetry 
signal.  A  post-run  calibration  was  done  on  the  system  to 
insure  no  signal  degradation  occurred  during  the  run. 
The  system  passed  all  calibration  tests  and  documentation 

results  of  the  Ramjet  Sled  qualification  are  contained  in 

1  8 

the  in-house  report  presently  being  written. 
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5.1.2  SLED  FORCES 


Using  the  results  of  the  mini-lift  study  mentioned 
earlier,  it  was  known  that  the  lift  and  center-of-lift 
vectors  would  have  to  be  well  defined  as  a  quasi-steady- 
state  force  in  order  to  accurately  simulate  the  Ramjet 
test  results.  The  quasi-steady  loads  are  those  that  are 
slowly  time  varying  and  are  composed  of  thrust,  drag, 
braking,  transmitted  thrust,  lift,  and  the  downtrack 
inertia  required  for  equilibrium.  Dynamic  forces  are 
characterized  by  oscillatory  forces  with  a  frequency 
content  in  excess  of  six  to  ten  hertz.  The  dynamic 
forces  are  generally  composed  of  the  rail  roughness 
forcing  function,  unsteady  aerodynamics,  and  thrust 
transients  like  those  seen  at  motor  burnout. 

5. 1.2.1  LIFT  FORCE  CALCULATION 

In  February  of  1984,  the  Ramjet  sled  was  run  on  the 
track  (the  1 9 Y  —  C 1  Mission)  to  a  peak  velocity  of  about 
2000  feet  per  second  (M=1.81).  An  extensive  amount  of 
data  was  taken  and  analyzed,  among  which  were  the  four 
slipper  forces.  The  data  from  the  slippers  were  passed 
through  a  two-hertz  filter  to  remove  the  "dynamic" 


portion  of  the  forces 


The  remaining  data  were  digitized 


and  plotted  and  consisted  of  the  cumulative  slipper 
forces  caused  by  the  "quasi-steady-state"  loads.  The 
results  of  this  operation  for  this  run  are  shown  in 
Figure  5.3* 

Based  on  a  given  velocity,  the  instantaneous  slipper 
quasi-steady  force  is  known  from  the  run  data.  Using  the 
trajectory  analysis  for  the  same  velocity,  the 
instantaneous  drag,  thrust,  and  inertias  are  known.  The 
lift  and  its  center  can  be  calculated  by  solving  the 
equilibrium  equations  around  the  sled  freebody.  This 
procedure  was  followed  for  nine  (9)  thrusting  conditions 
and  for  eight  (8)  coasting  conditions  observed  during  the 
1QY-C1  run.  The  results  of  the  quasi-steady 
instantaneous  loads  are  shown  in  Table  5-1. 

Notice  from  the  table  there  are  six  center s-of-1 i ft 
(L)  that  are  flagged.  These  centers  are  physically 
unlikely,  and  they  happen  to  occur  at  an  instant  where 
extreme  transients  are  present.  These  transient 
conditions  are  explained  in  the  table.  Therefore,  these 
points  were  not  useH  in  any  further  analysis.  As  a  point 
of  interest,  notice  the  difference  in  the  c e n t e r -of -1 i f t 
locations  between  thrust  and  coast  conditions  for  the 
same  Mach  number.  The  center  of  lift  under  thrust  is 
forward  on  the  sled  from  that  under  coast.  This  is 
explained  by  the  bogie  pusher  beam  used  for  this  sled. 
See  Figure  5.^  for  a  schematic  of  the  sled  train 
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Figure  5.3  19Y-C1  Mission,  QSS  Loads  Vs.  Time 
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configuration.  Films  of  the  run  showed  the  bogie  push 
beam  creating  a  standing  shock  wave  at  the  rear  slippers 
of  the  Ramjet  sled  which  forced  the  lift  center  slightly 
forward . 

Now  that  all  valid  lift  conditions  were  known,  the 
simulation  for  the  Ramjet  could  continue.  This  analysis 
was  done  by  looking  at  the  rigid  body  response  at 
constant  velocity,  and  then  by  examining  the  modal 
participation  again  at  constant  velocity. 

5.2  RIGID  BODY  SIMULATION 

The  concept  of  rigid  body  simulation  is  actually  a 
misnomer:  it  is  in  reality  confining  the  modal 
participation  to  the  first  vertical  translational  sled 
modes.  Typically  for  dual  rail  sleds,  the  rigid  body 
modes  are  defined  by  the  slipper  beams'  first  vertical 
bending  modes.  However,  with  the  Ramjet,  the  rigid  body 
modes  are  defined  by  the  belleville  spring  system  used 
for  i3ol ation  . 

The  simulation  consisted  of  running  all  eleven  valid 
lift  conditions  for  a  track  distance  of  4000  feet  for 
each  defined  rail  function,  resulting  in  99  simulation 
runs  for  the  rigid  body  study.  The  subsequent  analysis 
compared  the  predicted  peak  slipper  forces  to  the  actual 
peak  slipper  forces  scaled  from  oscillograph  records. 
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Table  5-2  contains  the  scaled  slipper  peaks  from  the 
oscillograph  records. 

Notice  that  for  slippers  three  and  four,  the  aft 
beam,  the  values  peak  at  about  half  of  what  they  are  for 
the  front  beam.  As  a  result,  the  aft  beam  forces  are 
subject  to  a  larger  error  in  comparison  results  due  to 
seal ing  factors  . 

Each  rail  roughness  simulation,  consisting  of  the 
eleven  conditions,  was  statistically  averaged  for  each 
slipper  and  beam  combination.  The  results  of  this 
sampling  are  presented  as  a  percentage  deviation  between 
actual  and  predicted  forces,  and  also  the  standard 
deviation  (1-sigma  value)  about  the  average  deviation. 

5.2.1  DATA  ANALYSIS  — AMPLITUDE  DOMAIN 

5 . 2 . 1 . 1  PEAK  FORCES 

The  simulation  consisted  of  99  computer  runs  of  the 
eleven  conditions  for  each  of  the  nine  rail  profiles. 
Results  of  the  peak  force  simulations  are  shown  plotted 
in  Figures  5.5  through  5.8  for  slippers  one  through  fou'-  , 
respectively.  The  solid  line  is  the  actual  measured  peak 
force  for  that  condition.  The  interesting  thing  to  note 
is  that  slipper  forces  tend  to  exceed  the  measured  value 
for  rank  roughness  in  excess  of  0.020  inches. 
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Figure  5.5  Slipper  1  Peak  Force  Vs.  Predicted  Force 
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Figure  5.6  Slipper  2  Peak  Force  Vs.  Predicted  Force 
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Figure  5.8  Slipper  4  Peak.  Force  Vs.  Predicted  Force 


The  statistical  deviation  analysis  is  shown  plotted 
in  Figure  5.9  and  summarized  in  Table  5-3  for  the  front 
slipper  beam.  The  "  x  "  denotes  where  the  average 
deviation  for  the  eleven  simulations  occurs.  The  heavy 
arrowed  vertical  line  indicates  the  +2-sigma  value  about 
the  average  . 

Again  the  break  point  on  rank  roughness  occurs  at 
about  0.020  inches  of  roughness.  This  break  point  is 
where  the  average  starts  o v e r p r e d  i  c t i n g  the  forces 
instead  of  underpredicting. 

5. 2. 1.2  PEAK  FORCES  VERSUS  TIME 

One  output  feature  of  DASTARR  is  slipper  force-time 
history.  This  output  file  can  be  directly  plotted  to 
examine  what  is  actually  occurring  during  the  computer 
simulation.  Two  of  the  data  plots  generated  during  the 
1 9 Y  —  C 1  mission  were  slipper  force  versus  time  plots  for 
slippers  two  and  four  during  a  thrusting  condition  at  an 
average  velocity  of  1650  feet  per  second.  This  condition 
happened  to  coincide  with  an  actual  computer  simulation 
condition,  so  a  direct  comparison  was  possible.  Figures 
5.10a  and  5.11a  show  the  actual  force-time  history  for 
slippers  two  and  four,  respectively.  Shown  in  Figures 
5.10b  and  5.11b  are  the  simulated  force-time  plots  for 
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the  baseline  rail  function.  The  other  eight  rail 
function  histories  are  contained  in  Appendix  C,  Figures 
C.1  through  C.16.  It  should  be  noted  here  that  the  sign 
convention  for  the  actual  run  data  is  in  an  opposite 
sense  to  that  of  DASTARR. 

Notice  for  slipper  two  the  peak  slipper  forces 
stayed  in  an  envelope  between  6  and  41  kips  and  were 
centered  about  20  kips.  This  indicated  that  during  this 
time  slice  there  was  no  flight  time  associated  with  this 
slipper,  or  the  slipper  was  in  contact  with  the  rail  at 
all  times.  Looking  at  the  plots  for  the  various  rail 
surfaces  for  this  slipper  indicated  a  similar  centering 
about  20  kips  with  various  envelope  sizes.  The  envelope 
size  appeared  to  be  directly  related  to  rail  roughness. 
The  rougher  the  rail,  the  larger  the  envelope.  The 
baseline  data  did  give  a  few  occurrences  of  flight  time; 
however,  it  was  minimal  compared  to  total  rail  contact 
time  . 

Examining  actual  force  data  for  slipper  four 
indicates  the  forces  fluctuating  by  about  _+ 6  kips.  The 
baseline  data  has  the  same  fluctuating  characteristics; 
however,  the  envelope  is  about  +_7  kips.  All  other 
profiles  are  one-sided  in  the  operating  envelope  (i.e., 
the  centering  about  0  kips  does  not  occur)  with  a  large 
portion  of  flight  time.  Forces  approached  or  exceeded 
the  measured  peak  for  rnugnesses  in  excess  of  0.020 
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inches. 


The  rough  rails  (those  in  excess  of  0.020  inches) 
had  significant  numbers  of  occurrences  in  the  vicinity  of 
the  peak  spike  for  both  slippers.  However,  the  number  of 
spikes  in  the  vicinity  of  the  peak  from  the  actual  run 
are  few  compared  to  the  rough  rails. 

An  interesting  siderote  came  from  inspecting  the 
plots  for  the  south  15000  rail  top  data  for  slippers  one 
and  two.  There  is  a  distinct  cycle,  or  spike,  occurring 
approximately  every  quarter  second,  with  the  first  spike 
occurring  about  50  milliseconds  into  the  simulation. 
This  phenomenon  repeated  itself  for  all  conditions  using 
this  rail  function.  The  quarter  second  interval  at  this 
velocity  corresponds  to  a  track  length  of  about  400  feet. 
This  particular  rail  profile  is  402  feet  in  length.  At 
first  it  was  thought  that  this  was  a  facsimile  created  by 
looping  the  data  end  to  end  with  a  built-in 
discontinuity.  Reviewing  the  rail  roughness  data  showed 
there  was  no  discontinuity  created  by  looping  of  the 
data. 

The  next  plausible  explanation  was  a  localized  rough 
spot  in  the  rail  that  was  looping  on  itself  at  the  400 


foot  cycle.  Taking  into  account  the  forced  starting 
point  for  the  rail  function  and  the  distance  to  the  first 
spike,  it  was  calculated  the  rough  spot  should  be  about 
275  feet  into  the  data  file.  Examining  the  data  file 


[»1 


showed  three  points  (  277  -27  9  feet  into  file)  that  had  a 


roughness  twice  that  of  the  surrounding  points  on  either 
side  of  it  and  a  roughness  three  to  four  times  that  of 
the  overall  sample.  This  theory  was  further  validated  by 
the  presence  of  a  secondary  spike  occurring  about  50  feet 
from  the  first.  Not  knowing  direction  of  rail 
generation,  it  was  determined  there  should  be  another 
rough  spot  at  either  data  point  225  or  325.  Further 
review  of  the  data  indicated  another  rough  spot  at  data 
points  224  and  225  that  were  about  250  percent  larger 
than  the  total  rail  roughness. 

Returning  to  original  survey  data,  it  was  found  that 
the  roughness  generated  at  the  277-279  data  points  were 
coincidental  with  track  stations  12027-12029.  The  time 
that  the  exceptionally  large  force  spike  occurs  in  Figure 
5.10a  for  the  front  slipper  corresponds  to  an  area  of 
track  in  the  region  of  TS  12020-12035,  which  contains  the 
measured  rough  spot. 

The  fascinating  aspect  of  this  phenomenon  is  the 
capability  of  identifying  a  definite  track  location  for 
being  ground.  Prior  to  grinding,  the  overall  effect  on 
slipper  forces  can  be  examined  based  on  a  predetermined 
rail  smoothness.  By  looking  at  Figure  C.1  and  removing 
the  peak  force  spikes,  it  appears  possible  to  lower  the 
peak  forces  by  approximately  10  percent  by  smoothing  the 


two  discussed  localized  rough  spots 


A  follow-on 


simulation  was  made  where  the  identified  rough  spots  were 
aribtrarily  smoothed  to  0.020  inches  to  verify  this 
hypothesized  peak  force  reduction.  This  plot  is  shown  in 
Figure  5.12  for  this  modified  rail  function. 

This  modified,  or  smoothed,  rail  only  reduced  peak 
forces  by  4  percent  instead  of  the  estimated  10  percent. 
The  smoothed  rail  did  reduce  the  dynamic  range  of  the 
forces  by  17  percent.  This  in  itself  is  significant  in 
reducing  the  fatigue  amplitude  effects  on  the  sled. 
Looking  at  the  forcing  function  data  again,  several 
additional  rough  spots  on  the  order  of  0.030  inches  were 
found  that  corresponded  to  the  predicted  peaks  shown  in 
the  figure. 

5.2. 1 .3  FORCE  HISTOGRAM 

Discussed  in  Chapter  Three  was  the  histogram  feature 
of  DASTARR.  The  histogram  will  take  the  force- time 
history  for  a  slipper  and  put  each  data  point  into  one  of 
the  30  distinct  bins.  The  bins  were  set  up  to  allow  for 
15  bins  to  hold  the  positive  force  occurrences  and  15  for 
negative  forces.  The  bin  size  was  set  to  be  an  integer 
multiple  of  the  predicted  maximum  force  from  all  eight 


0  . 


Mixon  proved  that  dynamic  forces  were  normally 
distributed  and  hypothesized  that  peak  forces  were 
Rayleigh  distributed.  Using  the  histogram  feature,  tests 

were  set  up  to  prove  or  disprove  this  hypothesis. 

1  Q 

Examination  of  literature  7  indicated  that  a 
Rayleigh  distribution  is  a  special  subset  of  a  normal 
distribution,  hence  could  be  tested  as  such.  The 
standard  test  for  a  normal  distribution  is  with  the 
chi-squared  (  “ )  goodness  of  fit  test.  The  chi-squared 

test  will  test  the  actual  discrete  distribution  of  the 
data  for  a  small  range,  such  as  what  is  contained  in  an 
individual  bin. 

Further  literature  search  revealed  a  second,  less 

widely  usqd  but  more  accurate  test.  This  is  the 

?  0 

Kolmogorov-Smirnov  test.  This  test  concentrates  on  the 
deviation  between  a  predicted  "cumulative"  distribution 
function  and  the  observed  one. 

Both  tests  were  programmed  to  give  results  in  the 
histogram  section  of  DASTARR.  Tables  5  -  4  a  and  5  -  4  b  show 
typical  DASTARR  output.  This  example  was  generated  using 
a  constant  velocity  of  1138  feet  per  second  under 
coasting  conditions  for  the  south  15000  feet  of  track  on 
the  top  rail.  Figure  5.13  shows  the  plotted  values  for 
slipper  one  for  the  histogram  (the  step  function)  and  the 
predicted  true  Rayleigh  distribution. 


I 

I 


%y. 


It  is  inherent  in  the  graph  that  the  trend  is  toward 
a  Rayleigh  distribution  except  for  the  first  two  bins. 
This  discrepancy  in  the  first  two  bins  made  the  analysis 
fail  both  tests  for  goodness  of  fit  of  the  data. 
Reference  to  statistical  acceptance  criteria  for  a  data 
sample  this  large  (approximately  1500  points)  indicated 
the  histogram  should  have  35  bins  to  handle  the 
distribution  instead  of  15  as  is  done  due  to  computer 
memory  limitations. 

Noting  the  inherent  Ra y le i gh- 1 i ke  trend  of  the  data, 
one  last  test  was  performed.  A  data  sample  is  normally 
distributed  when  9  9  percent  of  the  data  is  contained 
within  -►  3-sigma.  A  routine  was  added  to  the  program  to 
determine  how  many  standard  deviations  were  required  to 
encompass  99  percent  of  the  data  at  the  upper  bin  limit. 
As  can  be  seen  from  Table  5  -  9  b  ,  the  order  of  three  and 
one -half  (3.5)  sigma  encompassed  the  necessary  data. 

Making  note  of  the  Ra y 1 e i g h - 1  i  ke  trend,  the  bin  size 
limitations  that  caused  test  failure,  and  99  percent  of 
the  data  being  encompassed  by  3.5  sigma,  it  was  concluded 
that  the  distribution  for  peak  forces  was  truly  a 
Rayleigh  distribution.  Examination  of  the  other  loading 
conditions  indicated  the  same  results.  However,  the 
extent  of  failing  the  acceptance  testing  was  strongly 
dependent  on  the  magnitude  of  the  quasi-steady  loads.  The 
bin  sizes  became  very  large  and  all  the  data  were 
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have  to  be  able  to  encompass  the  actual  force  data  range 
only.  Each  slipper  would  require  a  separate  bin 
distribution,  and  due  to  computer  memory  limitations  this 
distribution  is  impossible  to  do. 

5.2.2  DATA  ANALYSIS — FREQUENCY  DOMAIN 

In  the  previous  section,  the  data  output  had  been 
compared  to  actual  run  data  in  the  amplitude  domain.  It 
remained  to  be  analyzed  in  the  frequency  domain  to 
determine  whether  the  frequency  and  energy  content  of  the 
simulation  was  accurate. 

The  ff'rce-time  history  output  from  DASTARR  was 

passed  through  a  power  spectral  density  (PSD)  program  to 

examine  the  sled  response  in  the  frequency  domain.  The 

2 

y,  or  powc",  axis  for  standard  PSD’s  have  units  of  G  /Hz. 

p 

However,  for  this  analysis  the  units  are  lb^./Hz.  Over  a 
given  frequency  band,  the  square  root  of  the  area  of  that 
band  is  a  measure  of  the  root  mean  square  of  the  slipper 
force.  The  cumulative  area  of  the  PSD  is  a  measure  of 
the  total  power,  or  energy,  of  the  system. 
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Due  t  o  the  number  of  data  points  being  sampled  and 


n 

tue  time  differential  between  samples,  approximately  27  5 
hertz  is  the  maximum  frequency  that  can  be  generated. 
F-r  most  dual  rail  sled  systems,  frequencies  above  200 
hertz  generally  have  such  small  deflections  associated 
with  those  frequencies  that  little  or  no  mechanical 
damage  results. 

Examination  of  PSD  plots  for  the  given  run  condition 
under  thrust,  shown  in  Figures  5.14a  and  5.15a  for 
slippers  two  and  four,  shows  a  sharp  rolloff  after  30 
hertz  and  mostly  noise  above  70  hertz.  This  phenomenon 
is  typical  of  most  sled  structural  responses.  Total  power 
in  the  front  slipper  is  20374  pounds  force  rms  for  the 
frequency  range  0  through  300  hertz.  The  aft  slipper 
total  power  is  1745  pounds  force  rms  for  the  same 
frequency  spectrum. 


Figures  5.14b  and  5.15b  show  typical  wide  spectrum 
response  for  the  rigid  body  sled  system  for  the  same 
thrusting  condition  from  the  simulation.  Again  sharp 
rolloff  is  shown  above  30  hertz  with  noise  starting  about 
70  hertz.  The  local  PSD  routine  used  did  not  contain  the 
feature  to  accumulate  power,  so  a  total  system  rms  power 
comparison  is  not  possible. 

Figures  5.16a  and  5.17a  are  blowups  of  5.14a  and 
5.15a,  respectively,  in  the  0  to  80  hertz  range. 


Similarly  Figures  5.16b  and  5.17b  are  blowups  of  5.14b 


and  5.15b  for  the  same  frequency  spectrum,  only  clipping 

4  6  5 

the  power  axis  at  10  and  10  .  With  these  figures  it  is 
possible  to  isolate  individual  frequencies  that  were 
excited,  and  look  at  the  level  of  energy  contained  in 
that  frequency. 

Close  examination  of  the  actual  versus  the  predicted 
plots  show  frequency  content  that  are  almost  identical. 
The  actual  plots  have  frequency  content  not  seen  in  the 
predicted,  but  these  have  to  do  with  modal  participation, 
which  is  not  considered  in  this  portion  of  the  study. 
Frequency  content  is  summarized  in  Table  5-6,  in  the 
concluding  section  for  the  dual  rail  sleds. 

Power  levels  for  the  frequencies  fall  into  the  same 
order  of  magnitude  range  for  all  rail  functions. 
Examining  PSD's  for  the  individual  rail  functions  again 
shows  the  rail  roughness  break  at  0.020  inches. 
Roughnesses  greater  than  this  have  frequency  power  levels 
in  excess  of  the  measured,  and  roughnesses  less  than 
0.020  inches  have  power  levels  less  than  measured. 

Examination  of  the  frequencies  indicated  some 
unexpected  occurrences.  These  were  exhibited  by  all 
rails  unless  noted. 

The  frequency  centered  a^^und  J  2  hertz  on  the  data 
sample  sets  had  to  do  with  a  phenomenon  associated  with 
the  survey.  For  this  velocity  of  s  i  m  u  1  it,  i  n  n  ,  32  hertz 


sample  sets 

h  a  d 

to  do 

the  survey. 

For 

this 

corresponds 

t  ^  a 

track 

distance  for  the  north  and  south  15000  feet  was  50  feet. 
The  long  period  surveys  were  accomplished  in  100-foot 
intervals,  ''nly  repeating  the  last  50-foot  segment  of  one 
section  at  the  beginning  of  the  next  100-foot  section. 
This  resulted  in  a  looping  phenomenon  of  50  feet. 

The  frequency  spike  at  4  2  hertz  corresponds  to  a 
distance  of  39  feet,  or  the  rail  length.  Mixon  noted 
that  for  track  lengths  in  excess  of  30  feet  (and  less 
than  50  feet),  the  north  15000  was  smoother  (had  less 
power)  than  the  south  end.  This  is  confirmed  in  the 
simulations.  Energy  content  in  this  range  for  this 
velocity  (  i  .  e .  ,  32  to  46  hertz)  is  less  for  the  north  end 
of  track  than  for  the  south.  However,  for  all  other 
frequencies,  the  north  end  has  larger  energy  content  than 
the  south. 

The  spike  at  eight  hertz  for  TS  34000,  both  top  and 
bottom  surfaces,  corresponds  to  the  ionpin^  of  the  data 
period  (208  feet)  for  this  velocity.  For  :  ne  sled  to  t>e 
traveling  at  3800  feet  per  second  (  M =  3  .  5 )  ,  the  first 
fundamental  sled  suspension  mode  could  be  excited  by  this 
fundamental  looping  distance. 

5.2.3  RIGID  BODY  CONCLUSIONS  AND  SUMMARY 

This  section  concludes  the  rigid  body  simulation  an) 
analysis  phase  of  the  dual  -ail  sled  section.  it  will 
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address  all  questions  posed  in  the  rail  roughness  study 
directive,  from  a  rigid  body  point  of  view,  that  had  been 
left  unanswered. 

5.2.3. 1  CONCLUSIONS 

5. 2. 3. 1.1  PERIODICITY  VERSUS  NON-PERIODICITY 

By  examining  TS  44000  top  and  the  sf'uth  1  5000  top 

rail  function,  the  statistical  effect  of  periodicity  can 

be  examined.  Each  of  these  profiles  had  the  same  rank 

rail  roughness,  S  =0.012,  but  only  one  had  a  periodic 

y  *  * 

content.  Notice  from  the  slipper  beam  results,  the  same 
roughness  generated  approximately  the  same  standard 
deviation.  However,  the  periodic  rail  underpredicted  the 
mean  by  70  percent  m  o  r  e  than  the  nonperiodic  rail 
( u n de r pr ed  i  c t i o n s  of  15%  versus  24%)  for  the  front  beam. 
This  might  be  explained  by  the  periodic  rail  being  in 
phase  with  the  fundamental  sled  frequency  and  the  slipper 
just  "skimming"  the  top  of  the  smooth  period  rail  and  not 
being  excited.  The  sled  randomly  hits  the  rail  and 
generates  a  sufficient  quantity  of  larger  spikes  to  move 
the  mean  for  the  nonperiodic  rail. 

Looking  at  TS  J4000  top  and  bottom,  it  is  notice, able 
that  rough  periodicity  strongly  affects  the  mean  and 
standard  deviations.  These  profiles  overpredict  these 
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quantities  by  large  amounts.  From  the  frequency  and 
energy  plots  it  is  seen  that  the  fundamental  sled 
frequencies  contain  more  energy  for  rough  period  rails 
than  for  the  nonperiodic  rail  forcing  function. 

5.2.3 .1.2  RAIL  GRINDING 

Rail  grinding  effectiveness  can  be  quantitatively 
considered  by  examining  the  north  and  south  15000  feet  of 
track.  The  south  15000  is  smoother  than  the  north  due  to 
a  longer  period  of  time  to  work  at  smoothing  the  rail. 
The  simulated  results  indicate  the  same  thing.  The  south 
15000,  both  top  and  bottom,  tends  to  underestimate  the 
average  deviation  more  than  the  north  15000.  However, 
the  south  15000  has  a  larger  standard  deviation  than  the 
n  r  r  t  h  end. 

One  hypothesis  is  that  long  period  frequency  content 
(rail  lengths  of  30  to  50  feet)  controls  the  standard 
deviation  at  this  velocity,  but  the  overall  roughness 
controls  the  mean.  This  rail  length  corresponds  to  sled 
frequencies  in  the  30  to  45  hertz  range.  This  range 
does  not  contain  the  fundamental  frequency,  but  is 
typically  in  the  range  that  adds  the  largest  secondary 
effects,  that  is  to  the  standard  deviation.  For  the 
north  15000  this  range  of  track  is  smoother  than  the 
south  end,  and  hence  the  lower  standard  deviation.  The 
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overall  roughness  for  the  north  is  greater  than  the  south 
and  this  explains  the  higher  average  forces. 

It  is  apparent  that  rail  grinding  is  effective  in 
lowering  the  long  period  '■ail  effects  that  contribute  to 
the  dynamics  of  sled  forces.  This  dynamic  force 
reduction  in  turn  reduces  the  fatigue  aspects  of  the  sled 
design  process. 

5.2.3. 1.3  FINAL  RAIL  PROFILE  DESIGNATION 

Based  on  earlier  work  for  SLEDYNE  and  Mixon^’^  using 
monorails,  it  was  shown  that  modal  participation  absorbed 
about  thirty  (30)  percent  of  the  predicted  forces.  Using 
this  fact,  an  ideal  rail  would  be  one  that  overpredicted 
forces  on  the  average  of  less  than  30  percent,  with  a 
small  standard  deviation.  This  c r i t e r i o n  is  met  by  all 
but  one  of  the  r a  i  1  functions.  This  conclusion  is  based 
on  monorail  modal  participation  in  regard  to  "stiff" 
slippers.  The  validity  of  this  assumption  is  suspect  for 
the  much  softer  slippers  and/or  suspension  systems  of 
dual  rail  sleds. 

5.2. 3.2  RIGID  BODY  SUMMARY 

Using  all  data  examined  thus  far,  it  is  inherent 
that  a  rail  roughness  greater  than  0,02)  inches  tents  t  ' 
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overpredict  rigid  body  forces.  This  roughness 

Q 

corresponds  to  the  break  point  defined  by  Mixon  between 
the  roughest  rail  class  (class  A)  and  those  that  are 
smoother  (classes  B  to  D).  It  is  extremely  clear  from 
Table  5-3  this  value  is  the  break  point  for  both  the 

* 

average  and  standard  deviation,  for  both  the  front  and 
aft  beams . 

Based  on  this  rigid  body  work,  by  keeping  overall 
rail  roughness  below  0.020  inches,  the  peak  forces 
induced  into  the  slippers  can  be  r  e  p  u  c  e  d 
s  i  g n  i  f i c a n 1 1 y-- 1 h a t  is,  by  39  percent  from  the  roughest 
rail  (0.048  inches)  to  one  with  a  roughness  of  0.016. 
Across  the  spectrum  examined,  a  one  percent  reduction  in 
peak  force  appears  possible  for  a  0.001  inch  reduction  in 
total  roughness. 

Rail  smoothing  is  also  effective  from  a  localized 
roughness  standpoint.  Peak  forces  are  not  reduced 
drastically  (about  4  percent)  but  the  amplitude  swing  of 
the  dynamic  forces  can  be  reduced  by  as  much  as  17 
percent.  This  is  the  range  where  rail  grinding  appears 
to  be  particulary  effective. 

6.3  MODAL  SIMULATION 

In  order  to  accurately  simulate  the  >-eal  life 


performance  of  the  sled,  the  modal  participation  will 


have  to  be  superimposed  on  top  of  the  rigid  body 


response.  The  modal  participation  has  been  limited  to 
the  structural  modes  of  the  sled  not  directly  related  to 
the  slipper  suspension  system.  These  are  hereafter 
referred  to  as  the  secondary  structural  modes.  For  the 
Ramjet  sled  there  are  five  modes  identified,  and  these 
are  given  in  Table  5-6  in  the  concluding  section  of  this 
chapter  . 

Once  the  secondary  modal  parameters  had  been 
defined,  the  simulation  consisted  of  repeating  everything 
done  under  the  rigid  body  analysis.  In  other  worjs,  the 
eleven  valid  lift  conditions  were  repeated  for  a  track 
distance  of  4000  feet  for  a  constant  velocity  associated 
with  that  lift  condition.  This  was  done  for  each  rail 
roughness  function  for  a  total  of  99  simulations. 

Each  rail  roughness  simulation,  consisting  of  the 
eleven  valid  trajectory  conditions,  was  statistically 
averaged  to  determine  the  effect  of  the  modal 
p a r t i c i p a  i  o n  .  The  results  of  this  analysis  are  presented 
as  an  average  percentage  deviation  between  the  actual 
forces  (see  Table  5-2  for  the  actual  forces)  and  the 
predicted  values.  Also,  the  standard  deviation  about  the 
mean  deviation  is  determined. 

The  effect  of  modal  participation  was  alluded  to  in 
the  'inclusion  of  the  rigid  bn  j  y  work.  There  it  wan 


stated  that  for  one  given  condition  for  monorails,  using 


six  'indes  at  five  percent  damping,  a  reduction  in  the 
1 « ■  v  i  a  t  l  r  n  r.  f  about  30  percent  can  be  expected.  This 
amounted  to  reducing  the  overprediction  from  20  to  40 
p  f  >■  e  e  n  t  to  an  error  of  +  10  percent.  However,  it  is 
unknown  whether  this  is  applicable  to  dual  rail  sleds. 

‘,.3.1  DATA  ANALYS  1. 3 --AMPLITUDE  DOMAIN 

5 . 3 . 1 . 1  PEAK  FORCES 

The  simulation  consisted  of  99  computer  runs  of  the 
valid  lift  conditions  for  each  rail  roughness  function. 
T  :i  e  analysis  contained  the  participation  of  five 
structural  modes  superimposed  on  the  rigid  body  work  done 
previously.  Results  of  the  peak  force  predictions  are 
s  :i  o  w  n  plotted  in  Figures  5.18  through  5.21  for  slippers 
on.'  through  foUr,  respectively.  The  solid  line 
represents  the  actual  peak  forces  for  that  condition  from 
the  1 9  Y  —  C 1  mission.  The  interesting  thing  he-e  is  that 
only  one  rail  function  appears  to  consistantly  exceed  the 
measured  values.  This  is  TS  3^000  bottom  and  is  the  same 
function  hypothesized  in  the  final  rail  profile 
designation  section  for  rigid  body  work  that  would  fail. 

The  statistical  analysis  is  shown  plotted  in  Figure 
5.22  and  summarized  in  Table  5-5  for  the  slipper  beams. 
Tne  " x "  denotes  where  the  average  percentage  deviation 
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between  predicted  and  actual  forces  occurs.  The  heavy 
arrow  represents  the  _+  1  standard  deviation  about  the 
mean. 

Looking  at  the  table,  it  is  seen  that  TS34000  bottom 
severely  overpredicts  as  expected.  TS34000  top  also 
overpredicts  but  only  slightly.  Further  examination  of 
the  table  reveals  a  break  point  in  the  rank  roughness 
predictions  again  occurring  at  0.020  inches.  This  break 
point  holds  for  the  aft  beam  as  well  as  the  front. 

Comparison  of  Tables  5-3  and  5-5  will  show  the 
effect  of  modal  participation.  It  is  readily  seen  that 
the  modal  participation  does  absorb  energy  but  not  as 
much  as  originally  estimated  based  on  the  monorail  work. 
It  was  hypothesized  the  modes  would  absorb  up  to  30 
percent  of  the  forces;  however,  the  range  of  force 
reduction  is  from  5  to  20  percent  for  all  conditions 
examined.  The  standard  deviation  actually  increases  for 
the  modal  simulations  for  all  cases  except  the  smoothest 
rails,  and  then  it  remains 

5. 3.  1.2  PEAK  FORCES  VERSUS  TIME 

As  with  the  rigid  body  work,  the  force-time  history 
for  slippers  two  and  four  was  plotted  for  comparison  to 
the  real  measured  values.  Figures  5 . 2  j  a  and  5.24a  show 
th*1  actual  force-time  history  for  t  h  v  s  o  t  w  o  s  1  1  p  p  «.•  >-  s  . 
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Figure  5.24a  SLIPPER  FORCE  RIGHT  REAR  VS  TIME 
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Figures  5.23b  and  5.24b  contain  the  actual  simulation 
histories  for  the  baseline  rail  function.  The  remaining 
eight  rail  function  histories  are  contained  in  Appendix 
C,  Figures  C.17  through  C.32. 

From  the  actual  data,  slipper  two  forces  stayed  in 
an  envelope  between  6  and  41  kips  and  were  centered  about 
20  kips.  The  plot  also  indicates  a  low  frequency  content 
that  the  data  appears  to  be  modulating  about.  Looking  at 
the  modal  plots  for  all  nine  functions  for  this  slipper 
indicated  the  same  centering  around  20  kips  but  with 
envelopes  of  differing  sizes.  Generalizing,  the  rougher 
the  rail  the  larger  the  envelope  span.  Also,  all 
functions  exhibited  a  low  frequency  content,  yet  the 
smooth  rails  with  periodic  content  (i.e.,  TS44000  top  and 
bottom)  were  the  most  pronounced. 

Examining  the  aft  slipper  force  data  for  the  actual 
run  indicated  forces  fluctuating  by  _+  6  kips.  All 
forcing  functions  exhibited  a  fluctuating  component  of  up 
to  +  6  kips;  however,  they  are  centered  about  a  positive 
one  kip  instead  of  zero  as  for  the  actual  run  data.  The 
rail  functions  that  did  not  contain  periodic  forcing 
content  were  extremely  one-sided,  that  is,  they  ranged 
from  zero  kips  to  a  maximum  value  of  about  3.5  kips.  The 
periodic  content  rails  and  the  rough  rail  were  two-sided 
functions  that  resulted  in  a  low  frequency  content  that 


modulated  about  one  kip. 


Peak  spike  forces  all  were  under  the  actual  measured 
force  for  both  slippers  under  this  simulation  case.  The 
front  slipper  operating  envelope  stayed  within  +  6  kips 
of  the  quasi-steady-state  force  for  the  smooth  rail 
functions.  The  rough  rail  functions  (those  in  excess  of 
0.020  inches)  all  exceeded  this  dynamic  envelope  by  50 
percent  or  more. 

From  the  rigid  body  dual  rail  work,  it  was  noted 
that  the  south  15000  feet  of  track  exhibited  a  pronounced 
rough  spot  in  the  data  that  caused  distinct  force  spikes 
to  occur.  The  modal  simulation  also  exhibited  this  spike 
and  its  cycling  just  as  the  rigid  body  did,  only  it  was 
not  as  pronounced.  The  modal  participation  absorbed 
approximately  a  third  of  the  magnitude  of  this  spike. 
That  is,  for  the  rigid  body  simulation  these  rough  spots 
created  a  nine-kip  dynamic  spike,  whereas  with  modal 
participation,  this  was  only  a  six-kip  spike,  or  a  third 
less  than  the  rigid  body  case. 

5. 3. 1.3  FORCE  HISTOGRAM 

The  rigid  body  dual  ra i 1  simulation  concluded  t  n  a  t 
peak  forces  were  Rayleigh  distributed,  and  the  extent  of 
deviation  from  the  Rayleigh  trend  is  strongly  dependent 
on  the  qu a s i - s t e a d y - s t a t e  forces.  The  results  for  the 
modal  participation  ( 3ee  Figure  5.25)  exhibits  the  same 
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tendency  toward  a  Rayleigh  distribution.  However,  at  the 
expected  maximum  Rayleigh  occurrences  the  simulation 
occurrences  are  only  about  67  percent  of  what  they  should 
be.  This  is  a  result  of  a  small  standard  deviation 
forcing  a  large  expected  quantity  into  the  peak 
occurrence  range.  The  reason  fo<-  the  small  standard 
deviation  is  the  reduction  in  the  dynamic  forces  due  to 
modal  energy  absorption.  All  simulations  indicated  the 
same  phenomena.  The  extent  of  the  error  in  the  maximum 
number  of  occurrences  is  related  to  the  range  of  the 
standard  deviation.  The  smaller  it  is,  the  more  extreme 
the  difference  in  the  number  of  occurrences. 

It  is  concluded  that  the  peak  force  simulation, 
including  the  effects  of  modal  participation,  does  in 
fact  follow  a  Rayleigh  distribution.  The  amount  of 
deviation  is  a  direct  result  of  the  quasi-steady-state 
forces  and  the  dynamic  force  envelope. 

5.3.2  DATA  ANALYSIS--FREQUENCY  DOMAIN 

The  highest  frequency  content  obtainable  was  275 
hertz  from  the  rigid  body  simulation.  This  limit  had  to 
do  with  the  timer  integration  step  sizes  which  are 
related  to  the  highest  frequency  content  of  the 
simulation.  The  largest  frequency,  for  modal 
participation,  deals  with  the  highest  mode  used,  which  is 
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same  order  of  magnitude  as  the  measured  data  for  all  rail 
functions.  All  rail  functions  except  the  two  roughest 
(TS3^000  top  and  bottom)  fall  below  the  measured  levels, 
and  these  two  rough  rails  exceed  it.  At  the  rank 
roughness  of  0.020  inches,  again  a  distinct  separation  of 
power  levels  occurs.  Roughnesses  less  than  this  value 
have  power  levels  distinctly  below  those  of  the 
roughnesses  in  excess  of  0.020  inches. 

The  same  unique  frequency  phenomena  are  present  here 
as  in  the  rigid  body  simulations.  However,  the 
facsimiles  of  the  data  survey  are  amplified  in  this  case. 
The  frequency  at  32  hertz  is  the  50-foot  survey  data 
sampling.  It  also  coincides  with  a  sled  mode  at  34 
hertz.  This  coupling  induces  a  large  energy  spike  at 
this  frequency,  in  fact,  it  is  about  300  percent  of  that 
without  the  mode  shape  being  superimposed.  The  frequency 
at  42  hertz  couples  the  physical  rail  lengths  of  39  feet 
to  a  sled  frequency.  In  this  case  the  sled's  frequency 
mode  shape  does  not  contribute  much  to  the  power  level  as 
the  power  is  approximately  equal  to  that  without  the  node 
shape  . 

The  spike  at  4  hertz  on  the  S  1  5  0  0  0  top  data  is  the 
data  set  looping  end  to  end.  This  frequency  has  energy 
content  so  low  in  the  rigid  body  simulation  that  it  is 
apparently  nonexistent.  However,  this  spike  is  the 
predominant  ope  in  the  vertical  direction,  and  no 
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plausible  explanation  is  available  outside  of  possible 
modal  coupling  . 

The  spike  at  8  hertz  on  the  34  0  0  0  top  rough  periodic 
rail  is  the  looping  of  the  rail  period.  This  data 
facsimile  was  noted  on  the  rigid  body  work  and  is  again 
tne  controlling  energy  spike.  It  is  the  same  magnitude 
as  is  the  rigid  body  spike. 

5.3.3  MODAL  CONCLUSIONS  AND  SUMMARY 

This  section  concludes  the  modal  simulation  analysis 
for  a  dual  rail  sled.  It  will  address  all  questions 
posed  in  the  rail  roughness  study  directive  from  a  modal 
participation  point  of  view. 

5. 3.3.1  PERIODICITY  VERSUS  N ON- P E R 1 0 D I C I T Y 

As  for  the  rigid  body  simulation,  the  smooth 
periodic  rail  (TS44000  top)  generated  the  same  standard 
deviation  as  the  smooth  nonperiodic  rail  (TS  SI  5000  top); 
and  this  value  is  25  percent  less  than  the  rigid  body 
value.  Unlike  the  rigid  body  work,  the  periodic  rail 
underpredicted  by  10  percent  less  than  the  nonperiodic 
rail  (underpredictions  of  28  percent  to  31  percent)  for 
the  front  beam.  This  implies  that  the  modal  content  is 
moderately  excited  by  the  periodic  rail  content.  This 
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implication  is  further  substantiated  when  examining  the 
effect  of  modal  energy  absorption  between  the  rigid  body 
and  modal  cases.  The  modal  content  reduced  the  average 
deviation  underprediction  by  16  percent  for  the  periodic 
rail.  Comparably  the  nonperiodic  rail  modal  content 
reduced  its  average  underprediction  by  107  percent. 

The  comparison  between  rough  periodicity  and  smooth 
periodicity  is  similar  to  the  rigid  body  case  .  Rough 
periodicity  strongly  affects  the  mean  and  standard 
deviations.  However,  the  reduction  in  the  average  force 
predictions  is  significantly  more  for  the  rough  periodic 
rail  than  for  the  smooth  periodic  rail,  i.e.,  87  percent 
compared  to  16  percent.  Conversely,  the  rough  periodic 
rail  in  conjunction  with  the  modal  content  increases  the 
standard  deviation  by  20  percent. 

5. 3. 3.2  RAIL  GRINDING 

The  effectiveness  of  rail  grinding  can  be 
quantitatively  examined  by  looking  at  the  north  versus 
south  18000  feet  of  track.  The  south  end  of  the  track 
underpredicts  the  mean  when  compared  to  the  north,  as  is 
expected.  What  was  not  expected  was  the  standard 
deviation  to  be  smaller  for  the  south  end.  The 
interesting  phenomenon  here  is  that  the  content  for  the 
fundamental  heave  frequency  for  the  north  end  increased 
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by  67  percent  from  the  rigid  body  situation.  This 
increase  is  apparently  due  to  the  first  roll  structural 
mode  being  almost  superimposed  on  the  fundamental 
vertical  heave.  This  coupling  of  modal  roll  and  rigid 
body  heave  is  causing  a  large  dynamic  force  excursion  to 
occur.  The  other  secondary  structural  frequencies  have 
energy  content  considerably  less  than  the  south  end,  and 
these  frequencies  lie  in  the  smoother  north  end  energy 
content . 

Rail  grinding  must  still  be  considered  effective  for 
reducing  the  dynamic  force  range  as  long  as  the  forcing 
function  frequencies  lie  in  the  secondary  structural 
frequency  range.  This  statement  must  be  caveated  to 
account  for  structural  modes  acting  in  the  direction  and 
frequency  of  a  rigid  body  mode,  i.e.,  roll  structural 
mode  and  a  rigid  body  heave  coupling. 

5. 3. 3. 3  FINAL  RAIL  PROFILE  DESIGNATION 

As  hypothesized  in  the  rigid  body  section,  the  modal 
content  does  reduce  the  peak  forces  occurring  during  the 
simulation  but  not  up  to  j  0  percent  as  thought.  Based  on 
the  original  asssumption,  it  was  thought  all  rail 
functions  except  TS 34000  bottom  might  provide  adequate 
simulation  results.  This  assumption  was  subsequently 


verified  as  being  accurate,  with  T  S  3  11  0  0  0  top  only 


marginally  o v e r p r e d i c t i n g 


To  obtain  the  best  average 


prediction  with  the  smallest  standard  deviation,  the 
forcing  functions  to  use  should  be  the  north  1 3  0  0  0  bottom 
or  the  baseline  function.  Preference  shoal!  be  given  to 
the  north  end  bottom  funtion  . 

5 . j . j . 4  MODAL  SUMMARY 

From  all  data  reviewed  thus  far,  it  is  i  nr.  n  •  ‘  • ,  t 
rail  roughnesses  in  excess  of  0.020  inches  t e  n u 
overpredict  forces  even  with  the  modal  energy  absorption. 
This  break  point  holds  for  both  the  front  as  well  a  s  t  n  e 
aft  slipper  for  both  the  average  deviation  as  well  .s  ten- 
standard  deviation. 

The  modal  content  affects  peak  forces  in  varying 
degrees  but  in  all  cases  it  reduces  the  average  deviation 
from  the  rigid  body  case.  This  reduction  is  in  the  r 
of  ‘3  to  20  percent,  not  j  0  percent  as  originally 
hypothesized  . 

Based  on  the  modal  simulations,  the  peak  s  1  i  p  p •  • 
forces  can  be  reduced  significantly  by  reducing  t  r  -  r  j  r  ? 
rail  roughness  .  The  amount  of  reduction  is  s  l  n  i  .  - 
that  seen  for  the  rigid  body  case.  A  o n  e  (  1  )  per 
r  e  d  u  <•  t  i  o  n  in  peak  forces  is  obtainable  [or-  .  r. 
inch  reduction  in  rail  roughness. 

Frcm  a  localized  r  i  1  r  ^  u  •  p  p  --  i  n  •  ,  ■  ,  i  . 


smoothing  is  effective  in  reducing  slipper  forces  in  the 
dynamic  envelope  by  one-third  when  the  modal 
participation  is  utilized.  This  is  true  unless  a  mode 
shape  happens  to  coincide  with  a  rigid  body  frequency, 
then  this  dynamic  force  reduction  does  not  occur. 

b.4  DUAL  HAIL  NUMMARY  AND  CUNCLULIONL 

In  an  effort  to  verify  the  DALTArtR  program  for  dual 
rail  sleds,  an  extensive  number  of  computer  i-uns  were 
performed.  These  simulations  were  all  performed  using 
the  constant  velocity  option  lor  a  single  sled.  The 
variable  speed  simulation  has  Deep  performed  and  is 
discussed  now. 

•J  .  n  .  1  DUAL  KAIL  VARIABLE  VELOCITY  0  I  M  J  L  A  T  :  ... 

In  order  to  perform  variable  velocity  simulations,  a 
trajectory  analysis  has  to  exist.  It  i  s  file  is  one 
created  by  the  in-house  tracx  trajectory  program  called 
i  .1 1,  ED.  This  f  l  1  e  in  us  t  be  fund  e  d  i  t  e  .*  to  r  c  x  o  v  j  i  1 
li  i>  ,i  d  e  r  information.  The  resultant  file  contains 
,n formation  regarding  the  s  L  e a  '  s  d  o  w  n  t  r  a  o  k  velocity, 
t  r  a  c  k  station,  and  acceleration.  It  also  generates  the 
aerodynamic  drag,  t  h  e  thrust,  and  sled  weight  vectors 
acting  on  the  sled  at  a  given  instant  of  time  du^i ng  the 


trajectory.  This  trajectory  analysis  was  performed  to 
duplicate  the  velocity  profile  shown  in  Figure  j . j . 
Using  these  variable  inputs  from  TSLED,  the  program  will 
linearly  interpolate  the  lift  and  its  center  fr^n  user- 
supplied  tables.  The  simulation  will  t  ri  e  n  progress  •:  s 
for  the  constant  velocity,  using  only  t  he  variable 
parameters  instead  of  the  constant  ones. 

The  actual  simulation  for  the  variable  velocity  was 
terminated  when  a  speed  of'  1  5  1  7  feet  per  second  under 
thrust  was  achieved.  This  velocity  is  the  time  when  the 
on  board  propulsion  ignited  and  the  sled  mass  started 
changing.  The  results  of  this  simulation  are  shown 
plotted  in  Figures  5.28a  for  the  rigid  body  case  and  in 
5.28b  for  the  modal  case.  The  individual  slippers  are 
plotted  as  the  percentage  deviation  in  peak  force  at  the 
given  terminal  velocity  of  the  simulation. 

Notice  from  the  rigid  body  simulation  that  each 
slipper  except  for  number  two  tends  to  underpredict  by  10 
to  20  percent  in  all  cases  for  raij  roughnesses  less  than 
0.020  inches.  Slipper  two  overpredicts  by  about  10 
percent  on  the  .  ■  v  e  r  a  g  e  for  this  same  rail  classification. 
Koughnesses  above  this  defined  break  point  overpredict 
more  severely.  The  modal  simulation  results  (Figure 
5.23b)  mirror  those  of  the  rigid  body,  only  the 
overpre dictions  are  about  ,0  percent  larger  than  for  the 


rigid  body. 
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Examination  of  the  data  output  for  the  modal 
simulations  shows  all  slippers  peaking  at  times  that 
correspond  with  the  same  Mach  number  as  they  did  for  the 
rigid  body  simulation.  Since  the  rigid  body  contribution 
has  to  be  constant  between  simulation  cases,  the 
differences  must  be  due  to  the  modal  properties  being 


excited. 
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feet  must  be  present  to  excite  these  frequencies  at  this 
velocity.  This  encompasses  the  b  0  -  foot  survey  looping 
distance  common  to  all  rail  profiles.  It  is  evident  that 
this  survey  facsimile  is  exciting  the  predominant 
secondary  modes  and  is  adding  approximately  30  percent  to 
the  slipper  forces.  The  rejsop  this  did  not  occur  in  the 
constant  velocity  simulation  is  that  this  exact  velocity 
was  never  encountered. 


5 .  R  .  CONCLUSIONS 

It  is  apparent  that  DA  STARR  accurately  predicted  t h e 
peak  slipper  forces  for  a  dual  rail  sled  from  all  dat,. 
examined  thus  far-.  Depending  oP  the  r  a  i  1  function  used, 
velocity  option,  or  modal  configuration  the  results  are 
extremely  good. 

It  is  apparent  that  the  constant  velocity  function 
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for  a  single  sled  will  suffice  for  a  given  anlaysis 
s  i  t  u  j  t  1  r>  n  .  However,  t  o  use  this  option,  the  engineer 
will  have  to  define  the  actual  aerodynamic  loads  and  the 
transmitted  thrust  acting  on  that  sled.  This  option  does 
not  allow  for  the  transmitted  thrust  transients  to  be 
considered.  Trie  multiple  sled  option  with  variable 
velocity  should  be  used  to  counter  this. 

In  this  chapter,  it  was  proven  that  the  peak  forces 
for  a  dual  rail  sled,  either  as  a  rigid  body  or  coupled 
with  the  modal  participation,  are  indeed  a  Kayleigh 
distribution  as  Mixon  suggested.  However,  this 
distribution  is  strongly  affected  by  the  mean,  or  quasi- 
steady-state  forces,  and  tne  standard  deviation.  That 
is,  the  Gaussian  contribution  is  greatly  enhanced  by  the 
lower  the  quasi-steady-state  force  component.  Similarly 
the  skewness  of  the  distribution  is  dependent  on  the 
standard  deviation.  The  smaller  this  deviation,  the  more 
skewed  toward  the  mean  is  the  distribution. 

The  r  u 1  1  function  defined  by  the  baseline  rail  of 
the  bottom  of  north  end  of  the  track  yields  the  best 
slipper  estimates  from  a  statistical  standpoint.  These 
functions  apply  to  both  rigid  body  and  modal  analyses  for 
either  constant  o r  variable  velocities. 

Hail  grinning  is  indeed  effective  but  not  as 
o  r i s i n  a i  1  y  thought.  It  was  supposed  to  smooth  the  r  a i 1 
lor  pencils  of  length  less  tn.in  b  2  inches.  Instead  it 


smooths  the  long  period,  i.e.,  30  to  50  feet  in  length. 
These  long  period  effects  are  the  ones  that  typically 
affect  the  secondary  structural  modes  of  the  dual  rail 
systems.  This  smoothing  has  the  net  effect  of  reducing 
the  dynamic  oscillation  of  the  forces  about  the  quasi¬ 
steady-state  mean  force.  This  statement  has  to  ue 
flagged  because  rail  grinding  does  not  affect  the  dynamic 
envelope  if  a  mode  shape  is  superimposed  on  a  rigid  body 
mode,  i.e.,  first  structural  roll  mode  coupled  with  a 
rigid  body  vertical  translation. 

In  general,  a  one  percent  reduction  in  peak  forces 
is  possible  for  each  0.001  inch  reduction  in  rank  ''ill 
roughness.  In  order  to  significantly  affect  slipper 
forces,  the  rank  roughness  needs  to  be  held  at  or  below 
0.020  inches.  Going  below  this  demarcation  does  not 
reduce  si ippc  forces  as  drastically  as  just  attaining 
this  level. 

Verification  of  localized  rough  spots  op  the  track 
has  been  achieved.  The  effect  of  reducing  the  localized 
high  spot  can  now  be  studied  to  obtain  the  optimum  level 
to  which  grinding  should  be  performed. 

One  drawback  to  this  analysis  is  in  the  definition 
of  the  forcing  function.  Present  survey  techniques 
introduce  "spikes"  into  the  data  that  in  some  instance.-, 
adversely  affect  the  results.  These  survey  techniques 
must  be  improved  to  better  future  analyses.  However,  the 


present  definitions  are  more  realistic  than  trying  to 


assume  ,i  forcing  function  ,  as  has  been  done  previously. 

Verification  of  the  lateral  results  has  not  been 
attempted  since  no  lateral  force  data  exist  for  dual  rail 
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-  .  '  r  i  n  ♦  ,i  n  5  depend  i  p,;  o  n  direction  o  f  applied 

:  i  ’  . 

W  ■•••  n  ■••-'I.."'  ~  I  ]  ,bc  aPd  May  o  f  1  ,1  !  '<  ,  a  SClvS  o  f 

•-  r  ,, .  -  k  runs  war-.'  p  a  r  f  r  tr,  e  d  us  i  ri£  the  instrumented  mcncr  a  1 1 

it--  ill.  "  t  j  L  c  i  l  v  e  s  r  f  -  n  i  s  t  ••  s  t  r  i  e  s  a  r  e 

s  .  i -  1  If  1  b  ••  1  w  . 

1  ‘  'i  -  ■  U  r  'be  V  r  t  1  C  u  .  ,  "a  1  .  i.  - 1  t  ^  y  ,  :  n  ,1  ^  U  ..  S  1  — 

s  '  ■  . :  1  y  -  s  t  .  1 1  -  slipper  .  a  j  •  a  a  r  ,j  ..  p  i  R  .  “-p-rai  1  sled 

1  a  r  velocities  up  to  ii  0  d  0  feet  per  second. 

>■  f  :  i  i  n  e  t  h  e  <•  f  f  e  c  t  of  ,  n  1  s  o  ;  ,1 1  slipper  o  p 

ii.l 


FIGURE  h 


6.1.2  SLED  FORCES 

Based  o  n  figures  presented  u  y  M  1  x  "  r  ,  p  u  .  ,  . . 

information  was  available  t  n  duplicate  bis  efforts  orj 
see  if  better  correlation  would  be  ,i  v  j  i  1  j  b  1  e  b  >  .  •  g 

DASTAKR  .  The  critical  factor  at  this  stage  is,  .<  s 

mentioned  earlier,  to  accurately  :  ♦.?  f  i  r>  «.*  the  ..ft  for  , 
and  the  location  of  its  center  of  pressure. 

The  first  step  taken  to  e  u  1  c  u  1  at  the  ;  u  .  s  i  -  s  t  c  .in¬ 
state  lift  force  is  t  o  duplicate  t  h  e  t  r  ,  e  c  t  o  r  y  <-  f  t  he 

mission.  From  this,  it  will  be  possible  to  determine  the 
thrust,  dr  ug  ,  and  inertial  loadings  acting  o  r  trie  sled  .i  s 
quasi-steady-state.  Using  t  n  *.*  s  <•  v  a  i  u  e  s  ,  along  wit  n 
M  i  x  o  p  '  s  given  quasi-steuay  slipper  r  ,,  c  t  i  o  n  5  ,  t  n  e  lift 
force  can  be  found. 

Shown  in  Figure  a  .  2  a  is  the  a  e  t  u  1  *•  u  r  i  r  .i  j  e  c  t  °  r  , 

used  by  Mixon.  Figure  6.2  b  shows  the  generated  p  r  ~  •  i  ; 

that  attempts  to  mirror  the  profile  -  f  " ;  x  o  r  ,  :  :  ■  ,  • 

and  velocity  maximum  b  e  t  w  e  e  r  the  two  are  con  s  i  s  t  e  r-  *.  ; 

however,  the  time  of  p  e  a  k  velocity  o  ■  :J  r  r  r 

different.  This  is  due  to  j  faster  burning  r  ■"  g  -■  1  1  •  t 

n  o  w  u  s  e  d  i  r  t  h  a  t  r  o  c  w  <  r.  :  o  r.  o  r  .  1  h  .  -  -  :  r  ,  .  ■-  ,  a  .  a  ■ 

still  the  same  so  the  n  e  t  v  t;  1  "  t  i  t  y  e  f  f  e  c  t  is  t  h  e  s  i  •  ; 

the  tine  of  pe.iK  velocity  ;  a  :■  t  jr.  .  „r..r  ;  r  :  ■  .  r  . 

Shown  :  r>  Figure  ••  .  ..  r  >■  *.  n  e  si  .  p  j  •  ,,  u  a  s  l  -  r  *.  • 

s  t  a  t  e  for  e  u  sod  by  ".yon  i  n  a  i  :■  ,".u!  ,  h  ;  '  r  .  t  -  t  ;  ‘  . 
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the  static  freebody  equations  of  motion  about  the 
s  !  «.* o  ,  the  lift  .1  r'd  its  center  can  be  calculated  using  the 
scaled  qua si-st eudy-st ate  values  o  f  figure  6 .  j  a .  Figure 
t,  .  4  strws  the  resultant  lift  force  that  was  computed  by 
1  ~  i  i  r'  w  i  n  g  the  above  procedures  at  each  quarter  Mac  h 
number  between  1.0  and  the  maximum  velocity. 
,,'jpcr  imposed  in  circles  op  this  f  l  g  u r  .ire  the  estimated 
lift  forces  derived  from  wind  tunnel  testing  on  a  40% 
s  e  a  1  e  model  .  c  '  A  s  a  c  r.  e  c  n.  tne  quasi -  steady  slipper 
ads  were  recomputed  using  the  cululated  lift  force  and 
t. e  ire  shown  in  Figure  6  .  u  .  Trie  shape  and  magnitudes 
r  f  the  forces  match  those  of  Hire's.  Table  6-1  contains 
trie  instantaneous  quasi-st  e  ad  y-state  values  used  in  these 
calculations  f^r  tne  lift  forces. 

r.  .  KluILi  BODY  -  IM'JLATION 

T  h  e  simulation  consisted  of  running  e  a  c  h 
instantaneous  condition  shown  in  Table  6-1,  or  twelve 
t  n  r  ij  t  l  n  g  ,-opii  it  i  o  ns  fro  ~  y  a  cli  1  .  0  t  o  Mach  3  .  6  2  ,  t  h  e 
maximum  velocity.  There  were  also  seven  conditions  under 
c  r  a  ;i  t  from  maximum  velocity  bacK  t  o  Mach  2.0.  It  should 
b  ••  noted  that  the  lift  forces  between  Mach  1.0  and  2.0, 
un  !  ■■  r  •.  1,  rust,  a  r  e  suspect  sin  c  e  t  n  e  y  w  e  r  found  by 
e  x t  r  a  p  o 1  a  t i o n  from  the  wind  tunnel  results. 

T!i"  simulation  amounted  to  running  ail  1  )  conditions 
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for  each  rail  function  at  a  track  length  of  4000  feet  for 


a  given  constant  velocity,  ~r  a  total  of  171  iterations 
for  the  rigid  body  case.  The  subsequent  analysis 
compared  the  predicted  slipper  peak  values  to  th^se 
scaled  from  the  figures  generated  by  Mixon.  Table  6-2 
contains  the  scaled  average  force  values  as  well  as  the 
extremes  encountered  for  that  Mach  number.  Notice  that 
no  sign  is  given  so  only  the  absolute  values  may  be 
compared.  Also,  since  no  distinction  between  thrusting 
or  coasting  conditions  have  been  defined,  it  is  assumed 
that  the  value  given  applies  to  either  condition  at  that 
Mach  number. 

Comparison  of  the  lateral  to  vertical  slipper  peak 
forces  from  Table  6-2  indicate  that  lateral  forces  range 
between  40  and  60  percent  of  the  peak  vertical  frrce .  If 
the  quasi-steady  portion  of  the  peak  vertical  force  is 
removed  from  the  peak,  the  lateral  force  ranges  between 
50  and  70  percent  of  the  dynamic  vertical  force.  Present 
track  design  guidelines  call  for  the  lateral  dynamic 
forces  to  be  60  percent  of  the  vertical  dynamic  force  and 
are  a  direct  result  of  this  comparison. 

As  was  done  for  dual  rdil  sleds,  the  analysis  will 
statistically  average  all  conditions  for  a  given  rail 
function  for  each  slipper.  The  results  of  this  sampling 
are  to  be  presented  as  a  percentage  deviation  between 
actual  and  predicted  forces  and  the  standard  deviation 


TABLE  6-2  SCALED  PEAK  OSCILLOGRAPH 
FORCES  FROM  INSTRUMENTED  MONORAIL 
TEST  RUNS,  OPLAN  44X 


Mach 

Vertical-Front 

Vertical-Aft 

No . 

Ave . 

Max . 

Min . 

Ave . 

Max . 

Min . 

1  .00 

6.00 

6.50 

5.50 

7.40 

8.00 

1.25 

9.15 

10.50 

7.80 

10.00 

11.60 

8.40 

1.50 

9.25 

9.50 

9.00 

10.40 

11.60 

9.20 

1.75 

9.50 

9.50 

9.50 

12.40 

12.40 

12.40 

2.00 

9.38 

10.00 

8.75 

12.60 

13.60 

11.60 

2.25 

10.00 

10.75 

9.25 

11  .60 

11.60 

11.60 

2.50 

12.25 

14.00 

10.50 

11  .60 

12.20 

11  .00 

2.75 

11.75 

12.00 

11  .50 

10.60 

12.80 

8.40 

3.00 

13.50 

14.50 

12.50 

12.00 

14.00 

10.00 

3.25 

14.25 

15.50 

13.00 

15.20 

17.60 

12.80 

3.50 

16.00 

16.25 

15.75 

19.20 

19.40 

19.00 

Mach 

— 

Lateral-Front 

Lateral-Aft 

No . 

Ave . 

Max . 

Min . 

Ave . 

Max . 

Min . 

1  .00 

3.25 

_ _  __  _ 

3.45 

4.30 

'  1 

1.25 

2.83 

6.88 

4.50 

4.80 

1  .50 

4.08 

6.13 

6.40 

6.40 

6.40 

1.75 

4.02 

5.50 

2.00 

5.10 

6.00 

4.20 

2.00 

4.55 

7.00 

2.88 

3.55 

3.80 

3.30 

2.25 

5.14 

8.13 

3.38 

6.45 

6.60 

6.30 

2.50 

5.19 

7.38 

3.75 

4.50 

4.80 

4.20 

2.75 

6.14 

8.75 

4.25 

3.35 

3.70 

3.00 

3.00 

5.60 

9.00 

3.63 

4.90 

6.00 

3.80 

3.25 

6.38 

7.38 

5.38 

3.50 

3.50 

3.50 

3  .50 

6.41 

7.25 

5.38 

4.75 

4.90 

4.60 

about  this  mean. 


6.2.1  DATA  ANALYSIS-AMPLITUDE  DOMAIN 


6. 2. 1.1  PEAK  FORCES 


» 


The  simulation  consisted  of  171  computer  runs  of  1,  n  e 
19  conditions  for  each  of  the  nine  rail  profiles. 
Results  of  tiie  peaK  force  simulations  are  shown  plotted 
in  Figures  6.5  through  6.8  for  the  four  slipper 
conditions.  The  front  vertical  and  lateral  slipper 
forces  are  shown  in  Figure  6.6  and  6.6.  Figures  6.7  and 
6.8  represent  the  aft  slipper  vertical  and  lateral 
forces.  The  solid  line  is  the  average  force  shown  in 
Table  6-2,  and  the  dotted  lines  correspond  to  the 
extremal  values  for  that  given  Mich  number.  It  is 
interesting  to  note  that  the  lateral  force  predictions 
are  equal  to  or  exceed  the  measured  forces,  whereas  the 
best  Mixon  could  do  was  to  predict  approximately  40 
percent  of  the  value.  It  appears  that  for  monorail 
sleds,  track  roughnesses  in  excess  of  0.014  inches  tend 
to  grossly  overpredict  values;  and  roughnesses  less  than 
or  equal  to  this  value  tend  to  estimate  very  closely  the 
results  of  Mixon. 

The  statistical  deviation  analysis  is  represented  in 
Figure  6.9  and  summarized  in  Table  6-4  for  all  slipper 
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conditions.  The  "x"  denotes  where  the  average  deviation 
for  the  19  simulations  occurs.  The  heavy-arrowed 
vertical  line  indicates  the  +_  1  sigma  value  about  the 
mean  . 

6. 2. 1.2  PEAK  FORCES  VERSUS  TIME 

The  only  plot  of  force  versus  time  presented  by 
Mixon  was  a  short  data  clip  on  the  front  slipper  at  an 
average  velocity  of  approximately  3^00  feet  per  second. 
This  corresponds  to  a  Mach  number  of  3  •  0  0  .  The  plotted 
forces  are  the  total  vertical  and  lateral  forces  acting 
on  the  front  slipper  for  a  0.095  second  time  period. 
This  condition  coincides  with  a  data  case  being  simulated 
so  a  direct  comparison  is  available.  Figures  6.10a  and 
6.11a  are  reproduced  copies  of  Mixon's  data  clips  for  the 
front  vertical  and  lateral  slipper  fnrces,respectively 
(page  68  of  reference  6).  Shown  in  Figures  6.10b  and 
6.11b  are  the  simulated  force-time  plots  for  the  baseline 
rail  function  of  these  slipper  forces.  The  remaining 
eight  function  plots  are  shown  in  Appendix  C,  Figures 
C.33  through  C . 4 8 .  The  sign  convention  for  the  simulated 
output  of  DASTARR  is  opposite  that  for  the  actual  data. 

Examination  of  the  actual  run  data  for  the  front 
slipper  vertical  indicates  the  forces  stayed  in  an 
envelope  between  +8  kips  to  a  -12  kips,  and  centered 


about  a  negative  one  or  two  kips.  This  centering  is  to 
be  expected  as  the  quasi-steady  forces  on  this  slipper 
are  approximately  one  kip  negative.  Looking  at  the  plots 
for  the  various  rail  surfaces  for  this  slipper  indicates 
the  same  phenomenon,  a  centering  about  a  negative  one  or 
two  kips  and  envelopes  of  varying  size.  The  envelope 
sizes  were  all  larger  than  the  measured  forces  and  appear 
to  be  a  function  of  rail  roughness.  Of  interest  on 
monorail  sleds  is  the  amount  of  flight  time.  The  amount 
of  flight  time  appears  to  be  related  to  rail  roughness. 
The  smoother  the  rail  the  more  flight  time  involved. 

The  front  slipper  lateral  forces  for  the  actual  run 
data  centered  about  zero  kips  and  ranged  predominantly 
between  plus  or  minus  four  (4)  kips.  This  centering  is 
expected  since  there  is  no  lateral  quasi-steady  forcing 
function  involved.  These  results  should  be  strictly 
dynamic  forces  induced  by  rail  roughness.  Examination  of 
all  simulations  for  the  lateral  forces  indicated  the  same 
thing,  centering  about  zero  kips  and  modulating  an  equal 
magnitude  on  either  side  of  it.  The  interesting  thing 
about  the  lateral  simulation  is  the  ringing  seen  in  the 
output.  The  forces  increase,  peak,  and  decrease  just  as 
mechanical  systems  do  when  approaching  and  passing 
through  a  natural  frequency;  however,  this  ringing  is  not 
seen  in  the  actual  rail  data. 


6 . 2 . 1 . 3  FORCE  HISTOGRAM 


Discussed  in  Chapter  Five  is  the  histogram  feature 
of  DASTARR  as  it  applied  to  dual  rail  sleds.  It  was 
concluded  that  Mixon's  hypothesis  for  dual  rail  sled  peak 
forces  is  true;  peak  forces  do  follow  a  Rayleigh 
distribution.  It  was  noted  in  the  discussion  for  dual 
rail  sleds  that  the  magnitude  of  the  quasi-steady-state 
forces  strongly  influences  the  distribution.  The  larger 
this  value,  the  more  biased  toward  a  Rayleigh 
distribution  is  the  data  sample.  In  contrast,  the  smaller 
the  quasi-steady  component  the  more  closely  the 
distribution  approximates  a  normal  or  Gaussian  function. 
For  the  case  of  this  monorail,  the  data  indicates  a  small 
quasi-steady  value  and  plots  like  a  slightly  skewed 
normal  distribution  as  shown  in  Figure  6.12.  The 
distribution  in  reality  appears  to  be  a  cross  between  a 
normal  and  Rayleigh  distribution.  Mixon  in  his  monorail 
work  did  an  extensive  amount  of  research  in  determining 
how  much  of  a  distribution  is  Rayleigh  content  and  how 
much  is  due  to  a  normal  distribution.  He  did  this 
through  an  error  function  analysis  approach.  He  noted  in 
his  work  that  the  Rayleigh  component  was  dependent  on  the 
bias  introduced  by  the  quasi-steady  loading.  These 
results  indicate  the  same  thing. 


Figure  6.12  Histogram  Distribution 
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6.2.2  DATA  ANALYSIS — FREQUENCY  DOMAIN 


Due  to  the  nature  of  monorail  sleds,  fundamental 
natural  frequencies  generally  occur  in  excess  of  100 
hertz  in  the  vertical  direction.  It  is  not  uncommon  for 
frequencies  up  to  500  hertz  to  cause  structural 
excitement  and  possible  damage  to  the  sled.  There  are  a 
sufficient  number  of  occurrences  for  the  monorail 
simulation  that  it  is  possible  to  analyze  the  system 
response  spectrum  through  2500  hertz. 

An  examination  of  PSD  plots  for  the  given  data  of 
Mixon,  shown  in  Figures  6.13a  and  6.14a  for  the  vertical 
and  lateral  directions  of  the  front  slipper,  shows 
maximum  energy  content  in  a  frequency  range  of  100  to  150 
hertz.  Sharp  rolloff  occurs  above  this  value,  and  noise 
appears  above  600  hertz.  As  explained  above,  this  is  a 
typical  response  spectrum  for  monorail  sleds. 

Figures  6.13b  and  6.14b  show  typical  wide  spectrum 

response  for  the  rigid  body  monorail  sled  simulation 

sequence  for  the  same  given  velocity  condition  discussed 

in  Figures  6.13a  and  6.14a.  Again,  sharp  rolloff  is 

noticeable  above  120  hertz,  and  noise  is  encountered 

above  800  hertz.  Figure  6.13c  and  6.14c  are  blowups  of 

Figures  6.13b  and  6.14b,  respectively,  in  the  0  to  1000 

?  5 

hertz  region,  only  clipping  the  power  axis  at  10  J  and 
g  n 

13.  With  these  figures,  it  is  possible  to  isolate 
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individual  frequencies  that  were  excited,  and  look  at  the 


level  of  energy  contained  in  that  frequency  band. 

Close  examination  of  the  plots  for  the  actual  versus 
predicted  PSDs  indicates  frequency  contents  that  are 
almost  identical,  considering  that  the  actual  data  have 
been  smoothed  by  a  40-hertz  filter.  Due  to  this 
filtering,  only  predominate  power  levels  will  show  and 
will  peak  within  40  hertz  of  the  true  frequency.  Because 
of  this,  many  modes  and  rigid  body  frequencies  may  be 
impossible  to  isolate.  The  measureable  frequency  content 
is  summarized  in  Table  6-6,  in  the  concluding  section  for 
monorail  sleds. 

Power  levels  for  the  frequencies  are  all  in  excess 
of  that  from  the  actual  run  data;  however,  they  are  well 
within  an  order  of  magnitude  difference.  Examination  of 
the  PSD's  again  show  that  rail  roughnesses  in  excess  of 
0.014  inches  have  power  levels  that  greatly  exceed  those 
from  the  run,  while  roughness  levels  less  than  this  value 
have  energy  content  that  is  relatively  close  to  the 
actual  content. 

Reviewing  the  PSD's  indicated  the  following 
phenomena  associated  with  the  frequency  content.  The 
frequency  occurring  at  86  hertz  is  the  physical  rail 
length  of  39  feet.  The  tie  down  spacing  of  52  inches 
occurs  at  a  frequency  of  780  hertz  but  this  spike  is  hard 


to  observe  since  so  little  energy  resides  at  this 


frequency.  The  ringing  phenomenon  seen  in  the  lateral 
direction  occurs  at  about  5b  hertz,  which  is  the 
fundamental  lateral  frequency.  The  frequency  due  to  the 
50- font  data  survey  iteration  at  this  velocity  is  65 
hertz.  It  appears  that  the  ringing  is  being  caused  by  the 
fundamental  lateral  frequency  being  excited  by  the 
forcing  function  frequency.  Each  ringing  envelope 
repeats  at  a  frequency  of  about  7  hertz,  which  is  the 
track  data  set  looping  end  to  end. 

6.2. J  RIGID  BODY  CONCLUSIONS  AND  SUMMARY 

This  section  concludes  the  rigid  body  simulation  and 
analysis  phase  of  the  monorail  sled  verification  study. 
It  will  address  all  areas  posed  in  the  rail  roughness 
study  directive  as  it  pertains  to  monorail  sleds  in  a 
rigid  body  mode  of  operation. 

6.2. j.1  PERIODICITY  VERSUS  N ON P E R I  0 D I C I T Y 

As  with  dual  rail  sleds,  the  effect  of  periodicity 
may  be  examined  by  looking  at  the  statistical  results  of 
TSM000  top  and  the  south  1  5  00  0  bottom.  The  rails 
generate  the  same  standard  deviation  for  both  the 
vertical  and  lateral  responses.  The  periodic  rail, 
TS^tUOOO  top,  underpredicts  the  mean  by  60  percent  of  the 
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nonperiodic  rail  in  the  vertical  direction  and  by  19 
percent  in  the  lateral. 

Comparing  the  roughness  contribution  of  periodic 
rails  can  be  done  by  looking  at  TS  39000  versus  TS  99000, 
both  top  and  bottom.  Here  it  is  evident  that  rough 
periodicity  affects  both  the  mean,  by  as  much  as  200 
percent  o v e r p r e d i c t i o n  compared  to  the  smooth  rail,  and 
the  standard  deviation,  by  as  much  as  a  25  percent  larger 
prediction. 

The  extent  of  the  effect  of  periodicity  is  unknown 
since  a  proper  combination  of  rail  period  and  sled 
velocity  is  not  available  that  will  generate  the  250 
hertz  heave  frequency  of  the  sled.  Excitation  of  the  250 
hertz  frequency,  based  on  the  208- font  rail  period,  would 
require  a  sled  to  travel  at  an  unreasonable  velocity  of 
50000  feet  per  second.  On  the  other  hand,  if  a  periodic 
rail  length  of  39  feet  existed,  the  sled  could  be  excited 
at  its  fundamental  frequency  by  a  velocity  of  10000  feet 
per  second.  This  velocity  is  not  entirely  unreasonable 
since  the  Track  has  successfully  tested  sleds  at 
velocities  to  9000  feet  per  second. 


G.2.3.2  RAIL  GRINDING 


Comparison  of  means  and  standard  deviations  for  the 


north  and  south  15000  feet  of  track  can  be  done  t 


determine  rail  grinding  e f f e c t 1 v e n e s s  for  monorails. 
Locking  at  the  rail  roughness  functions,  it  is  readily 
seen  that  the  rougher  the  rail  the  larger  the 
o v e r p r e d  i  c t  i  o n  and  the  larger  the  standard  deviation 
spread  in  the  vertical  direction. 

It  was  hypothesized  that  the  long  period  frequency 
content  (rail  lengths  of  30  to  50  feet)  controlled  the 
standard  deviation.  If  these  rail  lengths  corresponded 
to  the  sled  frequencies  existing  above  the  fundamental 
frequencies,  i.e.,  the  secondary  frequencies,  a  lowering 
of  the  standard  deviation  could  be  expected.  For  a 
monorail  traveling  at  these  velocities,  these  rail 


lengths  correspond  to  forcing  function  frequencies  in  the 
range  of  6?  to  112  hertz,  well  below  the  fundamental 
vertical  sled  frequencies.  This  indicates  that  the 
secondary  frequency  effects  cannot  absorb  any  of  the  rail 
energy  and  thus  lower  the  standard  deviation. 

Looking  at  the  lateral  direction  it  is  seen  that  the 
means  are  relatively  close  to  each  other  and  do  not 
appear  to.  be  related  significantly  to  rail  roughness. 
However,  the  standard  deviation  for  the  north  15000  feet 
of  rail  is  less  for  both  rail  surfaces  than  is  the  south 
end  of  the  track.  This  should  be  expected  if  the 
hypothesis  is  true  about  secondary  effects.  The  forcing 
frequencies  for  the  long  period  rail  content  do  coincide 


with  frequencies  above  the  fundamental  resonances  in  this 


case. 

It  is  evident  that  rail  grinding  is  effective  in 
reducing  the  long  period  rail  content.  The  long  period 
rails  are  what  control  the  amplitude  swings  in  the 
dynamic  contributions.  This  fact  for  monorail  sleds  may 
not  necessarily  be  beneficial  unless  the  forcing  function 
frequencies  associated  with  these  long  period  rails 
coincide  with  the  sled's  secondary  structural 
frequencies.  This  frequency  overlap  is  typically  not  the 
case  in  the  vertical  direction,  but  can  be  so  in  the 
lateral  direction. 

6. 2. 3. 3  FINAL  RAIL  PROFILE  SUMMARY 

Examining  Table  6-3*  it  is  seen  that  all  rail 
functions  overpredict  actual  forces  by  quite  an  extreme 
range.  Noting  results  previously  cited  for  monorail 
work,  the  modal  participation  should  absorb  between  20 
and  30  percent  of  the  forces.  Therefore,  an  ideal  rail 
for  the  rigid  body  simulation  should  be  one  that 
overpredicts  the  actual  by  less  than  30  percent.  This 
condition  is  met  by  only  one  rail  function  in  the 
vertical  direction,  the  baseline  function.  However,  a 
second  rail  function  is  close  to  meeting  this  criterion 
and  should  be  considered.  It  is  the  top  of  the  rail  at 
TS  44000.  All  functions  meet  this  criterion  in  the 
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lateral  direction,  but  it  is  unknown  whether  the  modal 
participation  will  absorb  this  much  energy  in  the  lateral 
direction  since  lateral  prediction  has  never  before  been 
successful , 

Based  on  the  data  presented,  it  appears  that  rail 
roughnesses  in  excess  of  0.014  inches  tend  to  grossly 
overpredict  monorail  vertical  forces.  This  roughness 
value  corresponds  to  another  roughness  break  point 
classed  by  Mixon^.  It  is  the  class  existing  between  the 
second  and  third  defined  roughness  groupings. 

Based  on  the  rigid  body  work,  the  best  vertical 
simulation  is  made  by  rail  roughnesses  of  less  than  0.014 
inches.  Lateral  simulation  appears  to  be  relatively 
independent  of  roughness  although,  in  general,  rougher 
rails  tend  to  overpredict  more  than  smoother  rails.  Rail 
smoothing  on  monorails  does  not  appear  to  be  particularly 
effective  in  the  vertical  direction.  It  does  seem  to  be 
of  benefit  in  the  lateral  direction;  reductions  of  up  to 
2  5  percent  are  evident.  As  for  the  dual  rail  sleds,  rail 
grinding  affects  only  the  long  period  rail  content. 

0.3  MODAL  SIMULATION 

The  response  of  the  sled  system  in  the  rigid  body 
monorail  analysis  was  limited  to  the  direct  structural 
m~des  of  the  slipper  assemblies.  The  structural  modes 
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identifiable  with  the  sled  body  must  be  superimposed  onto 
the  rigid  body  system  to  accurately  reflect  the  true  sled 
behavior.  The  spike  monorail  sled  has  eight  secondary 
modes  identified,  two  of  which  are  second  lateral  bending 
modes  of  the  suspension  system.  These  frequencies  are 

summarized  in  Table  6-5  in  the  concluding  section  of  this 

\ 

chapter . 

Now  that  the  secondary  structural  response  has  been 
identified,  the  procedure  is  to  identically  repeat  the 
simulations  done  with  the  rigid  body  analysis.  The  modal 
analysis  amounts  to  running  19  different  trajectory 
conditions  for  each  of  the  nine  rail  functions,  or  a 
total  of  171  modal  simulations.  However,  it  was 
concluded  from  the  rigid  body  analysis  that  the 
conditions  existing  for  Mach  1  to  Mach  2  under  thrust 
were  invalid.  Therefore,  the  modal  simulation  consisted 
of  running  the  15  valid  trajectory  conditions  for  each  of 
the  nine  rail  functions,  or  1j5  total  computer  runs. 

Each  rail  roughness  simulation,  consisting  of  the  15 
valid  conditions,  was  statistically  averaged  to  determine 
the  effect  of  the  modal  participation.  The  results  of 
this  analysis  art  presented  as  an  average  deviation 
between  the  actual  forces  (see  Table  6-2  for  the  actual 
forces)  and  the  estimated  values  from  the  simulation 
sequence.  Also,  the  standard  deviation  about  the  mean 
deviation  has  been  calculated. 


The  effect  of  modal  participation  was  discussed  in 
Chapter  5  as  it  pertained  to  earlier  work  done  on  this 
monorail  sled.  It  was  concluded  that  six  modes  shapes 
and  five  percent  damping  would  reduce  the  percentage 
deviation  by  about  30  percent.  This  section  will  attempt 
to  verify  this  conclusion. 

6.3.1  DATA  ANALYSIS — AMPLITUDE  DOMAIN 

6. 3. 1.1  PEAK  FORCES 

The  modal  analysis  consisted  of  performing  135 
simulations  combining  the  monorail  rigid  body  system  with 
eight  structural  modes  at  three  percent  damping.  Results 
of  the  peak  force  predictions  for  this  analysis  are  shown 
plotted  in  Figures  6.15  through  6.13  for  the  four  slipper 
conditions  examined.  The  front  slipper  resultants  are 
shown  in  Figure  6.15  for  the  vertical  direction  and 
Figure  6.16  for  the  lateral  direction.  Figures  6.17  and 
6.18  contain  the  aft  slipper  vertical  and  lateral 
directions,  respectively.  The  solid  line  is  the  average 
force  described  in  Table  6-2,  and  the  dotted  lines 
correspond  to  the  maximum  and  minimum  values  for  that 
Mach  number. 

Examination  of  the  figures  indicate  the  vertical  and 
lateral  directions  for  the  monorail  sled  appear  to  be 
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o v e r p r e d i c t i n g  more  than  in  the  rigid  body  case.  The 
lateral  direction  overprediction  appears  to  be  at  least 
an  order  of  magnitude  larger  than  it  was  in  the  rigid 
body  case. 

The  statistical  deviation  analysis  is  presented  in 
Figure  6.19  and  summarized  in  Table  6-4  for  all  slipper 
conditions.  The  "x"  denotes  where  the  average  deviation 
between  actual  and  predicted  peak  forces  occurs.  The 
h e a v y- a r r o we d  line  represents  the  +  1  sigma  value  about 
the  mean  deviation. 

Looking  at  the  table,  it  is  seen  that  a  distinct 
break  point  in  o ve r pr ed i c t io n  in  the  vertical  direction 
occurs  at  a  rank  roughness  of  0.014  inches.  This 
distinct  break  point  applies  to  the  average  deviation  as 
well  as  the  standard  deviation  in  the  vertical  direction. 
A  defined  break  point  is  not  readily  apparent  in  the 
lateral  direction. 

Comparison  of  Tables  6-3  and  6-4  will  show  the 
effect  of  modal  participation.  It  was  earlier 
hypothesized  that  the  modal  participation  would  absorb  up 
to  30  percent  of  the  forces;  however,  it  is  seen  that  the 
average  peak  forces  increase  anywhere  from  16  to  60 
percent  in  the  vertical  direction.  The  standard  deviation 
increase  in  the  vertical  direction  remains  approximately 
the  same.  Lxamining  the  lateral  direction  indicates  a 
gross  increase  in  force  overprediction.  The  range  of 
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TABLE  6-4  STATISTICAL  SUMMARY 
of  MONORAIL  MODAL  PREDICTIONS 
VERSUS  RANK  ROUGHNESS 
FROM  MACH  2  THRUSTING  TO  MACH  2  COASTING 


Rail 

S 

y 

Vertical+++ 

Lateral+++ 

( in. ) 

Average 

1  Sigma 

Average 

1  Sigma 

34,000  Bottom++ 

0.048 

+323 .01 

83.0 

+398.81 

101.3 

34,000  Top++ 

0.043 

+126.61 

40.4 

+413.81 

135.5 

N1 5 , 000  Bottom 

0.020 

+264.71 

48.2 

+234.91 

98.2 

N1 5 , 000  Top 

0.016 

+213.51 

68.2 

+227.31 

78.7 

S15,000  Bottom 

0.015 

+167.31 

45.1 

+319.21 

134.7 

Baseline 

0.014 

+38.11 

13.8 

+345.21 

59.8 

44,000  Bottom++ 

0.013 

+92.11 

38.6 

+296.31 

75.4 

44,000  Top++ 

0.012 

+71.61 

33.9 

+287.11 

74.3 

S15,000  Top 

0.012 

+75.51 

28.1 

+316.21 

149.8 

Notes : 

++  Denotes  periodic  rail  content 

+++  Resultant  sum  of  both  slippers  acting  in  that  direction 

1.  All  values  listed  are  percentages 

2.  Positive  sign  denotes  overprediction 

3.  Negative  sign  denotes  underprediction 


cverprediction  increase  is  between  100  and  400  percent. 
The  standard  deviation  in  the  lateral  direction  increased 
also,  by  as  much  as  7  times  that  of  the  rigid  body  case. 

These  findings  are  in  direct  conflict  with  the 
monorail  modal  findings  expressed  by  MRI,  so  an  attempt 
was  made  to  find  out  why  this  discrepancy  occurred.  The 
MRI  work  was  reviewed  to  determine  exactly  what  procedure 
was  used  in  arriving  at  the  stated  results.  The 
following  excerpt,  found  on  page  36  of  Reference  7.  was 
found  and  explains  how  MRI  manipulated  the  data. 


..A  method  was  adopted  based  on  Mixon's  scheme  for 
reducing  data  from  the  Instrumented  Monorail  tests.  He 
extracted  the  peak  slipper  force  from  each  quarter  Mach 
number  velocity  range.  During  the  acceleration  phase, 
that  corresponded  to  about  .1  second  intervals.  This  was 
done  for  several  runs  giving  a  number  of  peak  amplitudes 
for  each  value  of  Mach  number.  These  were  averaged  and 
plotted  as  average  peak  versus  Mach  number.  The 
analogous  method  used  with  the  computer  simulation  is  to 
make  a  run  of  constant  velocity  and  look  at  the  peak 
amplitude  in  each  .1  second  interval.  A  run  of  1.0 
second,  for  example,  yields  ten  values  which  are  then 
averaged  . 


After  this  analysis  method  was  determined  to  have 
been  used,  a  routine  was  added  to  the  program  that  would 
average  the  peak  values  for  each  tenth  second  interval  of 
the  simulation  run.  A  simulation  was  performed  which 
matched  the  one  simulation  case  quoted  by  MRI.  These 
results  are  shown  in  Table  6-5. 


TABLE  6-5 


COMPARISON  OF  RESULTS 
SLEDYNE  VERSUS  DASTARR 


FRONT 

AFT 

Meas . 

Peak 

Average 

Meas 

Peak 

Average 

MRI-ver t ical 

16.00 

20.68 

16.61 

19.00 

23.28 

16.52 

DASTARR-vert 

22.76 

17.00 

17.89 

15.57 

DASTARR-lat 

7.25 

27.00 

13.63 

— 

4.90 

2  5.85 

13.26 

Notes : 

1.  Forces  are  in  kips 

2.  Condition  is  at  4000  fps  under  thrust 


As  cun  be  seen  from  the  table,  MRI's  results 
increased  when  modal  participation  was  added  and  the  peak 
averaging  was  not  done.  Results  in  the  vertical 
direction  compare  exceptionally  well  between  DASTARR  and 


SLEDYNE  for  both  peak  values  and  average  values. 


However,  the  peak  averaging  method  discussed  is  extremely 
questionable  and  should  be  avoided. 

6. 3. 1.2  PEAK  FORCE  VERSUS  TIME 

Contained  in  Figures  6.20a  and  6.21a  are  the  data 
clips  from  the  actual  monorail  sled  runs  discussed 
earlier.  Figures  6.20b  and  6.21b  are  the  simulated 
force-time  plots  for  the  front  vertical  and  lateral 
slippers,  respectively.  The  figures  represent  the  same 
velocity  condition  (3370  feet  per  second  under  thrust) 
for  the  south  15000  top  rail  function.  The  remaining 
eight  rail  function  histories  are  contained  in  Appendix 
C,  Figures  C.49  through  C.64. 

Actual  run  data  for  the  front  slipper  vertica. 
indicates  the  forces  stayed  in  an  envelope  between  +8 
kips  and  -12  kips,  and  centered  about  a  negative  one  or 
two  kips.  The  simulation  plots  for  the  various  rail 
profiles  indicate  the  same  centering  about  a  negative  one 
or  two  kips.  However,  the  envelope  size  varies,  and  in 
ail  cases  is  larger  than  the  measured  case. 
Additionally,  the  envelope  size  increases  with  increasing 
’-ail  roughness. 

The  front  slipper  lateral  forces  for  the  actual  run 
data  centered  about  zero  kips  and  ranged  between  plus  and 
minus  4  kips.  Examination  of  the  plots  for  the  front 
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lateral  slipper  exhibits  the  same  centering  about  zero 


kips,  and  modulates  an  equal  magnitude  on  either  side  o f 
zero.  The  range  of  modulation  is  up  to  six  times  as 
large  as  the  run  data,  and  appears  to  be  related  to  rail 
roughness . 

The  ringing  phenomenon  seen  in  the  lateral  direction 
for  the  rigid  body  case  is  again  present  for  the  modal 
simulations.  The  only  difference  is  that  after  the  first 
ringing  cycle,  the  amplitude  damps  to  the  level  of  the 
actual  lateral  measured  forces.  This  phenomenon  is 

referred  to  in  nonlinear  vibrations  as  s u p e r h a r m o n  i  c 

2  2 

response.  3 u p e r h a r mo n i c  response  occurs  when  the 
forcing  function  frequency  is  in  resonance  with  a  natural 
frequency,  and  is  also  an  integer  multiple  of  a  higher 
natural  frequency.  The  amplitude  of  the  higher  harmonic 
is  then  superimposed  on  the  resonant  frequency  amplitude. 
When  the  damping  is  relatively  small  the  superposition  of 
the  amplitudes  dies  out  within  one  or  two  cycles.  The 
rigid  body  simulation  did  not  contain  higher  frequencies 
that  could  be  superimposed  onto  the  fundamental 
frequency.  The  first  elastic  lateral  mode  occurred  at 
four  (  M  )  times  the  forcing  frequency  and  thus  became 
superimposed  upon  the  rigid  body  frequency  for  one  cycle. 

Examination  of  the  simulated  force-time  history  for 
the  south  15000  top  indicates  that  the  peak  forces  occur 


where  the  localized  rough  spots  exist  for  this  rjil 


function,  and  repeat  as  the  profile  loops  end  to  end. 
This  phenomenon  is  that  same  one  that  was  discussed  in 
the  rigid  body  dual  rail  work  simulation  for  this  rail 
function  . 

6.  j  .  1  . j  FORCE  HISTOGRAM 

The  discussion  on  force  histograms  for  rigid  body 
monorails  indicated  that  the  force  distribution  was 
strongly  dependent  op  the  mean,  which  in  reality  is  the 
quasi-steady-state  portion  of  the  forces.  The  smaller 
the  mean,  the  closer  the  approximation  is  to  a  normal 
distribution.  The  larger  the  quasi-steady-state 
contribution,  the  more  like  a  Rayleigh  distribution  the 
estimated  forces  become.  The  quasi-steady-state  content 
for  the  modal  monorail  is  small  so  the  expected 
distribution  should  be  a  slightly  skewed  Gaussian  one. 
Shown  in  Figure  6.22  is  a  typical  modal  monorail 
distribution,  and  it  is  readily  seen  that  it  closely 
approximates  a  normal  distribution  due  to  the  small  mean 
value  . 

6  .  j . 2  DATA  ANALYSIS--FREQUENCY  DOMAIN 

The  highest  frequency  content  for  the  rigid  body 
monorail  sled  simulation  was  about  2600  hertz.  Repeating 


the  simulation  sequence  for  the  modal  participation  adds 
a  sufficient  number  of  data  points  to  increase  the  upper 
frequency  limit  to  approximately  13,000  hertz.  This 
frequency  is  well  beyond  the  normal  structural  excitation 
range,  so  the  analysis  will  only  examine  the  response 
spectrum  through  2500  hertz. 

Shown  in  Figures  6.23a  and  6.24a  are  the  typical 
wide  spectrum  response  for  the  modal  monorail  simulations 
for  the  front  slipper  in  the  vertical  and  lateral 
directions,  respectively.  These  figures  are  directly 
comparable  to  those  shown  in  Figures  6.13a  and  6.14a  for 
the  actual  run  data.  Figures  6.23b  and  6.24b  are  blowups 
of  Figures  6.23a  and  6.24a,  respectively,  in  the  0  to 
1000  hertz  range  and  clipping  the  power  axes.  Using  the 
later  figures,  it  will  be  possible  to  isolate  the 
individual  frequency  content  of  the  system. 

The  response  of  these  slippers  indicate  typical 
monorail  frequency  content.  Maximum  energy  is  contained 
in  an  energy  band  between  100  and  200  hertz  with  sharp 
rolloff  above  this  value.  Noise  predominates  the  energy 
above  600  hertz. 

Comparison  of  the  actual  versus  predicted  PSDs 
indicates  frequency  contents  that  are  almost  identical. 
The  frequency  content  from  the  actual  run,  the  rigid  body 
simulation,  and  the  modal  simulation  are  summarized  in 
the  concluding  section  in  Table  6-6. 
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Power  levels  for  the  frequencies  match  up  quite  well 
with  the  values  seen  in  the  rigid  body  case,  yet  these 
levels  are  still  in  excess  of  the  actual  run  case. 
Examination  of  the  plots  again  show  the  distinct  break 
point  seen  at  a  rail  roughness  of  0.014  inches  in  the 
amount  of  excess  energy  contained.  Roughnesses  less  than 
this  break  point  have  energy  content  only  marginally  in 
excess  of  the  actual  run  data. 

The  frequency  content  associated  with  the  track 
survey  iteration  distance  (65  hertz)  and  the  39-foot  rail 
lengths  (86  hertz)  are  masked  by  the  data  filtering  of 
the  PSD  program.  The  program,  for  this  upper  frequency 
limit,  filters  the  data  in  25-hertz  increments; 
therefore,  these  particular  frequencies  are  combined  with 
actual  sled  frequencies.  The  ringing  seen  in  the  lateral 
direction  occurs  at  a  frequency  of  about  55  hertz,  which 
corresponds  to  the  fundamental  lateral  frequency  and 
appears  to  be  excited  by  the  50  foot  survey  iterations. 
Again  the  ringing  envelopes  repeat  at  the  corresponding 
frequency  for  the  data  looping  end  to  end. 

6.3*3  MODAL  CONCLUSIONS  AND  SUMMARY 

This  section  concludes  the  modal  simulation  analysis 
for  monorail  sleds.  It  will  address  all  questions  posed 
in  the  rail  roughness  study  directive  from  a  modal 
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participjtior  point  of  view. 

C.j.j.1  PERIODICITY  VERSUS  N 0 N P E R 1 0 D I C I T Y 

As  with  the  rigid  body  monorail  work,  the  smooth 
periodic  rail's  prediction  is  closer  to  the  actual 
measured  forces  than  is  the  nonperiodic  rail.  Both  rails 
experienced  an  increase  in  the  average  deviation  from  the 
rigid  body  to  the  modal  participation  case.  The  periodic 
rail  deviation  increased  by  3  5  percent  in  the  vertical 
direction,  whereas  the  nonperiodic  rail  increased  by  only 
13  percent.  The  standard  deviation  for  both  rails 
increased  by  approximately  10  percent. 

Comparing  the  roughness  contribution  of  periodic 
rails  can  be  done  by  looking  at  T  S  3  4  0  0  0  versus  T  S  4  4  0  0  0  . 
The  rough  periodic  rail  content  average  o v e r p r e d i c t  i  o n 
increased  by  about  50  percent  from  the  rigid  body  to  the 
modal  case.  The  rough  rails  still  grossly  overpredict 
compared  to  the  smooth  rails  by  as  much  as  a  factor  of 
two  on  the  average  deviation  for  the  vertical  direction. 
The  standard  deviations  for  the  rough  periodic  rail 
remained  relatively  unchanged  between  simulation  cases, 
yet  is  also  up  to  twice  as  large  as  the  smooth  periodic 
rail. 

It  appears  that  the  increase  in  force  prediction  in 
the  vertical  direction  is  due  to  the  sled  bending  modes 
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being  excited.  The  tie  downs  can  excite  the  first 
vertical  bending  mode  at  a  velocity  of  13  3  0  feet  per 
second.  Similarly,  the  second  bending  mode  in  this 
direction  corresponds  to  the  frequency  generated  by  the 
tie  downs  at  Mach  3.0.  Examining  data  output  showed  both 
of  these  sled  modes  exhibiting  significant  acceleration 
levels  in  the  vertical  direction.  Laterally,  the  tie 
downs  can  also  excite  the  first  and  second  bending  modes 
with  sled  velocities  at  Mach  1.0  and  3.5,  respectively*. 
The  first  lateral  bending  occurs  at  a  frequency  low 
enough  where  significant  structural  excitement  is 
possible  and  explains  the  nearly  doubling  of  the  lateral 
force  predictions. 

6. 3. 3. 2  RAIL  GRINDING 

It  was  concluded  in  the  rigid  body  section  for 
monorails  that  rail  grinding  was  virtually  ineffective  in 
the  vertical  direction.  The  results  of  the  modal 
analysis  arrive  at  the  same  conclusion. 

Laterally,  the  rigid  body  analysis  concluded  rail 
grinding  is  effective  in  reducing  the  dynamic  operating 
envelope  seen  by  the  sleds.  Dased  on  the  modal  results, 
the  same  conclusion  is  drawn.  The  dynamic  envelope  is 
smaller  for  the  ground  ra i 1  than  the  unground  Mil; 
however,  the  envelope  size  is  still  5  times  larger  than 


in  the  rigid  body  case 


The  increase  in  envelope  size  is 


a  direct  result  of  the  tie  down  spacing  exciting  the  sled 
Laterally. 

C.j.j.j  FINAL  RAIL  PROFILE  SUMMARY 

The  rigid  body  work  concluded  that  the  baseline  rail 
function  would  be  the  ideal  rai  1  to  use  in  the  vertical 
direction  since  it  only  overpredicted  by  2 5  percent. 


Consideration  of  modal  energy  absorbption  hypothesized 


this  overprediction  would  be  reduced  to  zero.  After 
performing  the  modal  analysis,  it  was  found  that  the  tie 
down  spacing  excites  the  sled  bending  modes  enough  that 
force  reduction  does  not  occur.  However,  the  modal 
results  indicate  the  baseline  is  still  the  best  rail 
profile  since  the  overprediction  only  increased  slightly, 
l  .  e  .  ,  12  percent,  and  the  standard  deviation  actually 
decreased . 

G.j.j.4  MODAL  SUMMARY 

Utilizing  all  data  examined  thus  far,  it  is  inherent 
In  at  rail  roughnesses  in  excess  of  0.014  inches 
significantly  increase  the  forces  subjected  to  the  sled 
when  modal  participation  is  added  if  the  model  is 
uccuratc.  This  breaK  point  holds  for  both  front  and  aft 


slippers  for  both  the  average  deviation  and  the  standard 
deviation . 

The  modal  content  for  monorails  affects  the  forces 
in  varying  degrees  but  in  all  cases  the  average  deviation 
increases  1  r  r  r,  the  rigid  body  case,  particularly  in  the 
lateral  direction.  Vertically,  the  increase  ranges  from 
' b  to  SO  percent  and  is  strongly  dependent  op  rail 
roughness  and  long  period  undulations. 

Based  o p  the  modal  simulations,  the  p  e  n  k  slipper 
forces  can  be  reduced  significantly,  even  ri'-re  than  seen 
for  dual  rail  sleds,  by  reducing  rank  roughness.  The 
reduction  for  monorails  is  not  as  straight  forward  as  the 
dual  rails  but  u  two  to  three  percent  reduction  in  the 
vertical  overprediction  is  possible  for  each  0.001  inch 
reduction  in  rank  rail  roughness. 

Rail  smoothing  appears  to  be  effective  in  the 
lateral  direction,  but  rail  alignment  at  the  tie  downs 
seems  to  be  more  critical. 

6.4  MONORAIL  SUMMARY  AND  CONCLUSIONS 

An  extreme  number  of  computer  runs  were  performed  in 
an  effort  to  verify  DA STARR  for  monorail  sleds.  The 
simulations  were  all  performed  using  the  constant 
velocity  option  for  a  single  sled,  and  it  was  found  that 
the  monorail  sled  was  particularly  sensitive  t  o  r-  i  1 
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nughness  data  transients  such  as  the  tie  dnwn  spacing 


and  survey  iteration  distances. 

6.4.1  MONORAIL  VARIABLE  VELOCITY  SIMULATION 

Making  note  of  the  monorail  sensitivites  tc  the  rail 
function  data  irregularities  and  the  uncorrelated  lift 
functions  through  Mach  2.0,  an  attempt  was  made  to 
perform  the  variable  velocity  simulations  knowing  the 
results  might  be  suspect.  The  results  of  this  attempt 
were  very  discouraging.  The  rigid  body  modes,  as  the 
sled  traveled  through  the  velocity  regime,  became  more 
and  more  excited  yielding  a  diverging  solution  that 
resulted  in  totally  unrealistic  slipper  forces.  Final 
force  predictions  for  the  rigid  body  analysis  approached 
b  x  10^  pounds,  which  indicated  the  calculated  slipper 
deflections  were  allowed  to  reach  impossible  proportions, 
i . e . ,  13b  inches. 

Results  for  the  modal  variable  velocity  simulation 
reached  the  same  conclusion,  only  the  peak  forces 
maximized  at  one  million  pounds.  The  times  of  peak 
occurrence  were  identical  in  both  cases  and  equated  to  a 
Mach  number  of  three.  At  this  time  no  viable  explanation 
for  the  divergence  of  the  solution  is  possible.  Further 
investigation  into  this  anomaly  is  recommended. 
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6.4.2  CONCLUSIONS 


It  is  apparent  that  DASTARR  can  predict  peak  slipper 
forces  for  monorail  sleds  in  a  manner  that  is  totally 
acceptable  to  engineering  judgment.  The  results  .'or  the 
monorail  are  strongly  dependent  on  the  rail  function, 
velocity  option,  and  modal  considerations  used. 

It  is  concluded  that  the  constant  velocity  function 
for  a  single  rigid  body  monorail  yields  the  best  results. 
The  results  become  even  more  attractive  when  the  data 
averaging  technique  suggested  by  M R  I  is  applied. 
However,  it  is  strongly  recommended  that  this  option  not 
be  used  . 

It  was  shown  in  this  chapter  that  peak  forces  for 
monorail  sleds  fall  between  a  Rayleigh  and  a  normal 
distribution.  The  amount  of  contribution  of  each 
function  is  strongly  affected  by  the  quasi-steady-state 
forces,  and  the  dynamic  envelope.  The  smaller  the  quasi¬ 
steady-state  force  contribution,  the  closer  the  data 
approximates  a  Gaussian  distribution.  The  skewness  of 
the  distribution  is  then  controlled  by  the  standard 
deviation  . 

Rail  grinding  appears  to  be  marginally  effective  in 
the  lateral  direction.  The  simulation  results  indicate 
the  tie  down  spacing  alignment  is  more  critical  in 
reducing  forces  for  monorails  than  is  rail  grinding.  In 


general,  a  two  to  three  percent  reduction  in  forces  is 
"btdinjblo  f  o r  each  0.001  inch  reduction  in  rank  rail 
roughness.  Hail  smoothing  is  particularly  effective 
until  <j  rank  roughness  of  0.014  inches  for  monorails  has 
been  achieved.  Results  do  not  improve  significantly  by 
smoothing  to  less  than  this  value. 

Verification  of  localized  rough  spots  on  the  Track 
have  been  accomplished  for  monorails.  The  effect  of 
localized  polishing  can  now  be  studied  to  determine 
optimal  grinding  values. 

When  considering  rail  functions,  the  baseline 
function  generates  the  best  slipper  estimates  for  both 
the  rigid  body  and  rnodal  cases.  Peak  forces,  in  the 
lateral  direction,  have  been  correlated  for  absolute 
magnitude  values  for  the  first  time,  though  actual 
simulation  force-time  plots  do  not  mirror  real  run 
results.  The  simulation  plots  exhibit  a  ringing 
sensation  that  resembles  narrow  band  response.  Narrow 
band  response  is  when  a  single  degree  of  freedom 
oscillator,  or  a  sharply  resonant  system,  is  subjected  to 
a  forcing  function  whose  energy  is  relatively  constant 

over  a  substantial  frequency  range  when  compared  to  the 

2  2 

oscillator  frequency.  The  resultant  is  defined  as  a 
narrow  band  response  spectrum.  Laterally,  the  math  model 
is  idealized  as  a  single  degree  of  freedom  system  for  the 
rigid  body  case,  or  as  a  sharply  resonant  system  for  the 


modal  analysis.  The  forcing  function  is  relatively  flat 
over  a  frequency  range  well  in  excess  ^  f  the  system 
resonant  frequency.  Much  of  this  energy  is  at 
frequencies  caused  by  the  data  anomalies  generated  by 
survey  techniques  and  data  set  looping.  These 
frequencies  do  not  occur  in  reality  and  account  for  the 
lateral  direction  simulation  discrepancies  encountered. 

Frequency  contents  found  in  the  simulation  are 
summarized  in  Table  6-6.  Excellent  correlation  was 
obtained  in  the  vertical  direction.  The  lateral 
direction  exhibited  i n c o ns i s t an c i e s  that  have  just  been 
d i scussed  . 

Comparing  the  modal  analysis  to  the  rigid  body  case 
indicates  that  for  monorail  sleds  a  rigid  body  analysis 
will  adequately  define  the  slipper  loads.  Modal 
participation  does  not  refine  the  structural  loads  as 
hypothesized,  in  fact  the  loads  increase  due  t  ° 
structural  modes'  excitement  caused  by  the  rail  tie 
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CHAPTER  7 


CONCLUSIONS  AND  RECOMMENDATIONS 

7.1  GENERAL 

Detailed  in  Chapter  One  were  the  proposed  areas  of 
study  to  be  accomplished  during  this  research.  These 
were  expanded  under  the  rail  roughness  study  directive 
described  in  Chapter  Four.  These  objectives  are 
summ  arized  again. 

1)  Define  the  equations  of  motion  for  up  to  three 
coupled  sleds  for  six  degrees  of  freedom.  These 
equations  were  to  be  defined  such  that  they  would 
collapse  to  predict  one  sled  as  well  as  three. 
Finally,  the  equations  were  to  collapse  such  that 
pitch  plane  simulation  could  be  achieved  so  the 
results  would  mirror  the  existing  SLEDYNE  program. 

2)  Allowance  for  a  variable  velocity  input  instead 
of  being  limited  to  a  constant  velocity  simulation. 

3)  Allowance  for  variable  slipper  gap  definitions. 
U)  Allowance  for  an  easier  slipper  stiffness 
definition  that  will  eliminate  the  need  for 
modifying  the  N  AST  R  AN  math  model.  Presently  a 
secondary  modal  analysis  that  is  compatible  with 


SLEDYNE  requirements  is  necessary. 

5  )  Establish  a  final  and  fixed  rail  forcing 

function  to  use  for  all  simulations. 

6)  Establish  a  fixed  length  of  rail  for  optimizing 
force  predictions  for  the  constant  velocity 
simulation  cases. 

7)  Determine  the  effect  of  long  period  rail 
undulations  on  slipper  impact  forces. 

8)  Establish  if  rail  grinding  is  effective  in 
reducing  the  slipper  impact  forces. 

7.2  CONCLUSIONS 

Each  of  the  proposed  areas  will  be  discussed  and 
general  conclusions  drawn  about  them.  The  next  section, 
Chapter  7.3,  will  discuss  the  limitations  found  and  make 
recommendations  for  future  research  to  expand  on  the  work 
that  has  been  accomplished  here. 

7.2.1  DEFINITION  OF  THE  EQUATIONS  OF  MOTION 

The  equations  of  motion  were  derived  for  coupled 
sled  systems  in  Chapter  Two,  The  equations  were 
developed  in  matrix  form  to  ease  the  required 
computational  effort.  Since  the  equations  were  derived 
in  matrix  form,  the  matrices  were  allowed  to  collapse  to 
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handle  either  one,  two,  or  three  sleds  coupled  together. 
These  equations  were  defined  to  allow  the  system  of 
matrices  to  further  collapse  into  only  a  pitch  plane 
analysis  . 

The  validity  of  this  collapsing  capability  was 
verified  in  Chapter  Four.  A  sled  that  had  been 
extensively  analyzed  in  the  pitch  plane  using  the  bLEDYNt 
program  was  subjected  to  an  analysis  on  DASTARR  in  the 
one  sled,  pitch  plane  only  configuration.  Adequate 
correlation  of  this  configuration  was  obtained. 

The  equations  of  motion  for  a  single  sled,  either 
dual  rail  or  monorail,  at  a  constant  velocity  were 
correlated  in  Chapters  Five  and  Six.  This  correlation 
examined  the  difference  between  performing  a  ri  '  ’  body 
simulation  and  the  superposition  of  modal  properties  onto 
that  rigid  body  system. 

Not  discussed  previously  is  the  effect  of  sled 
coupling,  due  to  an  insufficient  number  of  appropriate 
sleds  having  adequate  math  models  to  allow  for  the  proper 
sled  train  configurations.  However,  this  option  of  the 
program  has  been  tested  through  a  series  of  fictitious 
rigid  body  sled  trains.  These  fictitious  rigid  b^d  y 
sleds  had  properties  defined  that  were  op  the  order  of 
magnitude  of  real  sled  properties  used  under  various 
track  coupling  configurations.  Results  of  these 
simulations  gave  values  that  appear  to  be  valid,  or  at 


least  fall  into  an  expected  range  of  forces  based  on 
trdck  expedience  for  similar  systems.  It  should  be 
cHveated  that  transmitted  thrust  transient  predictions 
cannot  be  correlated  until  these  math  models  nave  been 
developed  and  verified. 

1.2.2  VARIABLE  VELOCITY  OPTION 

The  program  was  written  to  allow  the  user  to  specify 
a  variable  velocity  function  for  use  in  the  simulation. 
This  velocity  function  is  defined  to  be  an  output  file 
from  the  Test  Track  trajectory  program,  TSLED.  Using 
specific  parameters  from  the  velocity  profile,  the 
program  can  integrate  the  equations  of  motion  to  obtain 
the  required  force  output.  The  validity  of  the  results 
using  the  variable  velocity  were  confirmed  in  Chapter 
five  for  dual  rail  sleds.  Due  to  rail  function 
idiosyncracies  modulated  by  velocity,  the  monorail 
solutions  tended  to  diverge  under  this  option  and  have 
not  yet  been  correlated. 

7.2.3  VARIABLE  SLIPPER  GAPS 

The  program  was  designed  to  allow  the  user  to  define 
a  slipper  gap.  This  gap  was  then  assumed  to  be  constant 
throughout  the  computer  simulation.  This  gap  assumption 
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does  not  meet  the  requirement  levied  for  a  variable 
slipper  gap;  however,  one  may  vary  a  constant  gap  between 
simulations  to  determine  the  effect  of  the  gap.  No 
attempt  has  been  made  to  determine  the  "ideal"  gap  or  to 
study  the  effects  of  varying  the  gap  thickness. 

The  need  for  a  variable  slipper  gap  stems  from  the 
variations  in  the  cap  thickness  for  the  railhead  and  to 
slipper  wear  occurring  during  the  sled  run.  This 
parameter  of  the  program  can  eventually  be  modified  to 
account  for  these  effects.  An  extensive  amount  of  survey 
data,  which  is  unavailable  to  date,  is  required  to 
account  for  the  cap  deviations.  Once  this  data  is 
obtained,  the  program  can  be  modified  in  the  way  it 
handles  the  rail  roughness  parameter  to  yield  a  variable 
slipper  gap.  The  slipper  wear  problem  cannot  be 
addressed  until  sufficient  data  is  available  to  define 
the  slipper  wear  as  a  function  of  velocity.  When  this 
happens ,  an  additional  table  may  be  incorporated  into  the 
program  to  handle  gap  increases  due  to  the  slipper  wear. 

7.2.4  SLIPPER  STIFFNESS  DEFINITION 

Under  the  existing  SLEDYNE  program,  the  modal 
participation  was  limited  to  the  secondary  structural 
modes.  Definition  of  these  modes  was  accomplished  by 
redefining  the  structural  constraints  of  the  model  to 


eliminate  the  "rigid  body"  contributions.  This  model 
redefinition  resulted  in  defining  two  math  models  to 
examine  the  sled,  and  still  did  not  address 
simplification  of  the  slipper  stiffness  definition. 

Based  on  DASTARH  programming,  the  user  is  allowed  to 
define  the  modes  to  be  used  in  the  simulation.  This 
modal  selection  lets  the  user  develop  one  math  model  with 
realistic  support  constraints  and  then  delete  the  defined 
rigid  body  modes  in  the  modal  definitions.  This 
procedure  eliminates  the  need  to  define  two  math  models, 
but  again  the  slipper  stiffness  definition  is  not 
simplified. 

One  can  take  the  approach  suggested  by  Mixon*'  for 
slipper  stiffness;  that  is,  develop  an  extensive  math 
model  of  the  structural  suspension  system,  statically 
load  it,  and  calculate  the  stiffness  based  on  the 
deflection.  This  approach  is  the  best  method  for 
monorail  sleds  since  the  slipper  itself  is  the  primary 
suspension  clement.  For  dual  rail  sleds,  one  may 
determine  the  slipper  beam  stiffness  based  on  the  unit 
dummy  load  method  to  calculate  a  deflection  under  a  given 
load.  This  method  is  relatively  simple  to  apply  for  dual 
rail  sleds  since  the  beup  spanning  the  rails  is  the 
primary  suspension  element.  The  approach  used  for 
monorails  is  best  for  special  case  dual  rail  sleds,  those 
without  standard  slipper  beams. 


An  alternative  method  to  the  above  procedures  would 
Di  to  include  all  suspension  elements  in  the  modal 
analysis.  "sing  the  calculated  frequency  for  the  rigid 
bojy  motions,  ope  can  derive  the  required  stiffnesses 
based  on  the  proper  apportionment  of  the  sled  mass.  This 
procedure  was  alluded  to  in  Table  5-6  for  dual  rail 
sleds,  under  the  rigid  body  modes  heading  that 
apportioned  the  mass  to  calculate  the  rigid  body  heave 
frequenc ies  . 


7.2.5  FINAL  AND  FIXED  RAIL  FUNCTION 


Based  on  the  discussions  of  Chapters  Five  and  Six, 
it  was  concluded  that  the  Baseline  function  yielded  the 
best  results  under  all  conditions,  with  the  function  for 
TS  it^OOO  top  being  a  secondary  choice.  While  these 
profiles  yielded  acceptable  results,  they  introduced 
facsimilies  into  the  data  that  may  or  may  not  exist. 
These  data  irregularities  have  been  previously  discussed 
and  are  a  result  of  survey  techniques  used  to  o bt jin  the 
preliminary  rail  roughness  measurements  and  the  looping 
of  the  data  set  to  generate  .1  i  1  of  sufficient  length. 


7.2.6  FIXED  RAIL  LENGTH 


B  a S"d  on  all 


ail  functions  studied. 


i  f 1  x  e  J  long 1 1 , 


o  f  4000  feet  appears  to  optimize  rail  force  predictions. 
R . ,  1  i  lengths  short  e r  than  4000  feet  have  not  reached  the 
peak  force  predictions,  and  lengths  in  excess  of  this 
v  a  1 ue  Jo  not  increase  the  for ce  prediction  estimates. 


LUNG  PLR10D  UNDULATION 


Large  amplitude,  long  period  rail  undulations  do 
increase  the  predicted  impact  forces  provided  the  proper 
p  e  r  i  o  j  length,  sled  velocity,  and  fundamental  frequency 
a r e  consistent.  That  is,  a  typical  dual  rail  sled  having 
a  fundamental  heave  frequency  of  i  ij  hertz  and  traveling 
at  a  velocity  of  jooo  feet  per  second  can  have  this  mode 
excited  by  a  r  o  u  g  hi  rail  period  of  approximately  200  feet. 
Similarly,  secondary  structural  modes  can  be  excited  by 
trie  proper  c~rn  binut  ion  of  sled  frequency,  sled  velocity, 
and  physical  rail  period  lengths.  These  periodic  lengths 
a  r  (.•  not  limited  to  actual  rail  undulations  ,  but  can  be 
any  physical  rail  parameters,  such  as  rail  length  (  j  9 


t  *  •  i.  t  )  rr  tie  irWp  spacing  (02  inches). 


HAIL  GRINDING  EFFECTIVENESS 


•  i  i  1  g  '  l  v  d  i  n  g  1.  a  s  teen  determined  to  be  effectiv  e 
nut  not  as  f ■  r  i  g  i  r  „ 1  1  y  envisioned.  It  was  thought  that 


.ii  r  i  n  j  i  n  ,  „  a  ;  :  m  o  -  t.  h  t  ft  e  r  a  i  1  between  tie  downs. 


to  existing  measured  sled  slipper  data  for  both  dual  rail 
and  monorail  sled  systems,  ir  both  the  vertical  and 
lateral  directions.  The  Test  Track  now  has  a  viable 
design  tool  in  DASTARK  for  predicting  sled  track 
interaction  forces.  Additionally,  the  program  was 
developed  to  allow  sleds  to  be  coupled  so  that 
transmitted  thrust  transients  could  be  examined  for  the 
first  time.  With  this  new  design  tool,  sleds  can  now  be 
designed  more  efficiently  to  reflect  an  optimum  strength 
to  weight  ratio.  This  optimum  strength  to  weight  ratio 
will  allow  for  using  existing  propulsion  to  accomplish 
attaining  tne  higher  velocity  requirements  safely. 

7  .  J  R  E  C  u  M  M  E  N  D  A  T  1  0  N  S 

The  rcsejrcl:  involved  in  developing  the  DASTARR 
simulation  program  resulted  in  a  viable  design  tool  for 
tne  Test  T  r  a  e  k  .  The  program  has  been  verified  as  to  the 
quality  of  the  results;  however,  there  are  still  areas 
that  need  to  be  investigated.  These  areas  are  discussed 
below  and  are  not  given  in  any  priority. 

7  .  j  .  1  TRAN3MITTED  THRUST  TRANSIENTS 

It  was  discussed  in  tne  previous  section  that 
transmitted  thrust  transients  we  *■  <•  not  ,  '  r  r  ^  I  j  i  c  j  ,  a  s  a 


direct  result  of  not  having  appropriate  math  models 
available.  Verification  of  the  transmitted  thrust 
feature  of  DASTARR  must  wait  until  these  models  have  been 
developed.  Additionally,  sufficient  track  run  data  o  n 
these  transients  must  be  taken  to  obtain  an  adequate  data 
base  for  comparison  with  the  predicted  forces. 

7.3.2  RAIL  ROUGHNESS  PROFILE 

It  was  concluded  that  the  Baseline  rail  profile 
yielded  the  best  force  prediction  results.  Even  though 
the  results  were  correlated,  the  survey  techniques  used 
to  obtain  the  data  introduced  errors  into  the  functions. 
A  forcing  function  should  exist  that  spans  at  least  4900 
feet  of  track  in  order  to  accurately  estimate  rail  forces 
in  the  future.  A  rail  of  this  length  will  eliminate  the 
cycling  effect  of  looping  the  data  end  to  end.  Also,  a 
new  survey  procedure  must  be  established  that  will 
eliminate  the  50-foot  sections  of  data  taking.  Something 
along  the  line  of  a  laser  survey  that  can  measure  the 
roughness  for  at  least  5000  feet  of  continuous  track  is 
desirable.  The  entire  50,000  plus  feet  of  track  should 
be  measured,  for  all  surfaces  on  each  of  the  rails  to 
completely  delineate  the  rail  profile.  When  the  complete 
rail  has  been  defined,  the  study  should  be  repeated  with 
the  larger  data  base  to  determine  whether  a  given  section 


of  rail  will  satisfactorily  predict  slipper  forces,  or 
whether  the  entire  rail  definition  is  required. 

Y.j.j  MASS  VARIATIONS 

Under  the  present  version  of  DASTARR,  mass 
variations  of  the  sled  cannot  be  studied  except  in  the 
constant  velocity  case.  Tru  t  1  m  e  -  v  u  r  y  i  r  g  mass  problems 
associated  with  pusher  sleds  as  the  propulsion  is 
expended  is  precluded.  Preliminary  studies  have 
indicated  that  varying  the  mass,  and  its  associated 
center  of  gravity,  have  a  pronounced  effect  on  the 
results.  The  program  can  be  modified  to  allow  a  time  or 
velocity  dependent  mass  and  center  of  gravity  to  be 
defined.  This  modification  will  take  care  of  the  rigid 
body  contribution  to  the  total  slipper  forces  due  to  the 
mass  changing.  However,  the  extent  of  the  modal 
contributions  cannot  be  determined  until  an  extensive 
parametric  study  has  been  performed  investigating  how  the 
mode  shapes  change  as  a  function  of  the  mass  change.  The 
mass  variation  is  an  important  area  of  interest  as 
certain  rocket  motors,  particularly  at  burnout,  become 
longitudinally  excited  and  this  excitement  in  turn  can 
excite  sled  Jowntrack  modes  and  then  induce  transient 
spues  into  the  transmitted  thrust  vector. 


7. 3-1*  LATERAL  FORCING  FUNCTIONS 


The  program,  as  written,  does  not  account  for 
possible  secondary  forcing  functions,  such  as  those  in 
the  lateral  direction.  The  lack  of  this  definition  can 
be  remedied  by  allowing  additional  forcing  function  table 
definitions,  such  as  was  done  for  the  lift  and  center  of 
lift  tables.  The  addition  of  this  feature  will  be 
appropriate  to  study  the  design  requirements  op  sleds 
subject  to  explosion  testing.  Explosion  testing  is  when 
explosives  are  detonated  beside  the  track  to  study  blast- 
created  pressure  waves  and  how  these  waves  affect  various 
test  platforms  as  the  vehicle  travels  through  the 
pressure  wave.  Another  potential  use  for  this  feature 
will  be  to  define  cross  track  aerodynamic  forces  when  the 
sled,  particularly  monorails,  yaws  with  respect  to  the 
rail  centerline. 

7.3.5  RAIL  GOUGING 

A  problem  encountered  on  high  speed  monorail  sleds 
is  a  phenomenon  known  as  rail  gouging.  Rail  gouging  is 
the  physical  removal  of  a  piece  of  track  when  the  sled 
travels  over  a  section  of  rail.  There  are  several 
theories  to  explain  gouging,  and  the  two  most  likely  are 
concerned  with  the  physical  interaction  of  the  slipper 
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with  the  rail. 


One  theory  holds  that  matter  particulate,  probably 
debris  from  slipper  wear,  is  caught  between  the  slipper 
and  the  rail.  Under  this  condition  of  infinitely  high 
pressure,  the  particulate  flashes  into  a  plasma  and 
literally  blows  a  piece  of  rail  out.  The  second  theory 
states  that  at  these  high  velocities,  the  slipper 
impacting  the  rail  creates  an  extremely  high  localized 
pressure  that  causes  an  instantaneous  weld  to  occur 
between  the  rail  and  slipper.  Then,  as  the  sled 
continues  down  the  track,  the  slipper  pulls  a  section  of 
r  a  1  L  oUt,. 

The  second  of  these  theories  has  the  potentiality  to 
be  investigated  with  DASTARR.  The  program  can  be 
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7 . J . 6  RAIL  GRINDING 


Although 

r  a i 1  grinding  is 

effective 

i  n 

reducing 

the 

dynamic  force 

envelope,  proper 

al  ignment 

O  f 

the  rail 

a  t 

L h e  tie  downs 

will  result  in 

a  greater 

for 

ce  reducti 

on  . 

It  is  recommended  that  rail  grinding  be  suspended  in  lieu 
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APPENDIX  A 


DERIVATION  OF  DAMPING  COEFFICIENTS 

To  compute  the  damping  coefficients  associated  with 
the  flexing  of  the  slipper  beams,  consider  the  following 
figure  for  bounce  conditions. 


FIGURE  A .  1  VERTICAL  BOUNCE  DAMPING  COEFFICIENT 

By  summing  forces  in  the  vertical  direction,  the  equation 
of  motion  can  be  derived  as 

m’i  +  Krz*K  z  +  C  z  =  G  (AT) 

I  ri  L 

who r  e  K  (Jnd  K .  a  o  e  the  sun  of  the  stiffnesses  for  the 

I  ti 

fore  and  aft  slipper  assemblies,  respectively.  Equation 


+  L 
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1  F  r  +  1  F  =  0 
f  1  a  a 


Substituting  equation  A  8  ini''  A  (J  and  dividing  u  y  I 


( ! . “ K _  +  l  "K  i 
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(  A  !  0  ) 


So  the  natural  frequency  for  pitch,.  ,  is  given  by 


K  r  +  1  “  K 

if  a  a 


(AIM 


Also  f  r  o  m  equation  A 1 0 


(  A  1  2  ) 


where 


is  the  r^tu^y  viscous  damping  fact 


Substitution  of  equation  AM  into  A  1  2  gives 
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P  t  t  a  z  ’ 
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the  pitch  damping  per  slipper,  C  is 


(  1  -  + 
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For  roll  damping,  consider  Figure  A  .  2  o  n  i  y  acting 


the  roll  plane.  So  instead  "f  1  there  will  be  1,  and 


For  small  angular  rotations 


replaces  1  and  ■  •  becomes  -  < . 
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the  associated  spring  forces  are  given  by 
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(  A  1  5  ) 


where  -  k  ^  +  k  ^  and  Kr  =  +  k  ^  .  By  repeating  the 
process  shown  for  pitch  damping,  the  roll  damping  factor 
can  be  derived  as 
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( 1  ,  2k  + 

r  1  1 


1  “K  ) 
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(A  16  ) 


Bo  the  roll  damping  coefficient  for  each  slipper  can  be 
given  as 


1  “  K  ) 
r  r 


(  A  1  7  ) 
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T  h  e  s  e 
equivalent 
s  y  m  m  e  t  r  i  c  . 


t  h  o e  e  damping  terms  represent  the  total 
damping  at  a  slipper  provided  the  sled  is 
In  other  words. 


d  (  i  >  -  C  b  ♦  d  r  ♦  C  i  =  1  ,  A  (  A  1  8  ) 


However,  since  a  sled  is  not  totally  symmetric,  the 
!  .i  ;  p  ,  r  g  .1  *  .  Si  ipper  will  be  assumed  to  be  proportional 

to  trie  -•illness  associated  with  that  slipper  compared  to 


t  < 


i  i  ■  : 


lor 


i  n  y  given  slipper  its  damping 


Fallowing  the  procedure  for  bounce,  the  equation  reduces 
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or  t  h  c  lateral  damping  per  slipper  is 
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For  y ju  damping  consider  Figure  A. 4 
in o t  l  o  n  for  phis  system  is 
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i  ^  (  K  .  +  K  k  )  ] 
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which  eventually  results  in  an  e  q  u  a  t  i  o 
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.  The  equation  of 


=  0  (  A  2  3  ) 

for  damping  of 
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inis  can  now  do  expressed  as  a  damping  per  slipper 


long  as  the  directional  distinction  in  spring  stiffness 
is  considered.  Positive  directional  damping  will  he 
assumed  to  act  in  the  direction  of  the  positive 
coordinate  system  and  negative  directional  damping  in  the 
opposite  sense. 

From  track  experience,  it  has  been  found  that  the 
recommended  value  for  critical  damping,  f  .  ,  should  be 


APPENDIX  B 


DERIVATION  OF  STEADY  STATE  SLIPPER  FORCES — COUPLED  SLEDS 


B  .  1  COUPLED  SLEDS  DEFINITION 


Since  sleds  are  coupled  in  the  downtrack  condition 
only  (Figure  B.1),  velocity  and  acceleration  of  the  sled 
train  are  considered  to  be  constant.  However,  due  to  the 
spring  between  sleds  (which  represents  the  downtrack 
coupling  stiffness)  the  sleds  may  have  relative  velocity 
and  accelerations  with  respect  to  each  other.  Therefore, 


sled  one  (1)  will 
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to  simplify 

matters  it 
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Figure  B.l  COUPLED  SLED  SCHEMATIC 


However 


if  the  coupling  is  assumed  to  be  massless 


[  the  iperii.ii  force  terns  associated  with  the  coupling  :na 

be  ignored.  Therefore  the  transmitted  thrust  on  th 
front  of  one  sled  is  equal  in  magnitude  t  o  the 
|  transmitted  thrust  on  the  aft  of  the  sled  in  front  of  ft 
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Now  that  all  forces  acting  oP  a  sled  are  known,  th. 
slipper  reactions  can  be  solved  for.  A  :  i  r  ■  ■  x  a  n  i  n  .  n 
Figure  B  .  2  for  sled  i  and  summing  moments  ..bout  the  lr,-r 
slipper  beam  give s 
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It  should  be  noted  that  the  h^irntal 
components  of  the  transmitted  thrust  vector 
vertical  slipper  1  o  0  <j  5  I  u  e  t  ~  being  above 
Examining  Figure  B.j  for  the  dotted  lines, 
forCe,  the  vertical  reactions  can  be  f^uno 
moments  about  the  left  slipper,  or 
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G  U  mm  i  n  £  f  °  r  c  e  s  laterally,  >  F  :  G  ,  gives 

aLU  ,  1)  r  -GL  (  i  ,  2  )  -  HTT (  1  ,  1  )  -  HTT (  i  ,2)  (  B26  ) 

Again,  G L ( l  ,  1)  and  o  L (  i  , 2)  represent  two  slippers. 
There  is  no  easy  way  to  separate  these  into  the  two 
individual  components.  So  it  will  be  assumed  that  each 
component  will  receive  one  half  of  the  lateral  load. 

Therefore  for  sled  one,  using  program  nomenclature 
for  steady -st ate  forces,  3SF  becomes: 
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b  .  C : )  U  P  L L N G  r  U H C E  a  K H I  V  A T I  UN 

Ail  that  now  remains  is  to  define  the  components  of 
t  r .  < •  coupling  force  that  are  acting  o  p  the  sleds.  by 
examining  Figure  B.l,  it  is  seen  that  a  term  ELCPL  is 
i  n  c  u  d  <•  d  .  T  his  roprosonts  an  extra  1  i  n  k  needed  to  define 
the  mo  no  bull  type  coupling.  if  the  coupling  type  is  a 
i  '  i  r  m  a  r  (  ,  a  n  i  v  e“  r  r.  a  1  /irt)  ^  a  push  pad,  t  L  C  P  L  becomes 


B.j.1  MO NO BALL  COUPLING 


B.3.1.1  SLIPPER  GAP  GEOMETRY 


By  examining  Figure  B  .  V  it  ear  i-  e  -  •  -  r  •.  •  , 
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;ln  (  - mf(  n 


(  B  J  8  ) 


By  similar  reasoning  the  maximum  1  a  t  a  e  r a  1  rotation 
of  the  sled  may  be  given  as 


(  i  ,  1  ) 


-  1  ,  •  U-  (  i  )  +  I.A(  i  ) 

in  (  - no ( i  )  - 


(  B  2  9  ) 


Assuming  no  rigid  body  translation,  the  displacement 
of  the  coupling  connection  will  be  tahe  length  times  the 
sine  of  the  enclosed  <j  n  g  1  e  .  However,  the  length  now 


becomes  the  sum  of  the  length  between  slipper  beams  plus 


the  length  from  the  nearest  slipper  beam  to  the  coupling, 
o  r 


(  i  ,  j  )  -  (  D_D  (JLJ  +  DC  P  L  (  i  ,  k  )  )  s  i  n  ( 1  (  i  ,  j  )  )  (  B  3  0  ) 

> 

where  i  =  sled  ,  j  =  lateral  or  vertical  direction,  and 
R:frnnt  of  aft  slipper  beam  distance. 

Now  that  the  relative  displacements  of  each  end  of 
the  coupling  column  have  been  established,  it  is  possible 
to  find  the  angle  described  by  these  displacements.  The 
relataive  end  displacements  may  be  given  as 


A(  e  ,  j  )  =  i  ,  J  )  -  i  +  1  ,  j  ) 


(  B  3  1  ) 


Division  of  the  relative  displacements  by  the 
coupling  length  then  describes  the  enclosed  coupling 
angle,  or 


( c  ,  j  ) 


-  1 

s  i  n 


■'.(  C  ,  j  )_  , 

F.I.C  IM.  (  i  ) 


(  B  3  2  ) 


B.3.1.2  MONOBALL  FORCE  ANALYSIS 

Now  that  all  geometric  angles  of  the  sled  in  free 
flight  have  been  defined,  it  still  remains  to  be  shown 


how  they  affect  the  thrust  that  is  t r a n sm i t t a b 1 e  to  the 


coupling.  Examining  Figure  B.5  it  can  be  seen  that  the 
thrust  vector  from  the  rear  of  the  sled  may  be  resolved 
into  a  vector  that  is  solely  in  the  downtrack  direction 


(X-axis  o p 1 y ) ,  and  a 


tangent  to  the  X-axis.  So  we 


T  '  x  =  T  1  /  (  c  o  s  (  ;  (  1  ,  L  )  )  c  O  S  (  :  (  1  ,  V  )  )  ) 


i  m  1  1  a  r  1  y  , 


(  33  3  ) 


'  (  T  I  )  “  +  (  T  1  ) 

v  z 


,v))  +  T  i  "a  j  n  ~~  (  ;  (1  ,  I.  )  ) 

c  o  s  “  (  v  (  1  ,  v  )  ) 


( B34  ) 


s  in"  (  ‘  (  1  ,  v  )  )  4  s  in  2  (  1  (  1  , 1. )  ) 


■os  (  ;  (  I  ,  v  )  ) 


In  o  r  d  e  r  to  put  the  dynamic  system  into  equilibrium, 
we  will  nave  to  resolve  the  actual  transmitted  thrust 
available  to  input  into  sled  one.  This  is  accomplished 
by  reducing  the  total  system  thrust  by  the  inertial  and 
drjg  forces  experienced  by  the  coupling  and  the  rear 
pusher  sled.  Typically  these  components  on  the  sled 
train  are  insignificant  compared  to  the  thrust  magnitude, 
therefore  these  parameters  will  be  ignored  (as  mentioned 
in  lection  B.l).  In  other  woods,  all  thrust  available  at 


any  given  time  from  the  pusher  sled  is  directly 


And  this  must  react  T1  for  equilibrium. 

Following  this  logic  for  the  forebody  sled,  the 
thrust  vector  applied  to  it.  can  be  found  to  be 


T  2  :  F  C 

x  x 


=  T 1 / ( c  O  S  (  I  (  1  ,  V ) ) c  C  S  U (  1  ,  L ) )  )  (  B  41  ) 


Hence 


T2 


T2x  x  c  o  s  (  :  (  2  ,  L  )  )  c  o  s  (  .‘(2,  V)) 


(  B42) 


=  T  1  ( c o s (  :(2,L))cos(  :(2,V)))/(cos(:(l,L))cos(;(1,V))) 


B.J.1  COLEMAN  COUPLING  AND  PUSH  PAD  ANALYSIS 

By  examining  Figure  B.4,  a  Coleman  Coupling  or  Push 
Pad,  can  be  envisioned  by  letting  the  coupling  length 
(ELCPL)  be  zero.  By  following  the  same  derivation  as  for 
the  Mono  ball  Coupling,  the  Coleman  Coupling  or  the  Push 
Pad  can  be  found.  Since  the  coupling  length  is  zero,  it 
cannot  be  divided  into  the  relative  coupling  displacement 
(equation  B  j  2  )  ,  and  this  forces  the  defined  coupling 


2hS 


angles  to  be  zero. 


Et.N  COUPLING  DISPLACEMENT  TERMS 


In  jL'cti'-n  B.3.1.1,  coupling  displacement  terms  wee 
used  to  derive  the  sled  geometry  for  defining  coupling 
forces.  These  coupling  displacements  are  founc  by 
solving  Newton's  Second  Law  (F  =  ma)  on  the  sled  for  the 
rigid  body  acceleration  terms.  The  acceleration  is  then 
integrated  into  a  velocity  term  and  once  again  into 
displacements.  This  displacement  term  is  then  what  is 
used  for  the  coupling  geometry  analysis.  It  remains  t 
define  the  interrelation  of  the  sled's  position  on 
solving  for  this  coupling  displacement. 

h. n . 1  MONGBALL  COUPLING 

r or  the  monopjii  coupling,  the  jawbone  link  of  the 
coupling  allows  for  independent  motion  of  either  sled 
a  pout  t  n  e  link.  Therefore,  each  sled  can  be  solved 
1  n  d  r  p  e  n  d  e  n  t  1  y  of  e  ,,  c  h  other. 

B.R.S  COLEMAN  COUPLING 

Due  to  tin'  nature  of  Coleman  Couplings,  each  joint 
is  required  to  displace  an  amount  equal  to  its  mating 
point,  on  the  .idjucent  sled.  ..nee  the  lateral  and 
vertical  coupling  farces  were  included  in  the  analysis  to 


solve  fir  the  acceleration  terms,  the  displacements  r'f 


the  joints  will  be  identical.  Therefore  t  each 
this  type  of  coupling  can  be  solved  independently 


sled  of 
of  each 


o  t  h  e  r  . 


B.4.J  PUSH  TAD 


Since  this  coupling  type  allows  for  relative  motion 
with  respect  to  each  other  (through  slippage  across  the 
push  pad  surfaces),  the  integration  of  the  acceleration 
terms  can  be  done  independently  of  sled  interaction. 
However,  the  tangential  direction  coupling  forces  may 
have  to  be  modified  on  the  forebody  sleds,  if  slippage 
occurs.  Slippage  will  ~ c  e  u  r  when  the  tangential  coupling 
force  exceeds  the  static  frictional  resistance  force. 
This  static  frictional  force  is  given  by 


F  3 ( i , j )  :  TT (  1  ,  j  )  x  •  (  i  ) 


(bUj) 


where  ■  is  t h  e  c  o e f f i e 1 e  n  t  of  static  friction  for  * h e 
mated  in.itcnjls  '  f  the  push  pads.  Typic.il!  y  for  steel- 
o  n - s  t  e  e !  contact  this  is  0 .  j  . 

When  slippage  does  °  c  c  u  r  ,  the  tangential  forces 
transmitted  forward  will  be  equal  t  °  •„  h  e  :  f  a  '  ;  • 

frictional  force,  and  the  equations  of  mot  l  o  p  s  h  o  u  1  d  be 
solved  a  r  u  o  u  n  d  this. 


;’7d 
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Figure  C.17  Model  Dual  Rail,  S15000  Top 
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Figure  C.18  Modal  Dual  Rail,  S15000  Bottom 
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Figure  C.24  Modal  Dual  Rail,  44000  Bottom 
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Figure  C.28  Modal  Dual  Rail,  N15000  Bottom 
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Figure  C.47  Rigid  Body  Monorail,  44000  Top 
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Figure  C.48  Rigid  Body  Monorail,  44000  Bottom 
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Figure  C.53  Modal  Monorail,  TS  34000  Top 
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Figure  C.54  Modal  Monorail,  TS  34000  Bottom 
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